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20.  abstract  ]n  a  comparison  of\eleseismic  P-wave  delays  in  the  vicinity  of  the  Central 
Nevada  Test  Site  in  Hot  Creek  Valley.,  Nevada,  with  P-delay  data  over  a  wider 
region  in  the  Great  Basin,  we  found  that  upper  mantle  speeds  under  Hot  Creek 
Valley  stations  are  higher  than  the  average  for  the  Great  Basin  as  a  whole,  but 
lower  than  those  beneath  Pahute  Mesa.  These  observations  indicate  that  the  cal¬ 
dera  complex  in  Hot  Creek  Valley  may  have  a  high  speed  root  similar  to  that 
proposed  to  exist  beneath  the  Silent  Canyon  caldera  at  Pahute  Mesa.  In  con¬ 
trast,  the  Hot  Creek  Valley  anomaly  is  not  as  strong  as  the  Pahute  Mesa  ano¬ 
maly.  The  shadow  zone  caused  by  the  Pahute  Mesa  structure  is  much  more  pro¬ 
nounced  and  consequently  magnitudes  of  Pahute  Mesa  explosions  can  be 
underestimated  relative  to  the  magnitude  of  Hot  Creek  V^laj^explo^oi^s 
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y  Spectral  amplitudes  for  24  events  of  the  Mammoth  Lakes  earthquakes 
sequence  have  been  determined  for  the  frequency  range  0.1-10.0  Hz.  including 
the  ]Mi  6  earthquake  at  1450  UT  on  May  27.  1980,  which  is  of  particular  interest 
J  as  a  result  of  the  controversy  surrounding  its  mechanism.  "We  have  found  noth¬ 
ing  in  the  spectra  of  this  event  nor  ..in  the  spectra  of  the'aftershocks  to  distin- 
’  *  guish  them  from  spectra  of  "tectonic”  earthquakes.  However,  the  spectra  them¬ 


selves  do  not  distinguish  between  various  possible  explanations  for  the  non¬ 
double-couple  source  mechanism  observed  in  moment  tensor  inversion  at.d  first 
motion  data  for  the  largest  events  of  the  Mammoth  Lakes  earthquake  sequence. 

The  attenuation  of  Great  Basin  crustal  phases  in  the  Wood-Anderson 
bandpass  is  examined.  The  shape  of  the  logA0  (A)  curve  which  results  suggests 
more  attenuation  of  energy  in  the  Wood-Anderson  passband  than  in  southern 
California.  For  Great  Basin  paths  we  find,  for  events  with  Mi  £5.5,  that 
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where  A  is  epicentral  distance.  R  =  v'Az+/i2  is  hypocent.ral  distance,  and  h  is  focal 
depth.  Our  data  also  suggest  a  magnitude  bias  in  the  correction  term.  Using 
the  revised  magnitude  scale  and  seismic  moments,  M„,  estimated  from  spectral 
analysis,  the  data  are  well  fit  by  the  straight  line 

log  M„  =  (1.20±0.05)Afi  +  (17.492:0. 19) 


for  1  SMi  £6. 

The  relative  excitation  of  low  frequency  (0.05  to  0.125  Hz)  and  high  fre¬ 
quency  (0.5  to  10.0  Hz)  energy  is  compared  for  several  moderate  earthquakes 
(Mi  4  -  5.5)  in  the  1980  Mammoth  Lakes,  California  sequence.  Comparison  of  mb 
and  M,  for  the  Mammoth  Lakes  earthquake  sequence  and  Nevada  Test  Site 
explosions  show  that  all  of  the  earthquakes  analyzed  discriminate  from  explo¬ 
sions.  Seismic  moment  (AO  is  related  to  Ms  by 


log  M„  =  M,  +  19.4(±0.06) 


We  have  empirically  related  the  gross  logarithmic  amplitude  parameters  Mi , 
To;,,  and  AO  and  the  seismic  spectrum  observed  at  near  epicentral  distance. 

We  have  determined  an  average,  frequency  dependent,  apparent  Lg  Q  (/ ) 
function  for  Great  Basin  paths,  by  measuring  the  spatial  decay  of  spectral  ampli¬ 
tudes.  Spectra  of  the  Lg  phase  were  inverted  to  simultaneously  recover  the 
source  spectrum  and  frequency  dependent  Q.  Averaging  the  results  for  all  the 
paths  yields 

<?(/)  =  212/053  (0.2tfz  £/  £4.0tfz). 

This  result  is  in  excellent  agreement  with  Singh  and  Herrmann's  Q  model  for 
this  same  area,  based  on  observations  of  coda  decay. 

Observations  of  mb  bias  at  selected  sites  within  the  Soviet  Union  are  looked 
at.  Magnitude  mb  is  determined  for  five  earthquakes  on  25  and  27  May  1980, 
from  recordings  at  eight  Soviet  seismic  stations  on  a  4300  km-long  profile  from 
eastern  Kazakh  to  eastern  Siberia.  A  comparison  of  tabulated  magnitude 
ffsiduais  for  Soviet  seismic  stations  from  previous  work  of  Ringdal  (1985),  North 
(1976)  and  Vanek  et  al.  (1978,  1980)  shows  excellent  agreement  between  these 


studies. 


the  U  S  Geological  Survey  in  the  vicinity  of  the  proposed  Nuclear  Waste  Reposi¬ 
tory  at  Yucca  Mountain  in  southern  Nevada.  Gravity  and  well  log  data  are  used  to 
correct  the  refraction  travel-time  data  for  delays  due  to  lateral  variations  in 
near-surface  structure.  The  travel-times  for  the  two  profiles  were  combined, 
reparameterized  in  terms  of  delay  time  and  ray  parameter,  and  inverted  for 
average,  extermal  velocity-depth  bounds  for  the  region.  These  show  a  steep  gra- 
shallow^  crust  (depth  less  than  5  km)  with  speeds  increasing  from 
fe4.oA;s  to  6.2  ks  and  a  shallow  gradient  below  5  km  where  soeeds  increase 
from  *6.2 ks  1  to  6.5  Ars"1  at  25  km. 
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TECHNICAL  SUMMARY 


Keith  Priestley  and  David  Chavez 


Here  we  summarize  the  seven  projects  completed  during  the  time  period 
October  1.  1982  and  June  30,  1985.  Detailed  discussion  of  each  project  is 
covered  in  the  appendicies. 

Appendix  A  discusses  a  comparison  of  teleseismic  P-wave  delays  in  the 
vicinity  of  the  Central  Nevada  Test  Site  in  Hot  Creek  Valley,  Nevada,  with  P-delay 
data  over  a  wider  region  in  the  Great  Basin.  This  shows  that  upper  mantle  speeds 
under  Hot  Creek  Valley  stations  are  higher  than  the  average  for  the  Great  Basin 
as  a  whole,  but  lower  than  those  beneath  Pahute  Mesa.  These  observations  indi¬ 
cate  that  the  caldera  complex  in  Hot  Creek  Valley  may  have  a  high  speed  root 
similar  to  that  proposed  to  exist  beneath  the  Silent  Canyon  caldera  at  Pahute 
Mesa,  however  the  Hot  Creek  Valley  anomaly  is  not  as  strong  as  the  Pahute  Mesa 
anomaly.  These  high-speed  upper  mantle  structures  cause  defocusing  for  outgo¬ 
ing  rays  resulting  in  a  shadow  zone  at  the  distance  range  corresponding  to  many 
European  stations.  The  shadow  zone  caused  by  the  Pahute  Mesa  structure  is 
much  more  pronounced  and  consequently  magnitudes  of  Pahute  Mesa  explosions 
can  be  underestimated  relative  to  the  magnitude  of  Hot  Creek  Valley  explosions 
of  similar  yield.  When  magnitude  versus  yield  relationships  developed  using 
Pahute  Mesa  data  are  applied  to  other  areas,  such  as  Hot  Creek  Valley,  they  can 
result  in  a  biased  estimate  of  yield. 

In  Appendix  B,  spectral  amplitudes  for  24  events  of  the  Mammoth  Lakes 
earthquakes  sequence  have  been  determined  for  the  frequency  range  0.1-10.0 
Hz.  The  Mi  6  earthquake  at  1450  UT  on  May  27,  1980  is  of  particular  interest 
because  of  the  controversy  surrounding  its  mechanism.  A  comparison  of  spec¬ 
tral  levels  determined  from  analysis  of  locally  recorded  strong  motion  and 
broadband  digital  data  and  spectral  levels  from  regionally  recorded  surface 
waves  extrapolated  back  to  the  source  yields  consistent  results,  indicating  flat 
spectra  in  the  band  0. 1-1.0  Hz.  The  spectra  observed  in  this  study  do  not  show 
pronounced  spectral  peaks  predicted  by  theoretical  studies  of  ground  motion 
due  to  the  jerky  extension  of  a  fluid-driven  tensile  crack  containing  a  low- 
viscosity  fluid.  We  have  found  nothing  in  the  spectra  of  the  Mi  6  event  that 
occurred  at  1450  UT,  on  May  27.  1980,  nor  in  the  spectra  of  the  aftershocks  to 
distinguish  them  from  spectra  of  "tectonic"  earthquakes.  However,  the  spectra 
themselves  do  not  distinguish  between  various  possible  explanation  for  the  non¬ 
double-couple  source  mechanism  observed  in  moment  tensor  inversion  and  first 
motion  data  for  the  largest  events  of  the  Mammoth  Lakes  earthquake  sequence. 

In  Appendix  C.  we  examine  attenuation  of  Great  Basin  crustal  phases  in  the 
Wood-Anderson  bandpass.  The  local  magnitude  distance  correction  term, 
logA0  (A).  has  been  estimated  for  paths  within  the  Great  Basin  by  converting 
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seismograms  of  earthquakes  in  the  vicinity  of  Mammoth  Lakes,  California  into 
equivalent  Wood-Anderson  seismograms  and  measuring  the  decay  of  peak  ampli¬ 
tude  with  distance.  The  shape  of  the  logi40  (A)  curve  which  results  suggests 
more  attenuation  of  energy  in  the  Wood-Anderson  passband  than  in  southern 
California.  For  Great  Basin  paths  we  find,  for  events  with  Mi  £5.5,  that 
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where  A  is  epicentral  distance,  R  =  Va2+/i2  is  hypocentral  distance,  and  h  is  focal 
depth.  Our  data  also  suggest  a  magnitude  bias  in  the  correction  term,  in  that 
for  events  above  5.5,  peak  amplitudes  recorded  at  epicentral  distances 
below  20  km  decay  less  rapidly  than  those  for  smaller  events  recorded  at  the 
same  sites.  This  observation  is  consistent  with  the  hypothesis  that  there  is  a 
near  source  saturation  of  Mi  due  to  the  effects  of  a  finite  source  size  for  large 
earthquakes.  Using  the  revised  magnitude  scale  and  seismic  moments,  M0, 
estimated  from  spectral  analysis,  the  data  are  well  fit  by  the  straight  line 


log  M0  =  (l.20±0.05)A/i  +  (17.49±0. 19) 


for  l^Mi  £6.  M0  versus  ML  values  from  model  calculations  which  assume  con¬ 
stant  stress  drop  predict  curved  moment-magnitude  plots  over  this  magnitude 
range.  The  fact  that  we  observe  a  linear  relationrhip  suggests  that  the  stress 
drop  of  these  events  is  not  constant  but  rather  increases  with  Mi,  in  particular 
for  events  above  A/&«5.5. 

Appendix  D  compares  the  relative  excitation  of  low  frequency  (0.05  to  0.125 
Hz)  and  high  frequency  (0.5  to  10.0  Hz)  energy  for  several  moderate  earth¬ 
quakes  (. Mi  4  -  5.5)  in  the  I960  Mammoth  Lakes,  California  sequence.  Such  stu¬ 
dies  have  importance  for  both  source  mechanism  and  seismic  decrimination 
studies.  We  have  determined  local  magnitude  {Mi),  teleseismic  body  wave  mag¬ 
nitude  (mft),  and  surface  wave  magnitude  (Ms)  for  a  large  number  of  events  in 
the  magnitude  range  1  to  6.2.  Comparison  of  long-period  spectral  amplitude,  fl0, 
derived  from  body  waves  recorded  at  near  epicentral  distances  with  that  derived 
from  surfaces  waves  recorded  at  regional  distances  are  not  significantly 
different,  indicating  that  the  low  frequency  and  high  frequency  radiation  are 
consistent  with  a  simple  source.  Comparison  of  mb  and  Ms  for  the  Mammoth 
Lakes  earthquake  sequence  and  Nevada  Test  Site  explosions  show  that  all  of  the 
earthquakes  analyzed  discriminate  from  explosions.  Seismic  moment  ( M0 )  is 
related  to  Ms  by 
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logM0  =  Ma  +  19.4(±0.06) 

We  have  empirically  related  the  gross  logarithmic  amplitude  parameters  Mi, 
mb,  and  Ms,  and  the  seismic  spectrum  observed  at  near  epicentral  distance. 
One  representative  event  of  this  group  was  selected  as  a  standard  and  Mi,  mb 
and  Ms  for  that  event  were  scaled  (at  their  observed  frequencies)  to  the 
observed  seismic  spectrum.  Corresponding  magnitudes  for  three  other  events 
were  then  scaled  to  the  seismic  spectrum  using  the  derived  relationships.  For 
events  of  similar  depth,  the  local  spectral  values  can  be  predicted  from  the  dis¬ 
tant  magnitude  measurements  to  within  10%.  For  a  deeper  event,  spectral 
values  were  predicted  to  within  about  25%.  This  confirms  that  the  spectral 
shape  inferred  from  distant  measurements  of  magnitude  is  congruent  with  the 
spectral  shape  observed  directly  in  the  near  field.  The  majority  of  Mammoth 
Lakes  events  examined  can  be  interpretated  as  simple  events  while,  in  com¬ 
parison,  several  of  the  San  Fernando  aftershocks  in  the  same  magnitude  range 
studied  by  Tucker  and  Brune  (1977)  could  be  interpreted  as  partial  stress  drop 
events.  The  near  equality  of  high  and  low  frequency  radiation  for  Mammoth 
Lakes  aftershocks  suggests  that,  for  each  event,  the  seismic  radiation  over  the 
frequency  band  studied  is  coming  from  a  simple  source,  possibly  the  breaking  of 
a  single  asperity.  The  differences  suggested  between  Mammoth  Lakes  aft¬ 
ershocks  and  the  San  Fernando  earthquakes  may  arise  because  most  of  the 
events  studied  for  Mammoth  Lakes  are  occurring  in  the  more  competent  (ie. 
higher  static  stress),  granitic  block  above  a  deeper,  through-going  fault  plane. 
The  data  suggest  that  partial  stress  drop  events  are  more  common  on  developed 
fault  systems,  whereas  full  stress  drop  events  are  more  common  in  the  sur¬ 
rounding  country  rock  as  it  adjust  to  the  stress  field  generated  by  movement  on 
the  developed  fault  plane. 

In  appendix  E  we  have  determined  an  average,  frequency  dependent, 
apparent  Lg  £?(/)  function  for  Creat  Basin  paths,  by  measuring  the  spatial  decay 
of  spectral  amplitudes.  We  examined  seismograms  for  several  earthquakes 
occurring  around  the  periphery  of  the  Great  Basin  recorded  by  wide-band  digital 
seismographs.  Spectra  of  the  Lg  phase  were  inverted  to  simultaneously  recover 
the  source  spectrum  and  frequency  dependent  Q.  Averaging  the  results  for  all 
the  paths  yields 


<?(/)  =  212/' 


(0.2  Hz  z  f  S4.0  Hz). 


This  result  is  in  excellent  agreement  with  Singh  and  Herrmann's  Q  model  for 
this  same  area,  based  on  observations  of  coda  decay. 

Appendix  F  discusses  observations  of  mb  bias  at  selected  sites  within  the 
Soviet  Union.  Magnitude  is  determined  for  five  earthquakes  on  25  and  27 
May  1980,  from  recordings  at  eight  Soviet  seismic  stations  on  a  4300  km-long 
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profile  from  eastern  Kazakh  to  eastern  Siberia.  The  records  were  hand- 
digitized,  and  magnitudes  were  determined  from  traces  corrected  for  instru¬ 
ment  response  as  well  as  the  uncorrected  traces.  Four  of  the  five  earthquakes 
were  at  Mammoth  Lakes,  California,  in  the  western  Great  Basin;  the  fifth  was  at 
Tonga  in  the  south  Pacific.  Magnitude  residuals  with  respect  to  network- 
averaged  m6's  listed  in  the  ISC  Bulletin  were  high  for  stations  at  Yakutsk  and 
Seymchan  in  Siberia,  low  for  raypaths  through  the  Baikal  rift  zone,  and  inter¬ 
mediate  for  a  station  at  Semipalatinsk,  about  100  km  from  the  East  Kazakh  test 
site.  A  comparison  of  tabulated  magnitude  residuals  for  Soviet  seismic  stations 
from  previous  work  of  Ringdal  (1985),  North  (1978)  and  Vanek  et  al.  (1978,  I960) 
shows  excellent  agreement  between  these  studies.  Our  results  are  slightly  more 
scattered  but  in  good  agreement  with  previous  work.  Recalculation  of  for  83 
events  recorded  on  granite  at  the  Nevada  Test  Site  provided  a  determination  of 
magnitude  bias  ( Jm*  =  -0.10  +-  0.35)  for  the  NTS  granite  site  with  respect  to  ISC 
magnitude.  Comparison  of  this  value  with  Jm6's  for  the  Soviet  stations  provides 
a  direct  measure  of  the  magnitude  bias  of  the  NTS  area  relative  to  areas  in  the 
USSR  in  which  the  stations  are  located.  This  bias  reflects  only  differences  in 
attenuation,  and  does  not  account  for  other  effects  such  as  differences  in  cou¬ 
pling,  focussing,  and  tectonic  release.  The  magnitude  bias  due  to  attenuation 
differences  between  the  NTS  granite  site  and  the  station  at  Semipalatinsk  is 
-0.20. 

Finally,  in  Appendix  G,  we  interpret  data  from  two  seismic  refraction 
profiles  recorded  by  the  U.S.  Geological  Survey  in  the  vicinity  of  the  proposed 
Nuclear  Waste  Repository  at  Yucca  Mountain  in  southern  Nevada.  The  first 
profile  consists  of  33  recordings  of  a  nuclear  explosion  at  the  Nevada  Test  Site 
and  extends  from  45  to  110  km.  The  second  consists  of  71  recordings  of  a  chem¬ 
ical  explosion  near  Beatty,  Nevada,  and  extends  from  0  to  65  km.  Gravity  and 
well  log  data  are  used  to  correct  the  refraction  travel-time  data  for  delays  due 
to  lateral  variations  in  near-surface  structure. 

The  travel-times  for  the  two  profiles  were  combined,  reparameterized  in 
terms  of  delay  time  and  ray  parameter,  and  inverted  for  average,  extermal 
velocity-depth  bounds  for  the  region.  These  show  a  steep  gradient  in  the  shallow 
crust  (depth  less  than  5  km)  with  speeds  increasing  from  Rs4.5fcs-1  to  8.2  ks~x 
and  a  shallow  gradient  below  5  km  where  speeds  increase  from  «36.2fcs~'  to  6.5 
ks~l  at  25  km. 

A  specific  velocity-depth  model  consistent  with  the  extremal  bounds  has 
been  found  by  two-dimensional  ray  tracing.  The  structure  of  the  upper  crust  is 
highly  variable  with  the  Paleozoic  basement  varying  between  the  surface  and  3.5 
km  depth.  A  thick,  low-speed  section  is  indicated  by  large  delays  across  the 
Crater  Flat-Prospector  Pass  caldera  complex.  A  prominent  arrival  1.5  to  2.0 


seconds  after  the  initial  P-wave  time,  observed  in  the  distance  range  40  to  75 
km,  is  interpreted  as  a  mid-crustal  boundary  at  12  to  16  km  depth.  The  profiles 
are  of  insufficient  length  to  see  Pn,  however  a  prominent  arrival  3  to  4  seconds 
after  Pgl  observed  at  distances  greater  than  80  km,  is  interpreted  as  PmP  from 
the  Moho  at  approximately  35  km  depth. 
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MAGNITUDE  BIAS  IN  THE  GREAT  BASIN  AND  ITS  IMPLICATIONS 
FOR  EXPLOSION  MAGNITUDE  VERSUS  YIELD  ESTIMATES 

Keith  F.  Priestley  and  David  E.  Chavez 

Seismological  Laboratory,  University  of  Nevada 


Abstract.  Comparison  of  teleseismic  P-wave  delays  in 
the  vicinity  of  the  Central  Nevada  Test  Site  in  Hot 
Creek  Valley,  Nevada,  with  P-delay  data  over  a  wider 
region  in  the  Great  Basin,  shows  that  upper  mantle 
speeds  under  Hot  Creek  Valley  stations  are  higher  than 
the  average  for  the  Great  Basin  as  a  whole,  but  lower 
than  those  beneath  Pahute  Mesa.  These  observations 
indicate  that  the  caldera  complex  in  Hot  Creek  Valley 
may  have  a  high  speed  root  similar  to  that  proposed  to 
exist  beneath  the  Silent  Canyon  caldera  at  Pahute 
Mesa,  however  the  Hot  Creek  Valley  anomaly  is  not  as 
strong  as  the  Pahute  Mesa  anomaly.  These  high-speed 
upper  mantle  structures  cause  defocusing  for  outgoing 
rays  resulting  in  a  shadow  zone  at  the  distance  range 
corresponding  to  many  European  stations.  The  shadow 
zone  caused  by  the  Pahute  Mesa  structure  is  much 
more  pronounced  and  consequently  magnitudes  of 
Pahute  Mesa  explosions  can  be  underestimated  relative 
to  the  magnitude  of  Hot  Creek  Valley  explosions  of 
similar  yield.  When  magnitude  versus  yield  relation¬ 
ships  developed  using  Pahute  Mesa  data  are  applied  to 
other  areas,  such  as  Hot  Creek  Valley,  they  can  result 
in  a  biased  estimate  of  yield. 

Introduction 

The  acceptance  of  a  Threshold  Test  Ban  Treaty 
requires  the  ability  to  accurately  estimate  the  yield  of 
nuclear  explosions.  Such  estimates  are  based  primarily 
on  seismic  measurements,  and  the  bulk  of  the  data  for 
relating  seismic  amplitudes  (or  magnitudes)  to  yield 
comes  from  the  Nevada  Test  Site  (NTS).  The  near 
surface  structure  at  NTS  is  complicated  by  extensive 
Cenozoic  volcanism,  and  it  is  possible  that  this  has 
affected  structures  as  deep  as  the  upper  mantle  low 
velocity  zone  beneath  NTS.  Thus,  it  is  important  to 
evaluate  how  well  the  empirical  magnitude  versus  yield 
curves  derived  for  NTS  can  predict  the  yield  of  explo¬ 
sions  at  other  sites. 

The  underground  nuclear  explosion  Faultless  was 
detonated  at  the  Central  Nevada  Test  Site  (CNTS)  in 
Hot  Creek  Valley,  Nevada,  approximately  150  km 
north  of  NTS.  The  International  Seismological  Centre 
(ISC)  body  wave  magnitude,  mb ,  of  6.3  is  more  than 
0.3  units  (a  factor  of  two  in  amplitude)  greater  than 
that  predicted  from  empirical  magnitude-yield  curves 
for  NTS  explosions  (Everden,  1970],  given  the  max- 
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imum  yield  of  1000  kilotons  [Springer  and  Kinnaman, 
1971],  CNTS  is  physiographically  similar  to  Pahute 
Mesa  at  NTS,  and  in  view  of  the  close  proximity  of  the 
two  sites  it  might  be  expected  that  the  magnitude-yield 
curves  for  NTS  sould  apply  to  CNTS.  Here  we  provide 
evidence  that  this  is  not  the  case. 

Spence  [1974)  used  travel-time  residuals  for  NTS 
explosions  recorded  at  teleseismic  distances  to  infer  the 
presence  of  a  high  speed  structure  in  the  mantle 
beneath  the  14  m.y.  old  Silent  Canyon  caldera  at 
Pahute  Mesa.  Hedge  and  Noble  [1971]  argue  that 
Silent  Canyon,  and  in  fact  most  of  the  basalts  of  the 
southern  Great  Basin  region,  originate  in  the  upper 
mantle.  Thus,  the  high  speed  material  proposed  by 
Spence  [1974]  may  represent  upper  mantle  material 
largely  depleted  of  its  partial  melt  component  by 
differentiation  and  associated  eruption  of  the  partial 
melt  phase.  Differences  in  upper  mantle  structure 
between  Hot  Creek  Valley  and  Pahute  Mi  t  may  be 
the  cause  of  the  Faultless  magnitude  anomaly. 

Inferences  on  Upper  Mantle  Structure 

Chavez  and  Priestley  [in  preparation]  have  studied 
teleseismic  P-delays  for  the  Great  Basin.  Among  the 
data  used  in  this  study  are  recordings  from  seven 
seismographs  operated  by  the  U.  S.  Geological  Survey 
in  and  around  Hot  Creek  Valley  (Fig.  1).  Here,  we 
summarize  the  results  for  the  Hot  Creek  Valley  data. 

Figure  1  shows  the  generalized  geology  in  the  vicinity 
of  Hot  Creek  Valley.  The  area  is  dominated  by  a 
series  of  nested  calderas  which  range  in  age  from  25  to 
35  m.y.,  i.e.,  10  to  20  m.y.  older  than  Silent  Canyon 
caldera  at  Pahute  Mesa.  To  the  southeast  of  Hot 
Creek  Valley  is  the  Lunar  Craters  volcanic  field,  a  site 
of  Quaternary  volcanism  [Stewart,  1980]. 

We  first  searched  for  lateral  variations  in  structure 
under  Hot  Creek  Valley  by  computing  teleseismic  P- 
wave  residuals  at  each  station  relative  to  the  average 
residual  across  the  CNTS  network.  Figure  2  shows 
relative  residual  as  a  function  of  azimuth.  Station  loca¬ 
tions  are  indicated  by  circles  centered  at  each  site.  The 
circles  have  a  radius  of  1.0  second,  and  the  circumfer¬ 
ence  represents  a  zero  residual.  Lines  drawn  from  the 
circle  are  scaled  according  to  the  size  of  the  residual. 
Those  inside  the  circle  are  negative  residuals  (early 
arrivals)  and  lines  outside  the  circle  are  positive  residu¬ 
als  (late  arrivals).  The  azimuth  of  each  line 
corresponds  to  the  azimuth  of  approach  of  the  wave- 
front.  Figure  2a  shows  the  residuals  at  each  of  the 
CNTS  stations,  after  the  data  have  been  corrected  for 
known  variations  in  sediment  thickness  and  for 
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Fig.  1.  Generalized  geologic  map  of  the  Hot  Creek  Valley,  Nevada  area  including  station  loca¬ 
tions  and  the  Faultless  shot  point. 
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inferred  Moho  dip.  The  data  show  that  the  residuals 
across  the  Hot  Creek  Valley  network  are  small,  indicat¬ 
ing  no  strong  small-scale  lateral  variations  in  structure. 

Residuals  relative  to  event  average  will  not  reveal 
any  large-scale  anomalies  under  Hot  Creek  Valley 


which  affect  the  entire  network.  We  therefore  com¬ 
puted  P-wave  residuals  relative  to  a  regional  event 
average  using  stations  throughout  the  Great  Basin. 
The  results,  shown  in  Figure  2b,  indicate  that  Hot 
Creek  Valley  is  relatively  fast  compared  to  the  Great 
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Fig.  2.  Maps  of  travel  time  residuals  versus  azimuth  (format  described  in  the  text):  a)  relative  to 
the  event  average  for  Hot  Creek  Valley  stations  after  corrections  for  surface  geology  and  inferred 
Moho  dip;  b)  relative  to  event  average  for  Great  Basin  seismic  stations;  c)  for  Faultless  and 
Greeley  relative  to  Dumont  at  Hot  Creek  Valley,  Pahute  Mesa,  and  Yucca  Flat,  respectively 
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Fig.  3.  a)  Ray  paths  through  a  hypothetical  high 
speed  body  beneath  Pahute  Mesa  as  postulated  by 
Spence  11974],  b)  Ray  paths  through  a  hypothetical 
high  speed  body  beneath  Hot  Creek  Valley  postu¬ 
lated  by  Chavez  and  Priestley  [in  preparation] 

Basin  as  a  whole,  suggesting  that,  as  beneath  Silent 
Canyon,  a  high  speed  vestige  of  a  root  to  the  Hot 
Creek  Valley  caldera  complex  exits  at  depth. 

In  Figure  2c  we  compare  travel  time  residuals  for 
Faultless  in  Hot  Creek  Valley,  and  Greeley  at  Pahute 
Mesa,  to  the  Dumont  explosion  at  Yucca  Rat.  Com¬ 
pared  to  Dumont,  speeds  under  Faultless  are  slow 
while  those  under  Greeley  are  fast.  This  suggests  that 
Yucca  Rat  travel  times  are  also  being  affected  by  the 
high  speed  material  beneath  Pahute  Mesa,  or  prehaps 
by  another,  similar  structure.  A  plausible  explanation 
for  these  observations  is  that  high  speed  mantle 
material  exists  under  both  Pahute  Mesa  [as  proposed 
by  Spence,  19741  as  well  as  Hot  Creek  Valley.  The  Hot 
Creek  Valley  caldera  complex  is  considerably  older 
than  that  at  Pahute  Mesa  however,  and  consequently 
the  anomalous  structure  under  Hot  Creek  Valley  has 
had  more  time  to  equilibrate  with  the  surrounding 
mantle.  As  a  result,  the  magnitude  of  the  Hot  Creek 
Valley  anomaly  is  less. 

Effects  of  Upper  Mantle  Structure 
on  Teleseismic  Amplitudes 

Figure  3  compares  ray  paths  from  hypothetical  explo¬ 
sions  above  the  proposed  upper  mantle  structures 


beneath  Pahute  Mesa  (Fig.  3a)  and  Hot  Creek  Valley 
(Fig.  3b).  The  models  shown  are  qualitative  and  are 
only  intended  to  illustrate  a  possible  mechanism  for  the 
amplitude  anomaly  observed  for  Faultless.  In  our  cal¬ 
culations,  the  high  speed  upper  mantle  bodies  shown 
by  the  stippled  regions  were  defined  by  a  gradation 
from  a  speed  of  7.7  km/ sec  in  the  "unperturbed'  upper 
mantle  to  8.2  km/sec  in  the  interior  of  the  high  speed 
zone.  The  location  of  the  high  speed  material  under 
Pahute  Mesa  is  taken  from  Spence  [1974],  and  that 
under  Hot  Creek  Valley  is  from  Chavez  and  Priestley 
[in  preparation].  The  density  of  rays  shown  in  Figure 
3  is  proportional  to  the  expected  amplitude.  In  both 
Figures  3a  and  3b,  there  is  a  smooth  variation  in  ray 
density  (and  hence  amplitude)  from  20°  to  approxi¬ 
mately  45°.  At  distances  greater  than  60°,  there  is  a 
substantial  defocusing  of  the  rays  leaving  Pahute  Mesa, 
resulting  in  a  shadow  zone.  The  shadow  zone  resulting 
from  the  structure  at  Hot  Creek  Valley  is  weak  in  com¬ 
parison.  Based  on  this  figure,  we  would  expect  that  at 
distances  beyond  60°,  the  amplitudes  for  events  in  Hot 
Creek  Valley  would  be  greater  than  those  for  similar 
sized  events  at  Pahute  Mesa. 

Most  mb  measurements  for  Nevada  explosions 
reported  by  the  ISC  are  heavily  biased  by  amplitude 
measurements  made  at  European  stations.  For  Pahute 
Mesa  explosions  (Fig.  3a),  these  stations  fall  in  the 
shadow  zone  while  Figure  3b  suggests  that  events 
detonated  in  Hot  Creek  Valley  will  be  affected  to  a 
much  smaller  extent.  Figure  4  is  a  comparison  of  the 
variation  in  amplitudes  of  Faultless  and  Greeley  with 
distance  along  a  north-east  azimuth  towards  Europe. 
The  values  plotted  are  ratios  of  the  Faultless  or  Gree¬ 
ley  P-wave  amplitude,  to  the  amplitude  of  Dumont  at 
Yucca  Rat.  (Spence  [1974]  used  Dumont  as  a  refer¬ 
ence  in  his  travel  time  study,  assuming  that  it  occurred 
above  more  normal  Great  Basin  mantle.  Our  P-delay 
studies  indicate  that,  compared  to  the  average  Great 
Basin,  Yucca  Rat  is  also  affected  by  high  speed 
material,  however  here  we  are  only  concerned  with 
differences  between  NTS  and  CNTS.)  The  amplitudes 
for  Faultless  and  Greeley  are  similar  at  distances  less 
than  40°.  Beyond  that,  the  Faultless  amplitude  ratios 
remain  close  to  a  constant  level  with  increasing  dis¬ 
tance,  while  those  for  Greeley  decrease.  At  distances 
between  60°  to  80°  (ie.,  European  stations).  Faultless 
clearly  has  larger  relative  amplitudes  than  Greeley  ,  the 
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Fig.  4.  Comparison  of  Faultless,  Greeley  and 
Kasseri  P-wave  amplitudes  relative  to  Dumont  as  a 
function  of  distance  along  a  northeast  azimuth. 
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Greeley  amplitude  ratio  having  decreased  by  a  factor  of 
2  to  3. 

Discussion 

The  variation  in  amplitude  with  distance  shown  in 
Figure  4  supports  our  hypothesis  proposed  here  that 
low  amplitudes  observed  at  European  stations  result 
from  defocusing  of  waves  due  to  anomalous  upper 
mantle  structure  beneath  Pahute  Mesa.  Lynnes  and 
Lay  [1984]  have  also  noted  low  amplitudes  for  Pahute 
Mesa  explosions  along  northeast  azimuths.  However 
the  amplitude  variation  shown  in  Figure  4  for  Greeley 
could  be  the  result  of  tectonic  release.  Wallace,  et  al. 
[1983]  studied  wave  forms  for  a  large  number  of 
explosions  at  Pahute  Mesa  and  found  that  a  number  of 
the  early,  large  explosions  (and  Greeley  in  particular) 
show  waveform  distortion  which  they  attribute  to  a 
significant  component  of  tectonic  release.  In  com¬ 
parison  to  the  1966  Greeley  explosion,  Wallace,  et  al. 
[1983]  found  that  the  1975  Kasseri  explosion  located  2 
km  from  Greeley  but  with  similar  yield,  showed  simple 
waveforms  indicating  only  a  small  component  of  tec¬ 
tonic  release.  In  Figure  4  we  also  compare  the  short- 
period  P-wave  amplitudes  of  Greeley  and  Kasseri  at 
common  stations  along  a  northeast  azimuth  and  find 
that  the  amplitude  versus  distance  relationship  for 
Kasseri  is  identic  '  to  that  shown  for  Greeley,  indicat¬ 
ing  that  the  Greeiey  amplitude  variation  is  not  solely 
the  result  of  tectonic  release. 

The  yield  of  nuclear  explosions  can  also  be  estimated 
from  surface  waves.  Lg  waves  are  a  superposition  of 
higher  mode  crustal  surface  waves,  and  as  such  they 
are  free  from  the  effects  of  anomalous  upper  mantle 
structure  beneath  the  site  of  the  explosion.  We  have 
compared  Lg  amplitudes  for  Greeley  and  Faultless  at 
two  short  period  World  Wide  Standard  Seismograph 
Network  stations  in  the  eastern  United  States,  and  find 
no  significant  difference  in  the  Lg  amplitudes  between 
the  two  events.  Based  on  these  measurements.  Fault¬ 
less  was  no  larger  (and  possibly  somewhat  smaller) 
than  Greeley. 

These  results  have  important  implications  for  the 
estimation  of  yields  of  explosions  based  on  teleseismic 
body-wave  magnitudes.  The  majority  of  magnitude 
versus  yield  data  are  derived  from  explosions  at  NTS, 
where  most  of  the  large  explosions  have  been 
detonated  on  Pahute  Mesa.  The  ISC  body  wave  mag¬ 
nitudes  for  these  events  are  biased  by  the  large 
number  of  stations  reporting  from  Europe.  If  our 


hypothesis  is  correct,  then  the  anomalous  upper  mantle 
structure  beneath  Pahute  Mesa  will  result  in  lower 
amplitudes  at  European  stations,  causing  the  magni¬ 
tude  versus  yield  curve  to  be  biased  towards  low  mag¬ 
nitudes  for  a  given  yield.  Application  of  the  magnitude 
versus  yield  curve  derived  from  NTS  explosions  to 
other  areas,  such  as  the  CNTS,  can  result  in  an  overes¬ 
timate  of  yield,  as  appears  to  have  been  the  case  for 
Faultless. 
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Spectral  amplitudes  for  24  events  of  the  Mammoth  Lakes  earthquake  sequence  have  been  determined 
for  the  frequency  range  0.1-10,0  Hz.  The  ML  6  earthquake  that  occurred  at  1450  UT  on  May  27,  1980,  is 
of  particular  interest  because  of  the  controversy  surrounding  its  mechanism.  A  comparison  of  spectral 
levels  determined  from  analysis  of  locally  recorded  strong  motion  and  broadband  digital  data  and 
spectral  levels  from  regionally  recorded  surface  waves  extrapolated  back  to  the  source  yields  consistent 
results,  indicating  flat  spectra  in  the  band  0.1- 1.0  Hz.  The  spectra  observed  in  this  study  do  not  show 
pronounced  spectral  peaks  predicted  by  theoretical  studies  of  ground  motion  due  to  the  jerky  extension 
of  a  fluid-driven  tensile  crack  containing  a  low-viscosity  fluid.  We  have  found  nothing  in  the  spectra  of 
the  Mt  6  event  that  occurred  at  1450  UT,  on  May  27,  1980,  nor  in  the  spectra  of  the  aftershocks  to 
distinguish  them  from  spectra  of  “tectonic”  earthquakes.  However,  the  spectra  themselves  do  not  dis¬ 
tinguish  between  various  possible  explanations  for  the  non-double-couple  source  mechanism  observed  in 
moment  tensor  inversion  and  first  motion  data  for  the  largest  events  of  the  Mammoth  Lakes  earthquake 


sequence. 

Introduction 

The  Mammoth  Lakes  earthquakes  are  located  at  the  inter¬ 
section  of  the  Sierra  Nevada  frontal  fault  system  and  the  Long 
Valley  caldera.  The  caldera  was  formed  0.7  m.y.  ago  by  col¬ 
lapse  and  subsidence  associated  with  the  eruption  of  the 
Bishop  tuff  [Bailey  et  al.,  1976],  and  volcanism  has  continued 
on  a  reduced  scale  into  the  Holocene.  The  dominant  fault  in 
the  area  south  of  the  caldera  is  the  Hilton  Creek  fault  along 
which  several  hundred  meters  of  pure  normal  fault  displace¬ 
ment  has  occurred  since  the  formation  of  the  caldera  [ Bailey 
et  al.,  1976],  Extension  of  the  Hilton  Creek  fault  into  the 
caldera  appears  to  have  occurred  as  recently  as  0.3  m.y.  ago, 
possibly  indicating  that  the  caldera  had  cooled  sufficiently  by 
that  time  to  support  stresses  large  enough  to  generate  earth¬ 
quakes. 

The  Mammoth  Lakes  earthquake  sequence  clearly  began 
with  an  WL  5.8  event  on  October  4,  1978,  located  beneath  the 
Wheeler  Crest,  30  km  northwest  of  Bishop,  California.  Fol¬ 
lowing  this  event,  activity  gradually  spread  to  the  north  and 
west.  The  main  energy  release,  which  began  on  May  25,  1980, 
was  preceded  for  over  a  year  by  ML  3.5-45  earthquake 
swarms  :n  the  region  of  the  M L  6  events.  These  swarms  gener¬ 
ally  clustered  near  the  southern  boundary  of  the  caldera  and 
migrated  eastward  toward  the  eventual  epicenter  of  the  first 
VfL  6  event  [Cramer  and  Toppozada,  1980].  The  main  episode 
of  strain  release  commenced  with  an  ML  6.1  (all  local  mag¬ 
nitude  estimates  for  the  main  events  are  from  Uhrhammer  and 
Ferguson  [1980])  event  at  1633  UT  on  May  25  (hereafter 
referred  to  as  May  25,  1980 — 1633)  and  was  followed  by  a 
second  ,Vf,  6.0  event  at  1649  UT.  These  initial  events  occurred 
along  an  east-west  trend  near  the  southern  boundary  of  the 
caldera.  Seismicity  then  migrated  southward  from  the  caldera 
boundary.  Cramer  and  Toppozada  [1980]  note  that  no  events 
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at  the  magnitude  5  level  occurred  south  of  the  southern 
boundary  of  the  caldera  until  after  the  ML  6.1  event  at  1944 
UT,  which  occurred  6-7  km  south  of  the  first  ML  6  event.  The 
last  of  the  large  events  (\fL  6  2)  occurred  about  9  km  south  of 
the  1944  event,  at  1450  UT  on  May  27  (hereafter  referred  to  as 
May  27,  1980—1450)  Lide  and  Ryall  [1984]  have  found  that 
in  the  2  months  following  the  main  events,  aftershocks  and 
continuing  swarms  occurred  along  the  southern  boundary  of 
the  caldera  and  along  several  NNE-SSW  trends  nearly  or¬ 
thogonal  to  the  caldera  boundary.  The  first  and  second  ML  6 
events  occurred  along  the  first  of  these  trends,  and  the  first, 
third,  and  last  of  the  ML  6  events  occurred  along  the  latter 
trend.  The  epicenters  of  nearly  all  of  the  aftershocks  and  all  of 
the  M L  6  events  arc  well  to  the  west  of  the  projection  of  the 
Hilton  Creek  fault  at  depth,  indicating  that  this  is  not  the 
causative  fault  for  the  earthquake  sequence. 

Considerable  controversy  has  arisen  regarding  the  source 
mechanisms  of  the  larger  events  of  the  Mammoth  Lakes  se¬ 
quence.  Focal  mechanisms  derived  from  local  and  regional 
first  motion  data  indicate  strike-slip  motion  with  a  NE-SW  T 
axis  and  NW-SE  P  axis  [Cramer  and  Toppozada.  1980;  Ryall 
and  Ryall,  1980],  This  implies  right-lateral  motion  on  east- 
west  trending  faults  or  left-lateral  motion  on  north-south 
trending  faults.  Unconstrained  moment  tensor  inversion  of 
long-period  teleseismic  body  and  surface  wave  data  result  in  a 
large,  non-double-couple  component  for  the  first  and  last  of 
the  M,  6  events,  and  the  inversion  of  the  teleseismic  waveform 
data  plus  teleseismic  first  motion  data  indicate  mechanisms  for 
the  first,  third,  and  fourth  of  the  M L  6  events  which  differ 
significantly  from  strike  slip  [Given  et  al.,  1982;  Barker  and 
Langston.  1983], 

The  observed  non-doublc-couple  mechanism  may  arise  in  a 
number  of  ways.  Given  et  al.  [1982]  have  constrained  the 
moment  tensor  to  be  a  double  couple  and  interpreted  the 
results  in  terms  of  left-lateral  oblique  slip  on  NNE  trending 
planes  dipping  eastward.  They  suggest  that  the  discrepancy 
between  the  regional  and  teleseismic  mechanisms  could  arise 
from  either  structural  effects  which  distort  the  teleseismically 
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Fig.  1.  Location  maps  showing  the  events  studied  and  the  stations  from  which  seismograms  were  analyzed.  In  the 
lower  right-hand  map  the  solid  circles  denote  locations  of  events,  triangles  denote  location  of  strong  motion  accelero- 
graphs.  and  squares  denote  location  of  digital  seismographs.  The  upper  left-hand  map  shows  the  locations  of  the  long- 
period  ..smographs  used  for  the  surface  wave  analysis. 


observed  radiation  pattern  or  complex  rupture  during  the 
event.  Distortion  of  the  outgoing  rays  may  result  from  com¬ 
plex  structure  beneath  the  Long  Valley  caldera,  for  example 
the  presence  of  magma  chambers,  which  were  first  proposed 
by  Bailey  et  al.  [1976]  using  geological  evidence  and  Hill 
[1976]  using  seismological  evidence.  Similar  anomalies  in 
focal  mechanisms  have  been  reported  for  mid-oceanic  ridge 
earthquakes,  and  these  have  been  attributed  to  defocusing  of 
the  ray  paths  due  to  magma  beneath  the  rise  crest  [Solomon 
and  Julian,  1974],  IVallace  [1984]  and  Lide  and  Ryall  [1984] 
have  interpreted  the  non-double-couple  mechanism  as  arising 
from  a  multiple  source.  Reversal  of  polarity  observed  on 
short-period  and  long-period  instruments  at  the  same  site 
[IVallace  et  al.,  1982]  suggest  a  complex  source  time  function 
with  the  short-period  arrivals  representing  the  failure  of  an 
asperity  (possibly  with  a  strike-slip  mechanism)  and  the  long- 
period  arrivals  due  to  the  overall  response  of  the  faulting  epi¬ 
sode.  Similar  non-double-couple  mechanisms  for  events  in 
Greece  have  been  shown  to  be  the  result  of  events  with  rup¬ 
ture  on  multiple  planes  (D.  MacKenzie,  personal  communi¬ 
cations,  1984). 

On  the  other  hand,  Julian  [1983]  has  suggested  that  the 
results  of  inverting  the  surface  wave  and  long-period  body 
wave  data  and  taking  the  regional  and  teleseismic  first  motion 
data  as  a  whole  are  better  fit  by  a  compensated  linear  vector 
dipole  (CLVD)  mechanism  [Knopoff  and  Randall,  1970]  than 
they  are  by  a  double-couple  mechanism.  Julian  has  noted  that 
the  CLVD  mechanism  is  consistent  with  all  but  four  and  two 
of  the  first-motion  data,  respectively,  for  the  May  25,  1980 — 
1633  and  May  27,  1980 — 1450  M L  6  events,  as  opposed  to  32 
and  16  inconsistent  data  for  the  double-couple  mechanism, 
and  has  suggested  that  the  May  25,  1980 — 1633  and  May  27, 
1980 — 1450  events  are  due  to  injection  of  a  fluid  into  an  ex¬ 
panding  crack,  i.e.,  dike  formation.  Aki  [1984],  utilizing  the 
calculations  of  Chouet  [1981],  has  suggested  that  the  non¬ 


double-couple  mechanism  results  from  the  rapid  opening  of  a 
channel  between  two  preexisting  cracks  filled  with  magma. 

In  this  study  we  compare  spectral  amplitudes  of  24  Mam¬ 
moth  Lakes  earthquakes  ranging  in  magnitude  from  ML  3.5  to 
M L  6.2,  over  the  frequency  band  *0.05-20  Hz.  Spectral  am¬ 
plitudes  from  approximately  0.5  to  20  Hz  are  derived  from 
spectral  analysis  of  locally  recorded  body  wave  data.  Spectral 
amplitudes  from  approximately  0.05  to  0.125  Hz  are  derived 
from  analysis  of  regionally  recorded  surface  wave  data.  The 
wide  bandwidth  of  the  spectral  amplitude  measurements 
allows  us  to  compare  the  observed  spectral  shapes  with  theo¬ 
retical  spectral  shapes  for  ground  motion  resulting  from  “tec¬ 
tonic”  earthquakes  and  ground  motion  predicted  for  a  propa¬ 
gating  tensile  crack. 

Data 

This  study  is  based  on  seismograms  of  the  Mammoth  Lakes 
earthquakes  recorded  on  the  seismic  stations  shown  in  Figure 
1.  The  events,  shown  as  solid  circles,  comprise  only  a  small 
fraction  of  the  total  number  of  events  which  have  occurred  in 
the  Mammoth  Lakes  area  since  October  1978.  Locations  after 
May  27,  1980,  are  largely  from  C.  S.  Lide  (written  communi¬ 
cation,  1984).  Locations  for  some  of  the  earlier  events  are  from 
the  University  of  Nevada  or  University  of  California.  Berkeley 
bulletins.  The  events  selected  for  analysis  were  those  recorded 
on  at  least  one  local  strong  motion  or  digital  seismograph  and 
one  long-period  regional  setsmograph.  The  SMA-1  strong 
motion  accelerographs,  shown  as  triangles  in  Figure  1,  were 
installed  by  the  California  Division  of  Mines  and  Geology 
[Turpin,  1980]  or  by  the  University  of  Southern  California 
(USC)  [Moslem  et  al.,  1983].  We  also  used  shear  wave  spectra 
supplied  by  R.  Archuleta  from  data  recorded  by  the  U.S.  Geo¬ 
logical  Survey.  Digital  stations,  shown  by  squares  in  Figure  1 
and  described  by  Archuleta  et  al.  [1982],  consist  of  either  force 
balance  accelerometer  or  velocity  transducers  and  Sprengne- 
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ther  DR-100  recorders.  For  events  along  the  southern  bound-  [ Thatcher  and  Hanks ,  1973].  At  some  azimuths  the  amplitude 

ary  of  the  caldera  the  SMA-1  and  digital  seismographs  cover  a  could  be  1.67  times  larger  than  the  amplitude  determined  as- 

large  fraction  of  the  focal  sphere  compared  with  many  pre-  suming  this  correction,  while  near  a  node  in  the  radiation 

vious  studies  of  a  similar  nature.  The  radiation  pattern  of  pattern  the  amplitude  could  be  quite  small.  Furthermore,  if 

events  south  of  the  caldera  are  not  as  well  sampled.  the  source  dimension  is  comparable  to  the  wavelength  in- 

The  inset  map  in  Figure  1  shows  the  locations  of  the  re-  volved,  source  propagation  might  distort  the  radiation  pattern 

gjonal  long-period  stations  from  which  seismograms  were  by  focusing  or  defocusing  energy,  thereby  changing  the  appar- 

used  to  determine  surface  wave  magnitude.  These  stations  ent  corner  frequency  for  different  directions  and  possibly  re- 

completely  surround  the  Mammoth  Lakes  region  and  are  all  ducing  the  calculated  moment. 

World-Wide  Standard  Seismograph  Network  (WWSSN)  sta-  No  Q  correction  was  made  in  determining  the  spectra.  Ar- 
tions,  with  the  exception  of  the  California  Institute  of  Tech-  chuleta  et  al.  [1982]  found  that  whole  path  Q  had  a  small 
nology  (CIT)  station  at  Pasadena  (PAS).  effect  on  the  measurement  of  ft0  for  Mammoth  Lakes  events. 


Analysis 

Determination  of  Body  Wave  Moment 

Displacement  amplitude  spectra  were  computed  by  Fourier 
transforming  the  S  wave  arrival  on  the  strong  motion  accel¬ 
erograms  and  then  dividing  the  acceleration  spectra  by  co'  2  to 
give  the  displacement  spectra.  The  time  window  was  chosen  to 
be  long  enough  to  include  the  entire  S  wave  arrival  but  short 
enough  to  minimize  the  inclusion  of  surface  wave  energy.  The 
spectra  have  been  normalized  to  a  distance  of  10  km.  Since  the 
precise  locations  of  the  foci  of  the  events  are  not  known,  we 
have  not  used  the  epicenter  and  seismograph  coordinates  to 
determine  the  distance  but  have  estimated  the  epicentral  dis¬ 
tance  R  from  the  S-P  times  assuming  an  average  S  wave  ve¬ 
locity  of  3.4  km  s' 1  and  an  average  P  wave  velocity  of  5.8  km 
s'  *.  In  addition,  we  have  not  rotated  the  data  into  transverse 
and  radial  components  but  have  computed  the  spectra  for 
both  horizontal  components  and  computed  the  vector  sum  of 
these  two  values. 

The  spectra  were  interpreted  in  terms  of  a  flat,  low- 
frequency  asymptote  fl0.  a  high-frequency  cu~2  falloff,  and  the 
corner  frequency  f{  at  which  the  two  trends  intersect.  The 
far-field  spectral  parameters  were  then  related  to  the  source 
parameter  seismic  moment  M0  using  the  relationship  given  by 
Brune  [1970]. 

v.  -  m 

where  p  is  the  density  ip  =  2.8  g  cm'3),  /?  is  the  shear  wave 
velocity  l/f  =  3.5  km  s'1).  R&,  is  the  rms  average  of  the  radi¬ 
ation  pattern  (taken  here  as  0.6),  and  k  is  a  correction  factor 
for  the  amplification  of  SH  waves  at  the  free  surface  (taken  as 
2)  The  correction  of  the  amplitude  of  the  S  wave  on  being 
reflected  at  the  free  surface  assume-  that  the  waves  are  hori¬ 
zontally  polarized  (SH).  Since  we  have  not  rotated  the  data 
into  transverse  and  radial  components,  the  data  are  in  general 
some  mixture  of  SH  and  ST  motion.  Abe  [1974]  discusses 
experimental  studies  of  amplification  on  reflection  and  con¬ 
cludes  that  using  k  =  2  applies  approximately  to  SH  wave 
signals  containing  some  Sk'  motion.  Helmberger  and  Malone 
[1975]  discuss  problems  due  to  multiple  arrivals  caused  by 
reflections  in  deducing  source  parameters  from  S  waves  re¬ 
corded  at  epicentral  distances  much  greater  than  the  source 
depth.  The  station-source  configuration  shown  in  Figure  1 
indicates  that  this  should  not  be  a  problem  in  this  study. 

A  source  of  error  inherent  in  the  estimate  of  and  conse¬ 
quently  in  the  estimates  of  M0  comes  from  the  assumption  of 
a  single,  average  correction  for  radiation  pattern.  If  sufficient 
instruments  were  available  or  if  the  source  parameters  includ¬ 
ing  direction  and  velocity  of  rupture  were  known,  this  source 
or  error  could  be  reduced.  A  value  of  0  6  was  used  for  R u 


Determination  of  Surface  Wave  Moment 

Tucker  and  Brune  [1977]  determined  the  seismic  moment  of 
seven  aftershocks  of  the  San  Fernando  earthquake  from  the 
surface  wave  spectra.  This  gave  an  estimate  of  seismic  moment 
independent  of  that  obtained  from  S  wave  spectra  and  one 
based  on  lower  frequencies  (0.03-0.10  Hz  in  their  case).  They 
made  a  careful  estimate  of  the  errors  involved  in  estimating 
the  seismic  moment  from  the  surface  wave  spectra.  The  errors 
included  procedural  errors,  path  structure  and  attenuation, 
lateral  refraction  (multipathing),  source  finiteness,  and  uncer¬ 
tainties  in  fault  parameters  (depth,  dip,  strike,  and  slip  vector) 
in  correcting  for  radiation  pattern.  They  concluded  that  these 
effects  could  lead  to  a  maximum  error  in  the  final  result  of  no 
more  than  a  factor  of  4. 

In  this  study  we  also  determined  the  seismic  moment  from 
surface  waves  for  comparison  with  seismic  moments  deter¬ 
mined  from  the  S  wave  spectra.  However,  we  used  a  different 
procedure  which  allowed  us  to  determine  the  seismic  moment 
from  a  large  number  of  events  without  the  time-consuming 
digitization  of  the  surface  wave  seismograms. 

Surface  wave  magnitude  (Ms) defined  by  the  formula  [ Mar¬ 
shall  and  Basham,  1972] 

Ms  =  log  A  +  B'( A)  +  P(T)  (2) 

was  determined  from  seismograms  of  long-period  seismograph 
stations,  primarily  in  the  western  United  States.  In  equation 
(2),  A  is  one  half  the  peak-to-peak  ground  amplitude  of  the 
maximum  amplitude  Rayleigh  wave  with  period  T  at  distance 
A.  B'(A)  corrects  for  average  attenuation,  scattering,  geo¬ 
metrical  spreading,  and  refraction,  and  P(T )  is  a  correction 
factor  for  dispersion  and  allows  measurement  of  the  surface 
wave  amplitude  at  any  period.  Marshall  and  Basham  [1972] 
have  tabulated  P{T)  for  continental  North  America  for 
periods  ranging  from  10  to  40  s.  In  some  cases  the  maximum 
surface  wave  amplitude  observed  for  the  Mammoth  Lakes 
earthquakes  was  in  the  8-  to  1 0-s  range.  Since  P(T)  is  a 
smoothly  varying  function,  we  have  extrapolated  P(T)  to  8-s 
period.  The  surface  wave  magnitudes  for  the  Mammoth  Lakes 
events  studied  here  are  tabulated  in  Table  1. 

The  surface  wave  moment  is  linearly  related  to  the  surface 
wave  amplitude  and  thus  to  the  surface  wave  magnitude.  We 
have  determined  a  relationship  between  the  surface  wave  mag¬ 
nitude  and  the  seismic  moment  in  the  following  manner.  Sur¬ 
face  wave  synthetics  were  generated  by  a  computer  routine 
(written  by  G.  Masters)  for  an  earthquake  with  M0  =  1021 
dyn  cm  at  a  number  of  regional  seismographs  from  which 
data  were  used  to  determine  the  surface  wave  magnitudes. 
From  the  synthetic  surface  wave  seismograms  we  determined 
the  surface  wave  magnitude  in  the  same  manner  as  discussed 
above  for  the  Mammoth  Lakes  earthquakes  and  thus  deter- 
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TABLE  1 .  Summary  of  Earthquake  Source  Parameters  and  Comparison  of  Short-Period  and 
Long-Period  Estimates  of  Seismic  Moment 


Focal  Parameters 


Body  Wave  Surface  Wave 


Mq,  M  Q, 


Date 

Time 

Latitude 

Longitude 

Depth 

Ml 

dyn  cm 

EM0 

M, 

dyn  cm 

EM0 

Dec.  8,  1979 

2138 

37.58 

118.88 

6.9 

4.2 

1.94  x  10“ 

2.23 

4.2 

X 

10“ 

1.66 

May  17.  1980 

0114 

37.62 

118.80 

1.3 

4.3 

5.18  x  102‘ 

3.01 

26  x  10“ 

May  26,  1980 

0436 

37.58 

118.80 

0.8 

4.2 

201  x  10“ 

283 

1.7 

X 

!0“ 

262 

May  26,  1980 

1020 

37.60 

118.78 

4.6 

4.5 

9.77  x  1011 

3.25 

4.5 

X 

10“ 

1.24 

May  26,  1980 

1224 

37.53 

118  87 

114 

5.1 

2.68  x  10“ 

1.22 

4.18 

3.8 

X 

10“ 

1.37 

May  26,  1980 

1857 

37.56 

1 18.83 

0.2 

5.7 

4.93  x  10“ 

1.17 

4.35 

5.6 

X 

10“ 

1.33 

May  27,  1980 

0201 

37.50 

118.90 

6.0 

3.6 

1.13  x  10“ 

231 

5.1 

X 

10“ 

May  27,  1980 

1450 

37.45 

118.79 

11.7 

6.2 

7.24  x  10“ 

1.81 

5.45 

7.2 

X 

10“ 

1.18 
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X 
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1 18.83 
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mined  a  relationship  between  the  surface  wave  magnitude  and 
seismic  moment. 

Factors  affecting  the  calculation  of  the  synthetic  surface 
wave  seismograms  are  the  earth  velocity  and  attenuation 
structure  and  the  fault  parameters  of  the  individual  earth¬ 
quakes.  Since  the  events  are  located  at  the  western  edge  of  the 
Great  Basin,  while  most  of  the  seismograph  stations  are  lo¬ 
cated  in  western  North  America,  we  have  used  the  velocity 
structure  determined  for  the  Great  Basin  from  surface  wave 
dispersion  measurements  [ Priestley  and  Brune,  1978],  The  at¬ 
tenuation  structure  for  the  Great  Basin  by  Patton  and  Taylor 
[1984]  was  used,  with  a  slight  lowering  of  Qt  from  that  of 
Patton  and  Taylor  in  the  surficial  layer.  Since  we  do  not  have 
fault  plane  solutions  nor  do  we  know  the  focal  depth  for  the 
majority  of  events  we  have  studied,  we  determined  in  average 
relationship  for  the  Mammoth  Lakes  area  in  the  following 
manner.  Vetter  and  Ryall  [1983]  have  determined  fault  plane 
solutions  for  a  large  number  of  events  in  the  Mammoth  Lakes 
area  and  found  that  the  mechanisms  can  generally  be  divided 
into  two  groups.  The  first  group  consists  of  earthquakes  with 
strike-slip  mechanisms  (average  plane  of  strike  N14°E,  dip 
80  ESE,  slip  N16  E  or  strike  N73°W,  dip  80°SSW,  slip 
N76'E)  which  are  most  commonly  at  depths  less  than  9  km. 
The  second  group  consists  of  earthquakes  with  oblique  or 
normal  mechanisms  (average  planes  of  strike  N03°E,  dip  64°E, 
slip  N40’E  or  strike  N50"W,  dip  38’SW,  and  slip  N86°W) 
usually  at  depths  greater  than  9  km  Synthetic  surface  wave 
seismograms  were  calculated  for  both  mechanisms  for  2,  5,  8, 
11.  and  14  km  depth.  Surface  wave  magnitudes  were  deter¬ 
mined  from  the  10  records  at  each  station  and  an  average 
conversion  factor  determined  for  both  focal  mechanisms  and 
all  focal  depths.  The  resulting  relationship  between  surface 
wave  magnitude  Ms  and  seismic  moment  M0  for  the  Mam¬ 
moth  Lakes  area  is  given  by 

log  M0  =  Ms  +  19.40  ( ±0.06)  (3) 


To  verify  the  validity  of  our  moment-surface  wave  mag¬ 
nitude  relationship,  we  can  compare  equation  (3)  with  other 
moment  magnitude  relations  established  for  the  western 
United  States.  Wyss  and  Brune  [1968]  originally  established  a 
moment-local  magnitude  relationship  for  events  less  than 
magnitude  6  in  the  western  United  States.  Seismic  moments 
were  determined  for  272  events  which  averaged  over  a  wide 
range  of  tectonic  regions  (California,  Nevada,  Arizona,  Utah, 
and  Baja  California).  These  were  found  to  be  best  fit  by  the 
curve  given  by 

log  M0  =  1.7 Ml  +  15.1  3  <  Ml  <  6 

Using  the  Ms  -  Mt  relationship  given  by  tVyss  and  Brune 
[1968],  this  becomes  in  terms  of  M^G)  (the  surface  wave  mag¬ 
nitude  defined  by  Gutenberg  [1945]): 

log  M0  =  Mj(C)  +  19.2  3  <  Ms  <  6 

This  definition  cam  be  compared  with  the  relationship  derived 
directly  from  Gutenberg's  [1945]  definition  of  surface  wave 
magnitude.  According  to  that  definition,  a  magnitude  6  earth¬ 
quake  produces  a  far-field  displacement  of  100  gm  at  distance 
22“  for  surface  waves  of  period  20  s.  In  terms  of  moment 
Gutenburg's  relationship  becomes 

log  M0  =  MjtG)  +  19.3 

This  point  thus  represents  average  of  the  numerous  observa¬ 
tions  on  which  the  surface  wave  magnitude  was  based.  As 
pointed  out  by  Richter  [1958],  the  scale  was  adjusted  to  agree 
with  the  local  magnitude  values  of  6-7.  Tucker  and  Brune 
[1977]  did  not  explicitly  give  a  relationship  between  surface 
wave  magnitude  and  seismic  moment.  However,  from  the 
values  given  in  their  Tables  I  and  3.  (he  average  relationship  is 

log  Vf0  =  Mj  +  19  3 

This  relationship,  determined  using  a  Gutenburg  continental 
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Fig.  2.  Spectra  from  strong  motion  recordings  at  the  sites  “McGee  Creek  Inn“  and  “Long  Valley  Fire  Station"  for  the 
Ml  4.9  event  occurring  at  1S16  UT  on  May  31,  1980.  Spectra  are  the  vectorial  sum  of  the  two  orthogonal  components  and 
have  been  normalized  to  10  km.  The  long-period  spectral  levels  picked  from  these  spectra  are  denoted  by  the  arrow,  and 
the  long-period  spectral  levels  determined  from  the  surface  wave  observations  are  denoted  by  the  triangle  with  error  bars. 


earth  model  and  the  Marshall- Basham  definition  of  surface 
wave  magnitude,  is  in  good  agreement  with  the  relationship 
based  on  the  Gutenburg  definition  of  surface  wave  magnitude. 
Comparing  these  various  relationships  between  Ms  and  M0, 
we  conclude  that  the  use  of  the  Great  Basin  velocity  [ Priestley 
and  Brune,  1978]  and  attenuation  [Patton  and  Taylor,  1984] 
model  has  increased  the  estimated  moment  of  the  events  by 
26%  over  that  determined  for  the  Gutenburg  continental 
model. 

Equation  (3)  allows  us  to  rapidly  determine  the  seismic 
moment  M0  from  the  surface  wave  magnitude  observation  in 
Table  I.  The  resulting  values  of  seismic  moment  (see  Table  1) 
were  converted  to  an  equivalent  value  of  by  solving  equa¬ 
tion  (1)  for  fl0  and  substituting  the  seismic  moment  derived 
from  the  surface  wave  observations.  This  gives  values  of  Q„  at 
frequencies  of  the  surface  waves  for  comparison  with  the  spec¬ 
tra  obtained  from  the  body  wave  record  on  the  acc^Wome- 
ters. 

Error  Analysis 

Errors  in  Q0  (and  consequently  in  M0)  were  estimated  by 
calculating  an  average  spectral  value  and  a  standard  deviation 
where  more  than  one  recording  was  available  for  an  event. 
The  average  values  of  fl0  was  determined  by 

<n0>  =  antilog  £  log  [Oq  K,/10]| 


where  O0,  is  the  long-period  spectral  level  at  the  ith  station,  R, 
is  the  distance  to  the  ith  station,  and  NS  is  the  number  of 
stations  recording  the  event.  The  standard  deviation  of  the  log 
average  value  of  Q0  was  determined  by 

s.d.dog  <no»  -  I  [log  (n0)  -  log  <n0>]2| 

A  multiplicative  error  factor  is  then  determined  from 
Eflo  =  antilog  {s.d.(log  <Q„»} 

As  pointed  out  by  Archuleta  et  al.  [1982],  when  calculating 
average  values  of  Qq  or  its  derivative  M0,  it  is  necessary  to 
compute  the  averages  using  the  equations  in  the  above  forms 
to  give  equal  weight  to  each  observation.  If  a  simple  arithmeti¬ 
cal  average  value  is  determined,  it  will  be  biased  toward  larger 
values,  since  the  errors  associated  with  are  lognormally 
distributed. 

Observed  Spectra 

Figures  2  and  3  are  displacement  spectra  for  two  events  of 
the  24  Mammoth  Lakes  earthquakes  we  studied.  All  spectra 
are  computed  from  strong  motion  accelerograms  of  the  USC- 
CIT  data  set.  The  plotted  spectra  are  the  vectorial  average  of 
the  two  orthogonal  components  and  have  been  normalized  to 
a  hypocentral  distance  of  10  km.  The  arrows  indicate  our 
picks  of  the  corner  frequency  and  long-period  spectral  level. 
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Fig  3  Spectra  from  strong  motion  recordings  al  the  sites  “Cash  Baugh  Ranch"  and  “Long  Valley  Fire  Station"  for  the 
M[  6.2  event  occurring  at  1450  UT  on  May  27,  1980.  The  notation  is  the  same  as  for  Figure  2. 
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Fig.  4.  Companion  of  the  body  wave  and  surface  wave  deter¬ 
mined  moments  for  the  24  events  studied.  The  line  is  for  reference  and 
is  of  unit  slope  and  zero  intercept. 

The  triangle  with  error  bars  denotes  the  average  spectral  level 
determined  from  the  surface  wave  measurements. 

Figure  2  is  two  of  14  spectra  for  an  ML  4.9  event  occurring 
at  1516  UT  on  May  31,  1980  (hereafter  referred  to  as  May  31, 
1980 — 1516).  The  focus  of  this  event  is  nearly  identical  to  that 
of  the  May  25,  1980 — 1633  event.  The  average  value  of  the 
long-period  spectral  level  <fi0>  measured  from  the  14  spectra 
is  1.03  x  10' 1  cm  s,  corresponding  to  M0  =*  1.34  x  10“  dyn 
cm.  Our  values  include  estimates  of  flofrom  the  data  of  Archu¬ 
leta  et  al.  [1982]  (sis  observations)  plus  eight  additional  spec¬ 
tra  from  the  USC-CIT  strong  motion  data.  Our  value  of  <fl0> 
docs  not  differ  significantly  from  their.  The  seismic  moment 
from  the  surface  wave  observation  is  1.1  x  10“  dyn  cm  which 
is  not  significantly  different  from  the  body  wave  value. 

Spectra  of  the  largest  earthquake  we  studied,  the  May  27, 
1980 — 1450  Ml  6  shock,  are  shown  in  Figure  3.  This  event  is 
of  particular  interest  because  of  the  controversy  concerning  its 
mechanism.  Four  values  of  moment  have  been  published  for 
this  event.  Archuleta  et  al.  [1982]  determined  a  moment  of 
2.33  x  10“  dyn  cm  from  the  spectral  analysis  of  two  locally 
recorded  seismograms.  Uhrhammer  and  Ferguson  [1980]  esti¬ 
mated  a  moment  of  5.0  x  10“  dyn  cm  from  the  spectra  of 
51.2  s  of  the  broadband  displacement  record  of  the  vertical 
component  seismogram  at  Jamestown,  California,  140  km 
west  of  Mammoth  Lakes.  Barker  and  Langston  [1983]  invert¬ 
ed  seven  long-period  teleseismic  P  waveforms  and  five  long- 
period  teleseismic  SH  waveforms  to  obtain  the  moment  tensor 
for  this  event.  They  found  a  moment  of  1.03  x  10“  dyn  cm 
and  36%  CLVD.  Given  et  al.  [1982]  inverted  long-period  sur¬ 
face  wave  from  this  event  and  used  this  result  to  model  the 
long-penod  body  waves.  They  found  a  moment  of  1.10  x  10“ 
dyn  cm.  Our  values  of  7.26  x  10“  and  7.20  x  10“  dyn  cm  for 
the  body  and  surface  wave  observations,  respectively,  are 
somewhat  higher  than  the  values  reported  by  Archuleta  et  al. 
[1982]  but  in  reasonable  agreement  with  the  other  reported 
values  of  moment. 

Discussion 

Figure  4  is  a  plot  comparing  the  surface  waves  moments 
and  body  wave  moments  of  the  Mammoth  Lakes  earthquakes 
from  Table  I  The  reference  line  through  the  points  has  unit 


slope  and  intercept  zero.  Those  points  lying  above  the  line 
have  larger  body  wave  moments  compared  to  the  observed 
surface  wave  moment,  while  those  below  the  line  have  larger 
surface  wave  moments  compared  to  the  observed  body  wave 
moment.  Vertical  error  bars  denote  the  spread  in  the  observed 
body  wave  moment,  while  horizontal  error  bars  denote  the 
spread  in  the  observed  surface  wave  moment.  It  is  clear  from 
this  figure  that  there  are  no  significant  differences  between  the 
moments  determined  at  8-20  s  from  surface  waves  and  the 
moments  determined  at  periods  near  1-s  from  body  waves. 
Where  we  have  analyzed  events  in  common,  our  results  from 
the  surface  wave  measurements  support  the  moment  estimates 
of  Archuleta  et  al.  [1982]  and  thus  their  findings  of  nearly 
constant  stress  drop  (A a  ~  50  bars)  for  events  with  seismic 
moments  greater  than  approximately  1  x  10“  dyn  cm.  All  of 
the  events  we  have  studied  are  larger  than  1  x  10“  dyn  cm, 
and  therefore  we  cannot  substantiate  the  dependence  of  stress 
drop  on  seismic  moment  which  Archuleta  et  al.  [1982]  ob¬ 
served  for  smaller  events.  However,  the  wider  bandwidth  of 
our  spectra  permit  us  to  address  the  problem  of  the  source 
mechanism  controversy  surrounding  the  Mammoth  Lakes 
events. 

Aki  et  al.  [1977]  and  Chouet  [1981]  derive  theoretical 
ground  motion  in  the  far  and  near-feld,  respectively,  due  to 
the  extension  of  a  fluid-driven  tensile  crack  embedded  in  a 
layered  half-space.  In  both  studies  the  source  was  a  jerky 
opening  of  a  channel  ahead  of  the  crack  tip  due  to  excess 
pressure  of  the  fluid  in  the  crack.  Seismic  waves  are  generated 
by  the  vibrations  of  the  newly  created  crack  walls  due  to  the 
rapid  application  of  excess  fluid  pressure.  Both  the  far-ficld 
and  near-field  studies  show  a  peaked  character  to  the  ground 
motion  spectra,  in  contrast  to  “tectonic”  earthquake  spectra, 
which  are  flat  to  some  corner  frequency  and  then  fall  off  at 
higher  frequencies.  Chouet  [1981]  varied  the  model  parame¬ 
ters  to  analyze  the  effects  of  fluid  compressibility,  source 
depth,  and  structure  of  the  medium  on  the  ground  motion.  He 
found  that  the  spectral  peaks  depend  on  me  source  geometry, 
the  medium  characteristics,  the  receiver  position,  the  compo¬ 
nent  of  ground  motion  being  studied,  and  critically  on  the 
crack  stiffness  factor  C  defined  as  C  =  bL/dft,  where  b  is  the 
bulk  modules  of  the  fluid,  L  is  the  full  length  of  the  crack,  d  is 
the  crack  width,  and  n  is  the  rigidity  of  the  solid.  The  crack 
vibrates  most  easily  when  the  crack  is  empty  (C  =  0)  [see  Aki 
et  al.,  1977,  Figures  6-9],  The  spectral  peak  becomes  broad¬ 
ened,  however,  when  the  crack  stiffness  factor  is  about  5, 
implying  that  the  bulk  modulus  of  the  magma  is  low  enough 
to  contain  bubbles  [_Aki,  1984],  Chouet's  [1981]  model  calcu¬ 
lations  show  a  pattern  similar  to  the  non-double-couple  mech¬ 
anism  observed  for  the  largest  Mammoth  Lakes  earthquakes 
(May  25,  1980—1633.  May  25.  1980—1450).  Utilizing 
Clwuet's  [1981]  calculations.  Aki  [1984]  has  proposed  that 
these  events  result  from  the  sudden  opening  of  a  channel  be¬ 
tween  two  preexisting  fluid-filled  cracks  due  to  the  higher  fluid 
pressure  in  one  of  the  cracks.  The  pressure  in  the  magma 
drops  due  to  the  increased  total  crack  volume,  and  the  crack 
tip  closes  as  a  result. 

However,  as  shown  in  Figure  4.  none  of  the  events  we  have 
studied  show  significant  differences  in  level  in  the  frequency 
band  0.1  -I  0  Hz  and  therefore  do  not  show  the  spectral  peak¬ 
ing  predicted  by  the  Aki  [1984]  model  for  a  low-viscosity 
fluid.  Most  of  the  strong  motion  and  digital  recording  seismo¬ 
graphs  for  the  May  27,  1980 — 1450  event  were  located  to  the 
north  of  this  event  thus  sampling  only  a  portion  of  the  focal 
sphere,  hence  it  may  be  argued  that  we  have  sampled  that 
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part  of  the  radiation  pattern  where  the  peaked  nature  of  the 
spectrum  is  least  prominent.  Figure  S  shows  the  distribution  of 
spectra  for  the  ML  4.9  May  31,  1980 — 1516  aftershock  which 
occurred  in  almost  the  same  location  as  the  May  25,  1980 — 
1633  Ml  6  event  which  also  had  a  non-double-couple  mecha¬ 
nism  [Given  et  al.,  1982;  Julian,  1983].  The  focal  mechanism 
from  local  and  regional  first-motion  data  for  the  May  31, 
1980 — 1516  event  is  from  U.  R.  Vetter  (written  communi¬ 
cation,  1984)  and  is  very  similar  to  the  solution  independently 
determined  by  Archuleta  et  al.  [1982],  Figure  5  shows  an 
upper  hemisphere  plot  with  the  northwest  and  southeast  being 
the  compression  quadrants.  This  solution  is  similar  to  the 
focal  mechanism  given  by  Cramer  and  Toppozada  [1980]  and 
Ryall  and  Ryall  [1980]  from  local  and  regional  first  motion 
data  for  the  May  25,  1980 — 1633  event.  The  location  of  the 
stations  from  which  the  spectra  were  computed  are  indicated 
on  the  focal  sphere.  Some  of  the  individual  spectra  show  a 
suggestion  of  peaking  near  1  Hz  but  not  of  the  magnitude 
predicted  by  Aki  et  al.  [1977]  and  Chouet  [1981].  We  find 
nothing  in  the  spectra  of  the  Mammoth  Lakes  events  we  have 
studied  which  distinguish  them  from  spectra  of  “tectonic" 
events  observed  in  numerous  previous  studies.  The  flat  nature 
of  the  spectra  does  not  rule  out  the  opening  of  a  tensile  crack 
as  the  source  of  the  elastic  waves  for  these  earthquakes.  It 
does  preclude  the  opening  of  a  tensile  crack  with  the  subse¬ 
quent  injection  of  a  fluid  of  low  viscosity. 

Summary 

Spectra  for  24  Mammoth  Lakes  earthquakes  greater  than 
1  x  10J1  dyn  cm,  including  the  ML  6  event  occurring  at  1450 
UT  on  May  27,  1980,  have  been  determined.  The  short-period 
portion  of  the  spectra  (0.5-10.0  Hz)  have  been  determined 
from  spectral  analysis  of  the  locally  recorded  strong  motion 
and  broadband  digital  data.  A  relationship  has  been  derived 
between  the  Marshall-Basham  surface  wave  magnitude  and 
seismic  moment,  allowing  us  to  determine  the  seismic  moment 
directly  from  time  domain  amplitude  measurements  of  the 
surface  waves  recorded  at  regional  stations.  This  gives  an  esti¬ 
mate  of  the  seismic  moment  independent  of  that  obtained 
from  the  S  wave  spectra  and  one  based  on  lower  frequencies 
(0.05-0.125  Hz).  A  comparison  of  the  seismic  moments  indi¬ 
cate  that  for  any  one  event,  no  significant  difference  in  the 
spectral  level  exist  between  the  values  determined  at  surface 
wave  frequencies  (0.05-0.125  Hz)  and  the  S  wave  corner  fre¬ 
quency  (1-5  Hz).  The  events  studied,  and  in  particular  the 
May  27,  1980 — 1450  ML  6  event,  do  not  show  peaked  spectra 
predicted  from  theoretical  analysis  of  the  ground  motion  re¬ 
sulting  from  a  fluid-driven,  propagating  tensile  crack.  We  have 
found  nothing  in  the  spectra  of  these  events  to  distinguish 
them  from  "tectonic"  earthquakes.  The  absence  of  a  peaked 
spectrum  does  not  rule  out  dike  injection  as  the  source  of 
these  earthquakes,  but  it  does  rule  out  a  particular  model  of 
magmatic  intrusion,  namely,  that  of  Chouet  [1981]  and  Aki 
[1984],  for  a  low-viscosity  fluid.  As  Aki  [1984]  has  pointed 
out,  the  opening  of  a  channel  between  two  preexisting  fluid 
filled  cracks  would  largely  alleviate  this  problem  by  involving 
the  movement  of  only  minor  amounts  of  viscous  material, 
which  could  be  accomplished  in  an  adequately  small  time. 

Note  added  in  proof.  Barker  [1985]  has  shown,  using  theo¬ 
retical  seismograms,  that  low-speed,  near-surface  material  can 
cause  significant  amplification  of  near-source  ground  motion. 
If  this  has  occurred  in  the  data  discussed  here,  it  would  tend 
to  strengthen  the  conclusions  of  the  paper. 
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M,.  OBSERVATIONS  IN  THE  GREAT  BASIN  AND  A/„  VERSUS  ML 
RELATIONSHIPS  FOR  THE  1980  MAMMOTH  LAKES,  CALIFORNIA 

EARTHQUAKE  SEQUENCE 

By  David  E.  Chavez  and  Keith  F.  Priestley 


Abstract 

The  local  magnitude  distance  correction  term,  log  A0(A),  has  been  estimated 
for  paths  within  the  Great  Basin  by  converting  seismograms  of  earthquakes  in 
the  vicinity  of  Mammoth  Lakes,  California,  into  equivalent  Wood-Anderson  seis¬ 
mograms  and  measuring  the  decay  of  peak  amplitude  with  distance.  The  shape 
of  the  log  A0( A)  curve  which  results  suggests  more  attenuation  of  energy  in  the 
Wood-Anderson  passband  than  in  southern  California.  For  Great  Basin  paths  we 
find,  for  events  with  Mt  S  5.S,  that 
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where  A  is  epicentrat  distance,  R  =  V A1  +■  h*  is  hypocentral  distance,  and  h  is 
focal  depth.  Our  data  also  suggest  a  magnitude  bias  in  the  correction  term,  in 
that  for  events  above  Mt  «  5.5,  peak  amplitudes  recorded  at  epicentral  distances 
below  20  km  decay  less  rapidly  than  those  for  smaller  events  recorded  at  the 
same  sites.  This  observation  is  consistent  with  the  hypothesis  that  there  is  a 
near-source  saturation  of  ML  due  to  the  effects  of  a  finite  source  size  for  large 
earthquakes.  Using  the  revised  magnitude  scale  and  seismic  moments,  M0, 
estimated  from  spectra.  -  ialysis,  the  data  are  well  fit  by  the  straight  line 


log  Mo  =  (1.20  ±  0.05 )Mt  +  (17.49  ±  0.19) 


for  1  £  Ml  s  6.  M0  versus  ML  values  from  model  calculations  which  assume 
constant  stress  drop  predict  curved  moment-magnitude  plots  over  this  magnitude 
range.  The  fact  that  we  observe  a  linear  relationship  suggests  that  the  stress 
drop  of  these  events  is  not  constant  but  rather  increases  with  Mt,  in  particular  for 
events  above  Mt  =  5.5. 


Introduction 

Local  magnitude,  ML,  was  defined  by  Richter  (1935)  to  be 

Ml  =  log  A  -  log  A0(A) 

where  A  is  the  maximum  zero  to  peak  trace  amplitude  in  millimeters  recorded  on  a 
standard  Wood-Anderson  torsion  seismograph  (static  magnification  =  2800,  natural 
period  =  0.8  sec,  and  damping  factor  =  0.8)  at  an  epicentral  distance  A.  The  log 
A„(A)  function  removes  the  distance  dependence  from  the  amplitude,  incorporating 
geometrical  spreading,  change  in  wave  type,  scattering,  and  anelastic  attenuation. 
The  nature  of  these  last  three  will  depend  on  the  structure  between  the  source  and 
receiver  and  in  general  will  vary  between  regions.  Richter  empirically  determined 
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the  shape  of  the  log  .4,,(  A)  curve  for  southern  California  and  calibrated  it  so  that 
log  .1,(100)  =  -3,  i.e.,  a  A//,  zero  earthquake  would  have  a  peak  Wood-Anderson 
trace  amplitude  of  1  ^m  at  an  epicentral  distance  of  100  km.  Initially,  the  curve  was 
defined  for  the  range  25  <  A  <  000  km,  however  Gutenburg  and  Richter  (1942) 
later  extended  it  to  zero  distance  using  data  from  low-gain  torsion  seismographs. 
Values  for  log  A, , (A)  from  0  to  000  km  are  published  in  Richter  (1958,  p.  342). 

The  Richter  (1958)  log  A(.(A)  values  have  been  used  in  estimating  ML  outside  of 
the  southern  California  region  for  which  they  were  developed.  As  pointed  out  in  a 
number  of  studies,  this  invokes  the  tacit  assumption  that  the  attenuation  function 
is  the  same  as  in  southern  California.  This  will  not  always  be  the  case  however,  and 
consequently  equal  numerical  Ml  values  for  earthquakes  in  different  regions  will 
not  necessarily  imply  the  same  seismic  energy  release  at  the  source.  This  is  of 
interest  in  engineering  seismology  because  ML  [and  hence  log  A0(A)]  plays  an 
important  role  in  the  analysis  and  prediction  of  strong  ground  motion. 

In  this  paper,  we  use  events  in  the  vicinity  of  Mammoth  Lakes,  California, 
recorded  both  near  the  source  and  at  regional  distances  to  develop  a  log  A,.(  A)  curve 
for  the  Great  Basin.  We  are  fortunate  in  that  ray  paths  from  Mammoth  Lakes 
eastward  travel  exclusively  in  the  Great  Basin  while  paths  westward  travel  through 
or  near  areas  where  the  log  A0(A)  function  is  known.  This  allows  us  to  directly 
calibrate  the  attenuation  curve  to  that  for  which  Ml  is  defined. 

For  many  applications,  the  seismic  moment,  M0,  is  often  a  preferable  measure  of 
earthquake  size  since  moment  is  directly  related  to  specific  source  properties,  namely 
the  fault  area  and  amount  of  slip.  Moreover,  Ml  will  saturate  for  events  sufficiently 
large  that  the  corner  frequency  is  below  the  Wood-Anderson  passband  while  A/0 
will  not.  Theoretical  source  models  and  the  widespread  use  of  digital  recordings 
have  made  possible  the  routine  estimation  of  M0  for  a  large  range  of  sizes.  Here,  we 
develop  a  Mn  versus  ML  ...ationship  for  events  in  the  1980  Mammoth  Lakes 
sequence  and  compare  it  to  those  for  other  areas. 

Previous  Work 

Few  iod- Anderson  attenuation  functions  have  been  determined  for  areas 
outside  of  southern  California  primarily  because  of  the  lack  of  a  sufficient  number 
of  standard  seismographs  with  which  to  constrain  the  curve.  The  increasing  avail¬ 
ability  of  digital  seismograms  for  calibrated  stations  now  offers  the  possibility  of 
doing  so  by  creating  “synthetic”  Wood-Anderson  seismograms.  Bakun  and  Joyner 
(1984)  used  such  data  to  develop  a  log  A0{ A)  curve  for  central  California.  They 
found  that  the  original  southern  California  curve  was  appropriate  except  for  some 
modification  at  near  distances,  where  their  data  indicated  that  Richter’s  (1958) 
curve  did  not  completely  remove  the  distance  dependence. 

The  original  log  A0(A)  curve  has  also  been  reevaluated  for  southern  California. 
For  instance,  Luco  ( 1982),  Jennings  and  Kanamori  (1983),  and  Munguia  and  Brune 
(1984)  have  used  strong  motion  and/or  digital  recordings  to  determine  Ml  at 
distances  from  0  to  100  km  for  events  with  Ml  >  6.  Similar  to  the  results  of  Bakun 
and  Joyner  (1984),  all  found  that  Mt.  is  systematically  underestimated  at  distances 
less  than  50  km,  relative  to  that  determined  at  more  distant  stations.  Luco  (1982) 
and  Jennings  and  Kanamori  (1983)  interpreted  this  as  a  result  of  deficiencies  in 
the  original  log  A0(A)  function,  and  each  developed  revised  attenuation  curves  to 
remove  the  bias.  Munguia  and  Brune  (1984)  however,  hypothesize  that  the  under¬ 
estimation  of  Ml  for  large  events  is  a  consequence  of  finite  source  effects.  They 
describe  a  numerical  experiment  which  suggests  that  at  large  epicentral  distances, 
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arrivals  from  various  portions  of  the  causative  fault  constructively  interfere  to 
produce  larger  amplitude  arrivals  and  hence  yield  a  greater  magnitude.  At  near 
distances,  the  arrivals  interfere  destructively,  leading  to  smaller  amplitudes  and  a 
lower  estimate  of  A//.. 

Munguia  and  Brune  noted  the  opposite  trend  when  considering  smaller  earth¬ 
quakes  (3.0  ^  Mi.  =  5.5).  Recordings  within  an  epicentral  distance  of  10  km  or  less 
yielded  magnitudes  which  were  up  to  one  unit  greater  than  those  at  larger  distances. 
They  suggested  that  this  was  due  to  site  amplifications  as  well  as  high-frequency 
source  spectra  not  adequately  accounted  for  by  Richter’s  (1958)  log  A„(A)  curve. 

Boore  and  Hutton  ( 1984)  used  regression  analysis  on  several  thousand  actual  and 
simulated  Wood- Anderson  data  to  determine  a  revised  southern  California  log  A„(A) 
curve  for  the  distance  range  10  ^  A  =  560  km.  They,  as  well  as  Bakun  and  Joyner 
(1984),  modeled  the  decay  in  peak  Wood- Anderson  amplitude  as  a  combination  of 
material  attenuation  and  geometrical  spreading  using 

e->* 

A.-C—  (I) 


where  y  and  t\  are  the  material  attenuation  and  geometrical  spreading  coefficients, 
respectively,  R  =  Va2  +  h2  is  hypocentral  distance,  h  is  focal  depth,  and  C  is  a 
constant.  [Although  Richter  (1935)  defined  IogA0(A)  in  terms  of  epicentral  distance, 
he  did  so  assuming  that  the  events  he  studied  had  similar  focal  depths.  It  becomes 
necessary  to  use  hypocentrJ  distance  in  order  to  eliminate  the  distance  bias  from 
near  source  ML  estimates  if  one  considers  earthquakes  at  a  variety  of  depths 
(Jennings  and  Kanamori,  1983)].  Boore  and  Hutton  (1984)  found  that  y  =  0.0046 
km-1  for  7?  =  1.  A  more  recent  solution  (Boore,  personal  communication,  1985)  has 
y  =  0.0035  km-1  and  7j  =  1.09.  In  this  paper,  all  reference  to  Boore  and  Hutton 
(1984)  should  be  taken  to  mean  their  more  recent  results,  given  above.  Their  results 
suggest  less  attenuation  for  distances  beyond  200  km  than  what  is  indicated  by 
Richter’s  (1958)  curve.  At  distances  closer  than  50  km,  they  obtain  a  curve  similar 
to  those  proposed  by  Luco  (1982)  and  Jennings  and  Kanamori  (1983).  Thus,  Boore 
and  Hutton  (1984)  and  Bakun  and  Joyner  (1984)  find  a  similar  near-source  bias  for 
both  large  and  small  earthquakes,  a  result  which  differs  from  that  of  Munguia  and 
Brune  (1984). 


Data  and  Analysis 

The  data  analyzed  in  this  study  consist  of  seismograms  of  earthquakes  in  the 
vicinity  of  Mammoth  Lakes,  California,  primarily  events  in  the  1980  sequence.  The 
data  at  distances  less  than  30  km  include  velocity  and  acceleration  records  obtained 
by  the  U.S.  Geological  Survey  (Archuleta  et  al„  1982),  and  strong  motion  accelero¬ 
grams  from  instruments  installed  by  the  California  Divisions  of  Mines  and  Geology 
(Turpen,  1980)  and  by  the  University  of  Southern  California  and  the  California 
Institute  of  Technology  (Moslem  et  al.,  1983).  For  distances  from  90  to  600  km,  we 
have  data  from  broadband  stations  operated  by  the  Lawrence  Livermore  National 
Laboratory.  Similar  broadband  instruments  installed  by  the  University  of  Nevada 
recorded  aftershocks  of  the  23  November  1984  ML  =  5.7  earthquake  near  Bishop, 
California.  We  have  also  used  seismograms  from  the  Wood-Anderson  seismograph 
at  Dugway,  Utah.  In  all,  479  seismograms  for  88  earthquakes  in  the  range  2.5  £  ML 
=  6  were  studied.  Figure  1  is  a  map  showing  the  source  regions  and  station  locations 
for  the  data  analyzed. 
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Following  Bakun  et  al.  (1978),  we  synthesized  Wood-Anderson  seismograms  by 
transforming  the  original,  horizontal  component  time  series  into  the  frequency 
domain  and  div  iding  by  the  instrument  response,  thus  converting  the  data  to  ground 
displacement.  This  was  then  multiplied  by  a  Wood-Anderson  response,  and  the 
result  was  transformed  back  into  the  time  domain  to  give  a  “synthetic”  Wood- 
Anderson  seismogram,  i.e.,  that  which  would  have  been  recorded  at  the  site  had  the 
original  instrument  been  replaced  with  a  Wood-Anderson  torsion  seismograph  with 
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°f  'h.e  we«ern  United  States  showing  the  Mammoth  Lakes,  California,  source  region 
'■  eP‘ce'Jters i  plotted  as  dots),  and  seismic  stations  used  in  this  study.  The  triangles  indicate 
of  I hlrJir™  Uivd  <Jet5rrmne  ,he  Great  Basm  attenuaticn  curve,  and  the  squares  give  the  locations 
?k.  nn»h  ,rfrlf  Wood-Anderson  seismographs  used  in  determining  the  reference  magnitudes  (UCB  to 
the  north,  CIT  to  the  south).  The  heavy  line  delineates  the  border  of  the  Basin  and  Range 

sufficient  dynamic  range.  We  determined  the  average  of  the  maximum  zero  to  peak 
amp  itudes  of  the  two  horizontal  components  at  each  station,  and  these  average 
amp  itudes  and  their  associated  epicentral  distances  constitute  our  data  set.  These 
a  a,  after  normalization  described  below,  are  plotted  in  Figure  2.  Figure  2a  displays 

Th  etlt,re  data  Set  While’  f0r  clarity’  Fisure  2b  shows  the  data  at  close  distan  es 
ine  shape  of  the  curve  implied  by  the  data  points  describes  the  Great  Basin  log 
/toiAJ  Junction. 
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Station  corrections.  The  amplitudes  at  close  distances  (A  <  30  km)  are  from  a 
variety  of  stations  with  differing  site  geology  and  presumably  different  site  responses 
(enlargement,  Figure  1).  We  attempted  to  remove  any  apparent  trends  from  the 
amplitude  versus  distance  plot  which  may  be  a  consequence  of  the  combination  of 
station  distribution  and  site  amplification  by  applying  corrections  determined  in 
the  following  manner.  Since  the  raw  amplitude  data  include  both  site  and  distance 
effects,  we  first  reduced  the  distance  dependence  by  correcting  the  data  using  Boore 
and  Hutton’s  (1984)  log  A0(±)  curve  with  the  assumption  that  the  shape  of  the 
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Epicentral  Distance  (km) 

Fic.  2.  Normalized  peak  Wood-Anderson  amplitudes  versus  distance  recorded  at  Great  Basin  sites. 
The  squares  indicate  amplitudes  for  events  with  reference  magnitudes  below  5.5,  the  stars  are  for 
magnitude  5.5  and  above  events.  Lines  give  the  indicated  southern  California  log  Ao(A)  functions. 

actual  log  /40(A)  curve  resembles  theirs  at  that  distance  range.  Then,  assuming  that 
a  station  on  granite  required  no  correction,  the  stations  were  grouped  by  site  geology 
(granite,  rhyolite,  rhyodacite,  glacial  till,  lake  deposits,  or  alluvium),  and  the  average 
magnitude  bias  was  determined  relative  to  the  granite  site.  These  corrections  were 
applied  to  the  raw  data,  and  a  log  /f0(^)  function  was  determined  for  the  corrected 
data  using  the  method  described  in  the  next  section.  The  above  procedure  was 
repeated,  this  time  using  the  new  log  A0(A)  curve  and  removing  the  constraint  that 
the  granite  correction  was  zero.  The  resulting  geology-based  corrections  ranged 
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from  -0.16  for  the  granite  site  to  approximately  -0.50  for  alluvium  and  other 
unconsolidated  sediments.  These  corrections  are  included  in  all  data  presented  here. 

As  it  turns  out,  the  shape  of  the  corrected  amplitude  versus  distance  curve  is  not 
significantly  different  from  the  uncorrected,  implying  that  distance  dependence  is 
responsible  for  the  major  trend  of  the  curve.  At  distances  beyond  90  km,  the  data 
from  each  station  plot  as  a  group,  isolated  from  the  others,  causing  individual 
station  biases  to  be  easily  recognized.  For  this  reason,  we  decided  it  was  not  necessary 
to  estimate  station  corrections  for  the  more  distant  sites  at  this  point  in  the  analysis. 

Normalization.  The  data  in  Figure  2  have  been  normalized  relative  to  the  average 
peak  Wood-Anderson  amplitude  at  California  stations  for  the  same  events.  This 
was  done  by  subtracting  the  average  ML  as  reported  by  the  University  of  California 
at  Berkeley  (UCB)  and  the  California  Institute  of  Technology  (CIT)  earthquake 
catalogs,  thus  at  the  same  time  calibrating  the  amplitude  versus  distance  data  to 
Richter’s  (1958)  log  A0(A)  curve  which  both  institutions  used  in  computing  Ml- 

Although  there  has  not  been  a  quantitative  evaluation  of  the  nature  of  the  log 
Ao(S)  curve  for  paths  between  Mammoth  Lakes  and  the  UCB  Wood-Anderson 
seismographs  (Figure  1),  our  data  suggest  that  the  apparent  -ttenuation  is  very 
similar  tc  that  along  paths  from  Mammoth  Lakes  to  the  CIT  Wood-Andersons.  If 
this  were  not  the  case,  then  one  would  expect  to  see  a  systematic  difference  between 
the  magnitudes  as  reported  by  the  two  institutions.  Figure  3  is  a  plot  of  the  reported 
UCB  versus  CIT  magnitudes  for  the  55  events  common  to  both  catalogs.  The  data 
scatter  about  a  line  of  unit  slope  and  zero  intercept,  indicating  that  no  path  bias  is 
present. 

The  scatter  about  the  line  in  Figure  3  shows  that  the  reference  magnitudes  (the 
average  of  the  two  catalog  values)  could  be  in  error  by  approximately  ±0.3  unit  [not 
taking  into  account  the  systematic  error  introduced  by  inaccuracies  in  Richter’s 
(1958)  curvej.  i  uis  uncertainty  in  the  reference  magnitude  is  transferred  to  the 
amplitude  data  when  they  are  normalized,  howevei  it  does  not  significantly  affect 
any  of  our  conclusions.  The  difference  between  Boore  and  Hutton’s  (1984)  revised 
log  Aa(A)  curve  and  Richter’s  version  is  about  0.3  at  350  km,  the  average  distance 
from  Mammoth  Lakes  to  the  UCB  and  CIT  Wood-Anderson  seismographs.  This 
difference  was  removed  from  the  reference  magnitudes  prior  to  normalizing  the 
data.  For  the  remainder  of  this  paper,  any  mention  of  “reference  magnitude”  will 
mean  the  UCB  and  CIT  average  minus  0.3. 

Estimation  of  Log  A0(A)  for  the  Great  Basin 

Figure  2  shows  the  normalized  data  including  the  final  station  corrections  de¬ 
scribed  in  the  previous  section.  For  comparison,  three  versions  of  the  southern 
California  log  A,(A)  curve  are  included  in  the  plot.  These  are  Richter’s  (1958) 
original  curve,  Jennings  and  Kanamori’s  (1983)  modifications  to  it,  and  Boore  and 
Hutton’s  (1984)  revised  curve.  In  Figure  2  (as  well  as  in  Figures  4  and  5),  stars  are 
used  to  indicate  data  from  events  with  a  reference  magnitude  above  5.5  while 
squares  denote  those  data  from  smaller  earthquakes. 

Figure  2  reveals  that  there  exists  a  difference  beween  the  southern  California  and 
Great  Basin  log  .4o(A)  functions.  At  distances  beyond  200  km  (Figure  2a),  peak 
Wood-Anderson  amplitudes  decay  more  rapidly  than  Boore  and  Hutton’s  (1984) 
curve  predicts.  The  data  at  that  distance  range  are  better  fit  by  Richter’s  (1958) 
curve.  At  near  distances  (Figure  2b),  the  data  points  at  distances  below  10  km  fall 
above  the  various  southern  California  curves,  but  merge  with  them  for  10  §  A  =  30 
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km.  The  trend  of  the  data  is  oblique  to  the  curves,  indicating  more  rapid  attenuation 
of  energy  than  predicted  by  the  southern  California  log  A0(A)  functions. 

The  data  also  suggest  that  the  amplitude  decay  at  close  distances  is  a  function  of 
magnitude.  The  points  in  Figure  2b  for  Mt.  =  5.5  events  (stars)  plot  below  the  main 
population  (squares)  and  show  little  variation  in  amplitude  with  distance  over  the 
interval  0  ^  A  Si  30  km.  These  points  tend  to  scatter  about  the  attenuation  curve 
proposed  by  Jennings  and  Kanamori  (1983).  This  observation  cannot  be  a  conse¬ 
quence  of  site  effects,  since  the  same  stations  recorded  both  large  and  small 
earthquakes.  This  is  consistent  with  the  hypothesis  of  Munguia  and  Brune  (1984) 
that  there  is  a  near-source  Ml  saturation  due  to  Finite  source  effects  for  large  events. 
In  view  of  this  apparent  magnitude  bias  in  the  log  40(  A)  curve,  only  those  data  from 
events  with  Ml  <  5.5  were  considered  when  developing  the  Great  Basin  Wood- 
Anderson  attenuation  curve. 

Figure  2a  shows  that  the  data  at  100  km  scatter  about  -3,  the  point  through 


L'CB  Magnitude 

Fig.  3.  Plot  of  Ml  as  reported  by  VCB  versus  that  reported  by  CIT.  The  dashed  line  has  unit  slope 
and  zero  intercept. 


which  the  southern  California  curves  were  constrained  to  pass.  Because  of  this,  we 
decided  to  also  constrain  the  Great  Basin  log  A0(A)  curve  in  a  like  manner.  (There 
is  no  a  priori  reason  to  assume  that  this  would  be  the  case.  Energy  from  a  magnitude 
zero  earthquake  propagating  through  a  medium  whose  attenuation  properties  differ 
substantially  from  those  of  southern  California  will  require  log  A0(100)  to  have  a 
value  other  than  -3.) 

We  used  the  model  given  by  equation  (1)  to  describe  the  observed  spatial 
attenuation.  There  is  a  trade-off  between  7  and  7,  so  we  assumed  that  7  =  1  for 
0  g  A  S  90  km  and  7  =  0.83  for  90  g  A  s  600  km.  The  value  of  7  =  1  represents 
body  wave  geometrical  spreading  for  sources  buried  in  a  homogeneous  medium, 
while  7  =  0.83  is  the  theoretical  time  domain  geometrical  spreading  coefficient  of  a 
dispersed  Airy  phase  which  Campillo  et  at.  (1984)  have  numerically  verified  to  be 
applicable  for  the  Lg  phase.  Lg  is  commonly  the  largest  amplitude  phase  on  our 
simulated  Wood-Anderson  seismograms  at  regional  distances. 
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In  Figure  4  are  the  same  data  as  in  Figure  2,  superimposed  on  curves  given  by  (1) 
for  a  variety  of  y,  and  constrained  to  equal  -3  at  100  km.  These  curves  indicate  the 
range  in  y  which  is  possible,  given  the  scatter  in  the  data.  At  distances  beyond  90 
km  (Figure  4a),  y  can  vary  between  0.004  and  0.009  km-1  with  a  (visually)  preferred 
value  of  0.006  km"1. 

At  the  closer  distances  (Figure  4b),  the  data  lie  within  lines  corresponding  to 
0.005  ^  y  ^  0.025  km"1,  and  are  best  fit  by  a  y  of  0.016  km'1.  These  curves  assume 
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Fir.,  4.  The  same  data  as  in  Figure  2,  superimposed  on  curves  given  by  equation  (1).  (a)  Assuming  a 
geometrical  spreading  coefficient  of  ij  =  0.83.  (b)  Assuming  rj  =  1  and  a  2-km  focal  depth.  The  material 
attenuation  coefficients,  are  shown  for  each  curve;  the  preferred  fits  to  the  data  ore  indicated  by  a 
heavy  line.  The  dotted  line  in  (b)  is  Boore  and  Hutton's  (1984)  curve  assuming  a  5-km  focal  depth. 


p  focal  depth  of  2  km,  chosen  since  the  trends  which  result  more  closely  match  that 
of  the  data  than  what  was  obtained  using  other  reasonable  focal  depths.  The  median 
depth  which  was  actually  computed  in  the  earthquake  locations  was  5  km,  however 
attenuation  curves  computed  using  that  depth  are  flatter  than  required  by  the  data 
for  epicentral  distances  below  10  km.  This  is  illustrated  by  the  dotted  line  in  Figure 
4b  which  assumes  a  5  km  focal  depth  (Boore  and  Hutton’s,  1984,  curve).  The 
implications  of  this  is  that  either  our  model  for  describing  the  attenuation  [equation 
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( 1 )]  is  :,i  adequate  at  very  near  distances  or  that  the  actual  focal  depths  are  shallower 
than  obtained.  The  former  being  the  case  would  preclude  our  ascribing  any  physical 
interpretation  to  the  constants  in  the  attenuation  function.  The  curve  empirically 
describes  the  spatial  attenuation  however,  and  can  be  used  in  computing  A//.  without 
regard  to  the  validity  of  equation  (1 ).  Based  on  Figure  4,  our  estimate  of  the  Great 
Basin  log  ,4„(A)  function  is 

!■  -O.UIfi/f" 

login  — - —  -  0.31  0  s  A  5  90  km 

A 

(2) 
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log,,,  -jpnu~  -  L08  90  g  A  g  600  km. 
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We  computed  a  ML(syn)  by  applying  equation  (2)  to  the  data,  and  Figure  5  shows 
the  differences  between  Mt.(syn )  and  the  reference  magnitudes,  called  M/.(wa),  as  a 
function  of  distance.  Jennings  and  Kanamori’s  (1983)  curve  was  used  for  the  data 
from  events  with  reference  magnitude  above  5.5  and  A  <  30  km.  The  points  scatter 
about  the  zero  reference  line  showing  that,  on  the  average,  the  distance  dependence 
has  been  removed  from  the  data.  The  Figure  has  a  suggestion  of  a  systematic 
overestimation  of  ML  at  distances  beyond  500  km,  however,  we  feel  this  is  an 
illus  >ry  consequence  of  the  need  for  a  station  correction  at  600  km  (Dugway). 
Dugway  is  the  only  station  beyond  90  km  which  is  situated  on  sediments;  all  the 
others  are  installed  on  competent  bedrock.  Although  we  cannot  exclude  the  possi¬ 
bility  that  these  other  stations  require  some  correction  (Tucker  et  ai,  1984),  they 
should  be  small  relative  to  that  at  Dugway.  Figure  5  indicates  a  correction  for 
Dugway  of  —0.41. 

Applying  Jennings  and  Kanamori's  (1983)  attenuation  curve  to  the  data  for 
events  above  reference  magnitude  5.5  causes  them  to  fall  within  the  main  population 
of  data  points.  Use  of  equation  (2)  instead  of  Jennings  and  Kanamori’s  log  A0(  A) 
values  would  ha’  e  preserved  the  difference  between  the  large  and  small  events 
evident  in  Figures  2b  and  4b.  We  have  considerably  less  data  for  the  larger  events 
and  so  are  unable  to  unequivocally  state  the  need  for  a  magnitude  dependent 
attenuation  curve.  Our  data  do,  however,  indicate  that  extra  care  must  be  taken 
when  determining  Mi  at  close  distances.  Given  the  data  currently  available,  we 
recommend  that  for  Great  Basin  paths  and  epicentral  distances  under  30  km,  the 
log  A,, (A)  function  given  by  equation  (2)  be  used  only  for  events  below  Ml  ~  5.5. 
Jennings  and  Kanamori’s  (1983)  curve  appears  more  appropriate  for  larger  events 
at  that  distance  range. 

Figure  6  compares  the  attenuation  curve  found  for  Great  Basin  paths  to  Richter  s 
(1958)  southern  California  curve,  Jennings  and  Kanamori’s  (1983)  modifications  to 
it,  and  Boore  and  Hutton’s  (1984)  revised  curve.  Based  on  the  results  of  Bakun  and 
Joyner  (1984),  Richter’s  (1958)  curve  at  distances  beyond  100  km  can  also  be  viewed 
as  representing  attenuation  in  central  California.  We  Find  that  in  the  Great  Basin, 
the  attenuation  of  energy  with  frequencies  in  the  Wood-Anderson  passband  is 
greater  than  in  southern  California  but  slightly  less  than  that  in  central  California. 
The  scatter  in  the  data  allow  for  equal  attenuation  in  the  Great  Basin  and  central 
California,  however.  The  trend  of  the  data  at  close  distances  also  suggests  stronger 
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Fig.  5.  The  difference  between  estimates  of  Ml  made  from  Great  Basin  equivalent  Wood-Ander*ons, 
A/tisjn),  and  the  California  reference  magnitudes,  ML(wa).  The  amplitudes  from  events  below  5.5 
(squares)  were  corrected  using  equation  (2);  those  from  larger  events  (stars)  were  corrected  using 
Jennings  and  Kanamori's  curve  for  epicentral  distances  below  30  km,  and  equation  (2)  for  greater 
distances 
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Fig  8  Comparison  of  the  log  A0(A)  curve  obtained  for  the  Great  Basin  with  Kichler’s  ( 1958)  southern 
California  curve.  Jennings  and  Kanamori's  (19831  modification  to  it.  and  Boore  and  Hutton's  (1984) 
revised  curve  I'he  question  marks  indicate  the  distance  range  for  which  no  Great  Basin  data  were 
available 
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attenuation  than  given  by  any  of'  the  southern  California  curves.  We  have  no  data 
for  the  disiar.ee  range  30  A  =  90  km;  t he  curve  shown  is  a  continuation  of  that 
found  for  A  =  30  km. 

'['he  material  attenuation  coefficients  are  related  to  the  seismic  quality  factor,  Q, 
through 


where  f  is  the  frequency  of  the  wave,  and  U  is  its  group  velocity.  Thus,  we  can 
estimate  the  Q  for  Lg  in  the  Great  Basin  by  assuming  U  =  3.5  km/sec  (the  typical 
Lg  group  velocity)  and  /=  1.25  Hz  (the  Wood-Anderson  natural  frequency).  Given 
the  range  in  y  indicated  in  Figure  4a,  we  find  that  the  Great  Basin  LgQ at  frequencies 
near  1  Hz  is  between  125  and  280,  with  a  preferred  value  of  190.  Boore  and  Hutton’s 
(1984)  curve  suggests  that  LgQ  in  southern  California  is  closer  to  our  upper  limit. 
For  comparison,  Patton  and  Taylor  (1984)  examined  surface  waves  to  conclude  that 
the  shear  wave  Q  in  the  Great  Basin  was  about  100  in  the  lower  crust  and  30  in  the 
upper  mantle,  while  Singh  and  Herrmann  (1983)  found  a  coda  Q  of  200  for 
California.  The  value  of  y  -  0.016  km-1  obtained  at  close  distances  implies  a  1  Hz 
shear  wave  Q  of  75,  assuming  U  —  3.3  km/sec. 

Our  conclusion  that  attenuation  in  the  Great  Basin  is  less  than  in  southern 
California  is  in  conflict  with  Evernden’s  (1975)  study  of  intensities  in  the  Great 
Basin  which  concluded  that  it  has  slightly  less  attenuation  than  California.  Part  of 
this  discrepancy  could  be  explained  by  the  fact  that  central  California  attenuation 
does  appear  to  be  greater  than  that  in  the  Great  Basin.  Whatever  the  case,  the  data 
require  a  systematic  increase  in  attenuation  relative  to  southern  California,  in  order 
to  bring  our  ML  estimates  into  agreement  with  those  made  using  the  UCB  and  CIT 
Wood-Anderson  recordings  of  the  same  earthquakes  and  Boore  and  Hutton’s  (1984) 
log  An( A)  function. 

A  recent  development  which  affects  studies  of  this  type  is  the  suggestion  that 
there  has  been  a  systematic  miscalibration  of  California  Wood-Anderson  seismo¬ 
graphs  (e  g.,  Bakun  and  Joyner,  1984).  It  appears  possible  that  these  instruments 
operate  at  a  gain  of  2100  rather  than  2800,  as  is  required  by  the  definition  of  Ml- 
This  presumably  was  also  the  case  for  the  instruments  used  by  Richter  (1935)  which 
would  indicate  that  the  “standard”  Wood-Anderson  gain  is  actually  2100.  Since  it 
has  not  yet  been  firmly  established  in  the  literature  that  this  problem  exits,  we  have 
used  a  gain  of  2800  when  converting  the  data.  Thus,  there  is  the  possibility  of  a 
systematic  error  of  approximately  0.12  in  the  constants  in  equation  (2).  We  also 
point  out  that  the  attenuation  curve  developed  for  the  Great  Basin  is  tied  to  the 
southern  California  curve  through  our  normalization  procedure.  Future  modifica¬ 
tions  to  the  southern  California  log  A0(A)  function  will  necessitate  similar  modifi¬ 
cations  to  the  Great  Basin  curve. 

Versus  Ml  Relationships  for  the  1980  Mammoth  Lakes  Sequence 

Using  the  log. AJ A)  curve  developed  here  for  the  smaller  events  (ML  «  5.5)  and 
Jennings  and  Kanamori’s  (1983)  curve  (at  close  distances)  for  the  larger  events,  we 
recomputed  A/;,  for  65  earthquakes  of  the  1980  Mammoth  Lakes,  California, 
sequence  for  which  estimates  of  the  seismic  moment  were  available.  The  data 
include  events  used  to  determine  the  revised  attenuation  curve  plus  17  others  with 
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0.5  £  A//.  £  2.8  which  were  not  used  above  since  no  reference  magnitudes  were 
available.  Figure  7  is  a  plot  of  these  revised  magnitudes  versus  the  seismic  moments 
taken  from  Archuleta  et  al.  (1982)  and  Priestley  et  al.  (1985).  These  moments  are 
based  on  spectral  analysis  of  the  same  local  digital  and  strong-motion  recordings 
from  which  the  synthetic  Wood-Anderson  seismograms  were  derived.  Priestley  et 
al.  ( 1985)  also  determined  the  moments  of  the  larger  {ML  S  4.0)  events  from  regional 
surface  wave  recordings  and  found  them  to  be  in  agreement  with  those  determined 
from  spectral  analysis  of  the  local  recordings. 

Two  straight  lines  were  fit  to  the  data.  For  events  in  the  magnitude  range  1  £ 
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Fig  7.  Log  \f„  versus  Mi.  for  earthquakes  in  the  1980  Mammoth  Lakes,  California,  sequence.  One 
standard  deviation  ol  the  mean  error  bars  are  given  for  events  with  multiple  recordings.  The  lines  are 
least  squares  fits  to  the  data,  and  the  stars  are  theoretical  values  taken  from  Hankb  and  Boore  (1984). 


Mi  S  6,  the  data  are  fit  by  the  line 

log  Af„  =  (1.20  ±  0.05) Aft  +  (17.49  ±  0.19).  (3) 

A  lew  of  the  smaller  events  (Af;  S  1.5)  deviate  from  this  line.  For  events  in  the 
magnitude  range  0.5  =  M,  =  3,  the  moments  scatter  about  the  line 


log  M0  =  (0.81  ±  0.08)Aft  +  (18.37  ±  0.18), 

however,  it  is  strongly  affected  by  the  three  smallest  events.  There  are  very  few 
recordings  for  these  three  events,  and  it  could  he  easily  argued  that  neither  their 
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nvmu'nts  n>ir  their  magnitudes  are  well  determined.  Furthermore,  they  were  re¬ 
corded  .it  vers  close  epicentral  distances,  and  it  is  possible  that  variations  m  local 
depth  are  affecting  the  results.  Hanks  (personal  communication,  1985)  has  also 

pointed  out  that  a  / . .  (Hanks,  1982)  of  10  Hz  would  cause  the  magnitude  of  those 

events  to  he  underestimated  by  about  0.8  magnitude  unit,  which  is  enough  to  bring 
the  points  i  lose  the  line  given  by  equation  (8).  For  these  reasons,  we  feel  that  the 
deviation  of  those  data  from  the  rest  is  not  particularly  significant.  If  we  neglect 
those  points,  then  the  moment-magnitude  data  are  linearly  related  over  the  entire 
magnitude  range  1  S  Af,  g  6,  and  do  not  exhibit  any  curvature  of  the  type  noted 
by  Bakun  (1984)  or  Hanks  and  Boore  (1984)  for  other  California  earthquakes. 

fhe  A/,,  versus  A/,,  relationship  given  in  equation  (3)  can  be  compared  to  that 
given  by  Archuleta  ef  al.  (1982)  who  used  data  primarily  in  the  range  3  =  A/;.  =  5 
with  only  one  event  above  ML  =  5.  Our  data,  on  the  other  hand,  span  the  range  0.5 
=  Ml  =  0  with  seven  events  above  Ml  —  5.  Archuleta  et  al.  (1982)  found  that  the 
M„  versus  M,  relationship  for  the  40  events  for  which  they  had  UCB  Mi.  determi¬ 
nations  was,  for  the  range  2.9  Si  Ml  =i  6.2, 

log  Mu  =  (0.96  ±  0.06)A/t  +  ( 18.14  ±  0.28) 

but  for  the  larger  events  with  more  reliable  ML  estimates  (3.5  £  ML  §  6.2),  it  was 

log  Mo  =  (1.05  ±  0.08) Aft  +  (17.76  ±  0.33). 

The  greater  slope  in  equation  (3)  is  due  to  the  fact  that  we  used  larger  moments  for 
events  with  Mi.  >  5.0  than  did  Archuleta  et  al.  (1982).  For  example,  for  the  27  May 
1980  event  at  14:50,  Archuleta  et  al.  (1982)  found  a  seismic  moment  of  2.33  x  1024 
dyne-cm  based  on  two  recordings  whereas  Priestley  et  al.  (1985)  obtained  a  value  of 
7.24  x  1024  dyne-cm  when  using  nine  recordings. 

In  comparing  this  result  with  other  regions  of  California,  the  slope  in  equation 
(3)  is  less  than  the  least-squares  fit  Archuleta  et  al.  (1982)  made  to  the  data  in 
Thatcher  and  Hanks  (1973),  but  is  very  similar  to  the  relation  found  by  Bakun  and 
Lindh  (1977) 


log  Mo  =  (1.21  ±  0.03 )Ml  +  (17.02  ±  0.07) 

for  the  Oroville,  California,  earthquake  sequence.  Oroville,  in  the  western  Sierra 
Nevada,  is  located  in  a  structural  environment  which  is  more  similar  to  that  of 
Mammoth  Lakes  than  the  southern  California  region  studied  by  Thatcher  and 
Hanks  ( 1973). 

In  a  recent  study  of  five  different  source  regions  (Parkfield,  San  Juan  Bautista, 
the  Sargent  fault,  Coyote  Lake,  and  the  Livermore  Valley),  Bakun  (1984)  detected 
an  infection  in  the  log  Mn  versus  ML  plot  at  about  Ml  =  3.  For  the  events  in  the 
range  1.5  s  Af/.  ^  3.5,  he  found  the  relationship  to  be 

log  Mn  =  l-2ML  +  17 


which  is  essentially  the  same  as  the  relationship  given  in  equation  (3)  based  on  data 
in  the  range  1  =  Ml  =  6.  However,  for  events  in  the  similar  range  3.5  ^  Mi.  ^  6.25, 
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Bakun  (1984)  found  the  relationship  to  be 

log  Mo  =  1.5  Ml  +  16. 

Bakun  (1984)  and  Hanks  and  Boore  (1984)  have  noted  that  moment  versus 
magnitude  plots  for  California  earthquakes  exhibit  positive  curvature  when  the 
range  of  magnitudes  considered  is  sufficiently  large.  Hanks  and  Boore  were  able  to 
model  this  curvature  using  the  method  described  by  Boore  (1983).  In  Boore  (1983), 
seismograms  of  hypothetical  earthquakes  are  generated  which  have  an  amplitude 
spectrum  equal  to  that  predicted  by  the  Brune  (1970,  1971)  source  model,  given 
specified  source  parameters.  Boore  assumed  a  constant  “stress  parameter”  of  100 
bars,  thus  restricting  the  spectrum  scaling  law  to  be  a  function  of  seismic  moment 
alone.  The  stress  parameter  controls  the  strength  of  the  high-frequency  radiation, 
and  Hanks  and  Boore  (1984)  have  referred  to  it  as  the  am,  stress  drop. 

Hanks  and  Boore  (1984)  computed  M0  and  ML  from  these  “seismograms”  and 
found  that  their  assumption  of  constant  stress  drop  allowed  them  to  reproduce  the 
observed  curvature  in  the  moment  versus  magnitude  plot  for  0  =  ML  ^  6.  Their 
theoretical  values  are  included  in  Figure  7.  They  show  only  slight  curvature  for 
magnitudes  below  ML  *  5,  but  beyond  that  the  seismic  moment  increases  ever  more 
rapidly  with  increasing  Ml-  This  phenomenon  is  not  present  in  the  observed  data, 
indicating  that  the  assumption  of  constant  stress  drop  independent  of  magnitude  is 
not  valid  for  the  1980  Mammoth  Lakes  earthquakes,  over  the  range  1  £  ML  5  6. 
However,  given  the  small  curvature  in  the  theoretical  values  below  ML  *  5,  together 
with  the  scatter  in  the  observations,  the  data  are  consistent  with  constant  stress 
drop  for  the  smaller  events.  The  data  fall  to  the  left  of  the  theoretical  values  for  ML 
<  5.  This  fact  suggests  that  the  stress  drops  for  these  events  are  ler~  than  the  100 
h  ’rs  assumed  by  Hanks  and  Boore  (1984)  in  their  computations.  Likewise,  Figure 
7  suggests  that  events  larger  than  Ml  ~  5.5  had  stress  drops  above  100  bars.  This 
is  in  agreement  with  results  from  spectral  analysis  of  Archuleta  et  al.  (1982)  and 
Priestley  et  al.  (1985)  who  found  that  the  largest  events  had  stress  drops  on  the 
order  of  1000  bars. 

Teleseismically  determined  estimates  of  seismic  moment  for  the  three  largest 
events  are  0.5  to  one  order  of  magnitude  greater  than  the  locally  determined 
estimates  shown  in  Figure  7  (Given  et  al.,  1982;  Barker  and  Langston,  1983).  While 
this  discrepancy  has  yet  to  be  explained,  it  does  indicate  the  possibility  that  the 
actual  moments  of  the  larger  events  are  greater  than  what  we  used.  If  this  is  the 
case,  then  the  moment-magnitude  plot  would  exhibit  curvature  at  larger  magnitudes. 

Summary 

We  have  developed  a  log  A„( A)  function  which  removes  the  distance  dependence 
from  local  magnitude  estimates  using  Great  Basin  recordings  of  Ml  =  5.5  earth¬ 
quakes.  The  data  indicate  greater  attenuation  of  energy  in  the  Wood-Anderson 
passband  than  in  southern  California,  but  equal  or  slightly  less  attenuation  than  in 
central  C alifornia.  Peak  Wood-Anderson  amplitudes  of  events  larger  than  ML  ~  5.5 
decay  less  rapidly  than  those  for  smaller  events  at  close  distances,  consistent  with 
near-source  saturation  of  Ml  due  to  the  finite  source  size  of  the  larger  earthquakes. 
A  revised  moment-magnitude  relationship  for  the  1980  Mammoth  Lakes  earth¬ 
quakes  is  similar  to  other  Sierra  Nevada  (Oroville)  earthquakes.  Comparison  of  the 
moment -magnitude  data  to  theoretical  values  suggests  that  events  below  Ml  ~  5 
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h;ul  constant  stress  drop  while  larger  events  had  increasingly  greater  stress  drops, 
if  we  use  locally  determined  estimates  of  seismic  moment. 
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SCALING  OFMi.mj  AND  Ms  TO  THE  SEISMIC  SPECTRUM 
FOR  THE  MAMMOTH  LAKES,  CALIFORNIA  EARTHQUAKES 
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SUMMARY  :  In  this  paper  we  compare  the  relative  excitation  of  low  frequency 
(0.05  to  0.125  Hz)  and  high  frequency  (0.5  to  10.0  Hz)  for  several  moderate 
earthquakes  {Mi  4  -  5.5)  in  the  1980  Mammoth  Lakes,  California  sequence.  Such 
studies  have  importance  for  both  source  mechanism  and  seismic  decrimination 
studies.  We  have  determined  local  magnitude  {Mi),  teleseismic  body  wave  mag¬ 
nitude  (m6).  and  surface  wave  magnitude  ( Ms )  for  a  large  number  of  events  in 
the  magnitude  range  1  to  6.2.  Comparison  of  spectral  amplitude  U0  derived 
from  body  waves  recorded  at  near  epicentral  distances  and  derived  from  sur¬ 
faces  waves  recorded  at  regional  distances,  are,  within  the  errors, the  same  indi¬ 
cating  that  the  low  frequency  and  high  frequency  radiation  are  consistent  with  a 
simple  source.  Comparison  of  mb  and  Ms  for  the  Mammoth  Lakes  earthquake 
sequence  and  Nevada  Test  Site  explosions  show  that  all  of  the  earthquakes 
analyzed  discriminate  from  explosions.  Seismic  moment  ( M0 )  is  related  to  Mi 
and  Ms  by  the  following  relationships: 

log  M0  =  (1.20  ±  0.05)  Af),  +  (17.49  ±  0.19)  , 
log M0  -Ms  +  19.4 (±0.06)  . 

We  have  empirically  related  the  gross  logarithmic  amplitude  parameters  Mi, 
m b,  and  Ms,  and  the  seismic  spectrum  observed  at  near  epicentral  distance. 
One  representative  event  of  this  group  was  selected  as  a  standard  and  Mi,  m t, 
and  Ms  (at  their  observed  frequencies)  for  that  event  were  scaled  to  the 
observed  seismic  spectrum.  Corresponding  magnitudes  for  three  other  events 
were  scaled  to  the  observed  seismic  spectrum  using  the  derived  relationships. 
For  events  of  similar  depth,  the  local  spectral  values  can  be  predicted  from  the 
distant  magnitude  measurements  to  within  10%.  For  a  deeper  event  spectral 
values  were  predicted  to  within  about  25%.  This  confirms  that  the  spectral 
shape  inferred  from  distant  measurements  of  magnitude  is  congruent  with  the 
spectral  shape  observed  directly  in  the  near  field.  The  majority  of  Mammoth 
Lakes  events  examined  can  be  interpretated  as  simple  events.  In  comparison. 


several  of  the  San  Fernando  aftershocks  in  the  same  magnitude  range  studied 
by  Tucker  and  Brune  (1977)  could  be  interpreted  as  partial  stress  drop  events. 
The  near  equality  in  the  high  and  low  frequency  radiation  for  Mammoth  Lakes 
aftershocks  suggest  that  for  each  event  the  seismic  radiation  over  the  frequency 
band  studied  is  coming  from  1  simple  same  source,  possibly  the  breaking  of  a 
single  asperity.  The  suggested  differences  between  Mammoth  Lakes  aftershocks 
and  the  San  Fernando  earthquakes  may  arise  because  most  of  the  events  stu¬ 
died  for  Mammoth  Lakes  are  occurring  in  the  more  competent,  higher  stress, 
granitic  block  above  a  deeper  through-going  fault  plane.  The  data  suggest  the 
idea  that  partial  stress  drop  events  are  more  common  on  developed  fault  sys¬ 
tems,  whereas  full  stress  drop  events  are  more  common  in  the  surrounding 
country  rock  as  it  adjust  to  the  stress  field  generated  by  movement  on  the 
developed  fault  plane. 

1.0  Introduction 

Earthquakes  with  similar  local  magnitude,  M i,  presumably  have  similar 
short-period  amplitudes,  and  yet  they  can  show  large  differences  in  their  long- 
period  amplitudes.  Regional  variations  in  the  relative  excitation  of  high  and  low 
frequency  seismic  waves  have  been  noted  in  a  number  of  studies  in  western 
North  America.  Brune.  Espinosa  and  Oliver  (1963)  found  a  wide  range  in  the  sur¬ 
face  wave  excitation,  for  a  given  Mi ,  of  earthquakes  in  the  California-Nevada 
region.  Wyss  and  Brune  (1968)  found  that  earthquakes  in  Nevada  often  gen¬ 
erated  smaller  amplitude  surface  waves  for  a  given  Mi,  than  do  earthquakes 
along  the  San  Andreas  fault  or  along  the  northern  California  coast.  In  a  later 
study.  Wyss  and  Brune  (1971)  noted  large  differences  in  the  ratio  of  Gutenburg 
energy,  Eq,  determined  from  ML,  with  seismic  moment,  M0,  determined  from 
surface  waves,  for  earthquakes  in  the  California  region.  In  particular,  they 
found  that  earthquakes  along  the  Sierran  front  to  the  north  of  Bishop,  California 
had  small  ratios  of  Ms  to  Mi  when  compared  to  other  areas  in  California.  How¬ 
ever  at  that  time  it  was  not  possible  to  assess  the  effects  of  variations  in  source 
depth  and  source  mechanism.  In  contrast,  Tucker  and  Brune  (1977)  found  that 
some  of  the  larger  aftershocks  (Mi  >4)  of  the  San  Fernando  earthquake 
sequence  had  spectral  amplitudes  at  0.1  Hz  (based  on  surface  waves)  which  were 
as  much  as  an  order  of  magnitude  larger  than  the  spectral  amplitude  at  1  Hz 
(based  on  body  waves),  i.e.  the  spectrum  was  nt-c  flat  between  0.1  and  1.0  Hz. 

These  differences  can  arise  due  to  the  effects  of  the  path,  focal  depth,  radi¬ 
ation  factor,  and  source  parameters.  Previous  studies  (Wyss  and  Brune,  1971; 
Thatcher,  1972)  suggest  that  the  observed  differences  between  low  and  high  fre¬ 
quency  radiation  result  primarily  from  variations  in  the  source  spectrum. 
Brune  p.t  al  (1985)  have  speculated  that  three  classes  of  earthquakes  exist: 
"partial  stress  drop"  events  (Act<<  100  bars);  "full  stress  drop"  events  (Act*  100 
bars);  and  "overstress"  events  (Act>  300  bars),  where  the  relative  excitation  of 
high  and  low  frequency  radiation  varies  depending  on  the  type  of  event.  Full 
stress  drop  events  are  those  with  a  stress  drop  comparable  to  the  effective 
stress,  i.e.,  on  the  order  of  a  hundred  bars.  This  value  of  stress  agrees  with  esti¬ 
mates  of  the  upper  limit  of  shear  stress  across  the  San  Andreas  fault  based  on 


the  lack  of  an  observable  frictional  heat  flow  anomaly  (Brune  et  ai. ,  1969; 
Lachenbruck  and  Sass,  1980).  Partial  stress  drop  events  may  occur  when  the 
fault  locks  soon  after  the  rupture  passes  so  that  the  average  slip  over  the  fault 
cannot  reach  a  value  corresponding  to  the  initial  dynamic  stress  drop  over  the 
whole  fault.  Alternatively,  they  may  occur  when  the  stress  release  is  not  uni¬ 
form  and  coherent  over  the  whole  fault  plane,  but  rather  is  more  like  a  series  of 
multiple  events  with  part  of  the  fault  remaining  locked.  Such  complexities  could 
occur  as  a  result  of  asperities  on  the  fault,  or  due  to  a  complex  fault  plane 
geometry.  Overstress  events  may  result  possibly  from  a  high  quasi-static  stress 
buildup  due  to  strength  inhomogeneities  (barriers  or  asperities)  which  are 
stressed  during  an  earthquake,  or  directly  from  the  failure  process,  i.e.,  fault 
slip  could  cause  a  rapid,  quasi-static  stress  buildup  to  much  higher  levels  than 
could  be  maintained  by  the  rocks  on  the  time  scale  of  long-term  tectonic 
processes.  In  the  latter  case,  the  very  high  stress  drop  events  are  not  neces¬ 
sarily  indicative  of  regions  which  were  under  high  stress  before  the  earthquake. 

Understanding  the  differences  in  the  relative  excitation  of  high  frequency 
(0.5-10.0  Hz)  and  low  frequency  (0.05-0.125  Hz)  seismic  waves,  and  the  propaga¬ 
tion  of  this  energy  to  regional  distances,  is  important  for  understanding  source 
mechanisms.  It  is  particularly  relevant  to  the  problem  of  discriminating  small 
explosions  from  earthquakes  at  regional  distance.  The  mb-Ms  discriminant, 
which  compares  the  relative  excitation  of  body  versus  surface  waves,  is 
presently  one  of  the  primary  techniques  for  distinguishing  explosions  from 
earthquakes. 

This  study  examines  the  relative  excitation  of  high  and  low  frequency 
seismic  energy  from  earthquakes  of  the  1980  Mammoth  Lakes,  California, 
sequence.  These  earthquakes  have  been  studied  more  thoroughly  than  perhaps 
any  other  earthquake  sequence  in  the  western  United  States.  A  dense  network 
of  high-gain  seismographs  provide  data  for  accurate  hypocentral  coordinates 
and  first  motions  for  most  events.  Moreover.a  large  number  of  broad-band  digi¬ 
tal  and  strong-motion  recordings  have  been  collected  in  the  epicentral  region, 
which  is  surrounded  by  a  network  of  long-period  seismographs,  and  is  20°  to  30° 
from  the  Canadian  Network  high  gain,  short-period  stations.  The  events  are  par¬ 
ticularly  well  suited  for  evaluating  the  mb-Ms  seismic  discriminant  since  they 
are  in  a  environment  similar  to  the  Nevada  Test  Site,  approximately  150  km  to 
the  southeast. 

2.0  Summary  of  the  Mammoth  Lakes  Earthquake  Sequence 

The  Mammoth  Lakes  earthquake  sequence  and  the  geology  and  geophysics 
of  the  Long  Valley  Caldera  with  which  they  are  associated,  have  been  discussed 
in  numerous,  recent  papers  [See.  for  example,  the  special  issue  of  the  Journal  of 
Geophysical  Research,  November,  1985],  The  earthquake  sequence  began  in 
October,  1979,  peaked  in  May,  1980  with  the  occurrence  of  four  ML  >  6  earth¬ 
quakes  over  a  two  day  period,  and  has  continued  to  the  present.  These  earth¬ 
quakes  are  associated  with  the  intersection  of  the  Sierra  Nevada  frontal  fault 
system  and  the  Long  Valley  Caldera.  The  Long  Valley  Caldera  was  formed  0.7 
million  years  ago  by  collapse  and  subsidence  associated  with  the  eruption  of  the 


Bishop  tuff  (Bailey  et  al.,  1976),  and  volcanic  activity  has  continued  on  a  reduced 
scale  into  the  Holocene.  The  dominant  fault  of  the  area  is  the  Hilton  Creek  fault 
along  which  several  hundred  meters  of  pure  normal  fault  displacement  has 
occurred  since  the  formation  of  the  Long  Valley  Caldera  (Bailey  et  al.,  1976). 
Extension  of  the  Hilton  Creek  fault  into  the  Long  Valley  Caldera  appears  to  have 
occurred  as  recently  as  0.3  million  years  ago,  possibly  indicating  that  the  cal¬ 
dera  had  cooled  sufficiently  by  that  time  to  support  stresses  large  enough  to 
generate  earthquakes. 

Lide  and  Ryall  (1985)  discuss  the  distribution  of  344  master  event  aft¬ 
ershock  locations  in  the  weeks  following  the  four  Mi  >  6  events  (the  time  period 
for  the  events  examined  in  this  study).  Most  of  the  earthquakes  they  located 
occurred  in  the  Sierrian  block,  south  of  the  caldera.  In  map  view,  they  lay  in  a 
triangular  shaped  zone,  of  which  the  north  boundary  is  approximately  formed  by 
the  southern  rim  of  the  caldera  and  of  which  the  east  boundary  is  approximately 
formed  by  the  Sierran  frontal  fault  system.  Three  trends  can  be  recognized  in 
the  aftershock  data:  two  NNE  trends,  one  of  which  forms  the  eastern  boundary 
of  the  aftershock  sequence,  and  a  WNW  trend  which  lies  along  the  southern  boun¬ 
dary  of  the  caldera.  A  cross-section  perpendicular  to  the  NNE  trends  shows  the 
aftershocks  in  the  Sierrian  block  occur  in  a  triangular  shaped  zone  the  base  of 
which  dips  at  about  45°  to  the  east.  Thus  in  the  west,  events  extend  down  to  4  -  5 
kms  depth,  while  in  the  east  they  extend  down  to  11  -  13  kms  depth.  The  Mi  >6 
events  lie  near  the  base  of  this  zone. 

3.0  Data 

This  study  utilizes  seismograms  of  those  earthquakes  whose  epicenters  are 
shown  in  Figure  1  and  listed  in  Table  1.  These  events  comprise  only  a  small  frac¬ 
tion  of  the  total  number  of  events  which  have  occurred  in  the  Mammoth  Lakes 
area  since  the  earthquake  sequence  began  in  October,  1978.  Events  selected  for 
study  consist  of  those  for  which  strong-motion  or  digital  local  seismograms, 
long-period  regional  seismograms,  and  short-period  Canadian  Network  body- 
wave  magnitudes  where  available.  All  events  have  been  located  using  arrival 
times  from  the  local  stations  shown  in  Figure  1,  with  additional  readings  from 
local  stations  operated  by  the  University  of  Nevada  and  the  California  Institute 
of  Technology.  Vetter  (1984)  has  determined  first  motion  focal  mechanisms  for 
many  of  these  earthquakes. 

Seismograms  used  in  this  study  are  from  the  seismograph  stations  shown  in 
Figure  1  and  listed  in  Table  2.  The  SMA-1  strong-motion  accelerographs  were 
installed  by  the  California  Division  of  Mines  and  Geology  (Turpin,  1980),  or  by  the 
University  of  Southern  California  (Moslem  e£  al.,  1983).  Digital  stations  operated 
by  the  U.S.  Geological  Survey  and  described  by  Archuleta  et  al.  (1982),  consist 
of  either  force-balance  accelerometers  or  velocity  transducers  and  Spreng- 
nether  DR-100  recorders.  The  SMA-1  analog  seismograms  were  digitized  at  a 
sample  interval  of  either  0.025  sec  or  0.005  sec;  the  digital  recordings  supplied 
from  the  U.S.  Geological  Survey  were  digitized  at  an  interval  of  0.005  sec.  All  the 
above  stations  are  denoted  by  triangles  in  Figure  1.  For  events  along  the  south¬ 
ern  boundary  of  the  caldera,  the  SMA-1  and  digital  seismographs  sample  a  large 


fraction  of  the  focal  sphere.  The  radiation  pattern  of  events  south  of  the  caldera 
is  not  as  well  sampled. 

Figure  1  also  shows  the  locations  of  the  regional  long-period  stations 
(denoted  by  squares)  whose  seismograms  were  used  to  determine  surface  wave 
magnitudes.  These  stations  surround  the  Mammoth  Lakes  region  and,  with  the 
exception  of  the  California  Institute  of  Technology  station  at  Pasadena  (PAS), 
belong  to  the  world  wide  network  (WWSSN).  Readings  from  the  Canadian  Net¬ 
work  stations  (denoted  by  stars)  were  used  to  determine  body  wave  magnitude. 
All  events  in  this  study  have  nearly  identical  paths  to  the  regional  and  telese- 
ismic  stations,  so  the  data  are  uncontaminated  by  path  effects. 

4.0  Data  Analysis 

4.1  Magnitude  Analysis 

Mi'.  Chavez  and  Priestley  (1985)  determined  Mi  for  82  events  of  the  1980  Mam¬ 
moth  Lakes  earthquake  sequence  (Table  2,  column  7)  by  converting  the  digital, 
horizontal  component  data  into  equivalent  Wood-Anderson  seismograms.  In 
determining  Mi  for  these  events,  they  developed  and  used  an  attenuation  func¬ 
tion  appropriate  to  the  Great  Basin.  "True”  Mi  based  on  Wood-Anderson  ampli¬ 
tude  measurements  for  the  main  shocks  and  larger  aftershocks  of  the  sequence 
have  been  reported  by  the  University  of  California-Berkeley  and  by  the  Califor¬ 
nia  Institute  of  Technology.  Chavez  and  Priestley  found  that  their  magnitudes 
were  in  close  agreement  with  Mi  determined  from  the  Wood-Anderson  seismo¬ 
grams  recorded  in  California. 

mft:  The  Earthquake  Data  Reports  of  the  U.S.  Geological  Survey  report  peak-to- 
peak  amplitude,  period,  and  individual  station  body-wave  magnitudes.  From  this 
listing  we  have  tabulated  the  body-wave  magnitudes  reported  for  Canadian  Net¬ 
work  stations  at  distances  greater  than  20°  from  Mammoth  Lakes,  and  from 
these  have  determined  an  average  mb  (Table  2,  column  8)  for  21  earthquakes 
studied. 

Ms:  We  have  determined  surface  wave  magnitude  for  34  events  of  the  Mammoth 
Lakes  sequence  using  amplitude  measurements  from  long-period  seismograms 
at  regional  distances.  Ms  (Table  2,  column  9)  was  determined  using  the 
definition  of  Marshall  and  Basham  (1972); 

Ms  =  logA  +B(A)  +  P(T)  (1) 

where  A  is  one-half  the  peak-to-peak  ground  amplitude  of  the  Rayleigh  wave  with 
a  period  T  at  a  distance  A,  B  (A)  is  a  correction  for  the  average  effects  of 
attenuation,  scattering,  geometrical  spreading  and  refraction,  and  P(T)  is  a 
correction  factor  for  dispersion  which  allows  measuring  the  surface  wave  ampli¬ 
tude  at  periods  other  than  20  seconds. 
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4.2  Body- wave  analysis 

Nature  of  the  Seismograms  at  Near  Distance:  The  seismograms  analyzed  display 
varying  degrees  of  complexity,  ranging  from  simple  pulses  to  wave  trains  many 
seconds  long.  This  variation  can  be  due  to  source  complexity,  position  of  the 
station  on  the  focal  sphere,  or  to  propagation  affects.  Some  of  the  complex 
seismograms  observed  are  easily  attributable  to  source  complexity  since  in 
those  cases  seimograms  with  both  complex  and  simple  nature  have  been 
observed  at  the  same  station  for  nearly  co-located  events.  In  other  cases,  how¬ 
ever,  both  complex  and  simple  seismograms  have  been  observed  at  different 
stations  for  the  same  event.  In  these  instances,  seismograms  showing  large, 
simple  S-waves  often  have  relatively  small  P-wave  amplitudes,  and  those  with 
complex  S-wave  have  large  P-wave  arrivals,  indicating  that  the  complexity  is 
related  to  the  position  of  the  station  on  the  focal  sphere,  i.e.,  if  the  S-wave 
amplitude  is  large  (away  from  a  node),  it  will  exceed  the  scattered  energy  and 
the  resulting  seismogram  will  have  a  simpler  character,  whereas  if  the  station  is 
near  an  S-wave  node,  the  scattered  energy  will  predominate  and  the  seismogram 
will  appear  to  be  more  complex  and  have  a  lower  amplitude  (Tucker  and  Brune, 
1977).  In  some  cases,  it  is  apparent  from  the  long  period  oscillations  following 
the  S-wave,  that  earth  structure  contributes  to  the  complexity  of  the  seismo¬ 
grams.  The  effects  of  earth  structure  on  seismic  propagation  will  be  discussed 
in  more  detail  below. 

Seismic  Spectra:  We  have  determined  the  seismic  spectra  for  the  events  indi¬ 
cated  in  Table  1.  We  have  compared  our  spectral  results  to  those  presented  by 
Archuleta  et  at.,  and  find  good  agreement  with  our  results,  which  also  include 
the  analysis  of  the  strong-motion  data.  Spectral  parameters  for  the  smaller 
events  in  Table  4  of  Archuleta  et  at.  are  consistent  with  the  spectral  parameter 
determined  by  Andrews  (  written  communications,  1985).  In  the  discussion 
below  we  have  used  the  spectral  parameters  for  small  Mammoth  Lakes  aft¬ 
ershocks  reported  by  Archuleta  ef  al.  (1982). 

Displacement  amplitude  specta  for  the  events  indicated  in  Table  1  were 
computed  from  the  locally  recorded  body-wave  data.  Time  windows  were  chosen 
to  be  long  enough  to  include  the  entire  S-wave  arrival,  but  short  enough  to 
minimize  the  inclusion  of  surface  wave  energy.  Rather  than  rotate  the  data  into 
radial  and  transverse  components,  we  computed  the  vector  sum  of  the  two  hor¬ 
izontal  component  spectra.  The  spectra  were  normalized  to  a  hypocentral  dis¬ 
tance  of  10  km,  and  plotted  as  log  spectral  amplitude  (fl)  vs.  log  frequency  (cj) 
and  interpreted  in  terms  of  a  flat,  low-frequency  asymptote,  fi0,  a  high  fre¬ 
quency  fall-off  with  a  slope  in  the  range  -1  to  -3,  and  the  corner  frequency  (/<-) 
at  which  the  two  trends  intersect.  In  almost  all  cases,  there  is  a  second  corner 
in  the  data  at  frequencies  between  10  and  16  Hz.  Below  this  frequency  the  spec¬ 
tra  drop  off  with  a  slope  of  -4  to  -5.  The  asymptotes  were  fit  by  eye  using  an 
interactive  computer  routine;  an  example  of  the  fiting  is  shown  in  Figure  2. 
Although  simple,  the  fitting  technique  used  here  is  felt  to  be  adequate,  given 
other  known  errors  involved  in  estimating  fl0  and  /c.  Where  we  have  events  in 
common,  our  spectral  fits  are  consistent  with  those  of  Andrews  (written  com¬ 
munication,  1985),  who  used  a  more  objective  technique. 


The  fits  were  then  corrected  for  average  radiation  pattern  and  converted  to 
full-space  spectra  by  accounting  for  the  free  surface  effect,  using  the  following 
eauation: 


Qa(fulL-space)l0km  = 


Qa(obseTved)iCkm 


k  R 


(2) 


where  k  is  the  correction  for  the  free  surface  SH-wave  amplication,  and  R$9  is 


the  correction  for  the  average  radiation  pattern.  Values  of  Q0  and  fc ,  and  their 
associated  errors  are  given  in  columns  10  through  13  of  Table  2. 

Estimation  of  Data  Errors-.  Where  more  than  one  recording  was  available  for  an 
event,  errors  in  H0  and  fc  due  to  scatter  in  the  data  were  estimated  by  calculat¬ 
ing  an  average  spectral  value  and  a  standard  deviation.  The  average  values  of 


n0,  <fi0>,  were  determined  by  the  formula 


<n0  > 
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=  antilog  —  £  log  n0. 


where  Q0.  is  the  long-period  spectral  level  at  the  i ^  station  ,  and  NS  is  the 
number  of  stations  recording  the  event.  The  standard  deviation  of  the  log  aver¬ 
age  value  of  n0  was  determined  by  the  formula 


s.  d.  (log  <  n0  > )  =  j-  2  [l°s  ( no )  -  log  <  nD  >  2  J 
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A  multiplicative  error,  EQ0  factor  is  then  determined  from 


EQ0  =  antilog 


s.d.  ( log  <  n0  >  ) 


</c>  and  Efc  were  computed  in  a  similar  manner.  As  pointed  out  by  Archuleta 
et  al.  (1982),  when  calculating  average  values  of  Q0  or  fc,  it  is  necessary  to  com¬ 
pute  the  averages  using  the  equations  in  the  above  forms  to  give  equal  weight  to 
each  observation.  If  a  simple  arithmetical  average  value  is  determined  it  will  be 
biased  toward  larger  values,  since  the  errors  associated  with  D0  and  fc  are  log- 
normally  distributed. 

Sources  of  error  in  estimating  the  seismic  spectrum-.  The  error  bars  calculated 
for  the  long  period  level  D0  and  the  corner  frequency  fc  in  the  above  section 
represent  only  the  variations  in  Da  and  fe  due  to  scatter  in  the  observed  record¬ 
ings.  There  are,  however,  a  number  of  sources  of  systematic  error  inherent  in 
our  estimates  of  ff0  and  fc.  These  arise  from  both  wave  propagation  and  record¬ 
ing  effects,  as  well  as  the  assumptions  and  procedures  made  in  analysing  the 
data. 

The  correction  of  the  S-wave  amplitude  for  reflection  at  the  free  surface 
assumes  that  the  waves  are  horizontally  polarized  (SH),  however  since  we  have 


not  rotated  the  data  into  transverse  and  radia'  components,  the  energy  is  in 
general  some  mixture  of  SH  and  SV-motion.  Abe  (1974)  discusses  experimental 
studies  of  S-wave  amplification  on  reflection  and  concludes  that  using  k  =  2  in 
equation  (2)  (which  we  have  done)  applies  approximately  to  SH-wave  signals  con¬ 
taining  some  SV-motion.  The  station-source  configuration  shown  in  Figure  1 
indicates  that  multiple  arrivials  due  to  reflection  should  not  be  a  problem  in  this 
study,  since  the  S-waves  are  recorded  at  close  epicentral  distance. 

No  Q  correction  was  made  in  determining  the  spectra.  Archuleta  et  al. 
(1982)  found  that  variations  in  whole-path  Q  had  a  small  effect  on  the  measure¬ 
ment  of  n0  for  Mammoth  Lakes  events.  More  recent  studies,  however,  suggest 
that  near  source  recordings  of  the  Mammoth  Lakes  aftershocks  (Peppin,  1985: 
Barker,  1985),  and  possibly  near  source  recordings  in  general  (Anderson,  1985: 
Cranswick  et  al,  1985)  are  strongly  contaminated  by  variations  in  the  structure 
very  near  to  the  recording  site.  Archuleta  (1985)  has  analysed  data  recorded  on 
a  borehole  seismograph  in  the  epicentral  region  of  the  Mammoth  Lakes  aft¬ 
ershocks,  to  estimate  the  effect  of  the  near  surface  structure  on  the  seismic 
recording  of  the  data  (and  consequently  on  the  seismic  spectra).  In  this  experi¬ 
ment,  recordings  were  made  at  three  levels  in  the  borehole  --  at  the  surface, 
just  below  the  sediment-bedrock  interface  (35  meters  depth),  and  well  into  the 
bedrock  (166  meters  depth).  The  most  obvious  difference  between  the  surface 
and  bedrock  time  series  was  the  increased  complexity  of  the  surface  recordings 
compared  with  those  in  bedrock.  Archuleta  found  little  difference  in  the  spec¬ 
tral  levels  of  the  two  bedrock  recordings,  but  approximately  a  four-fold  ampli¬ 
tude  increase  in  the  data  from  the  surface  instrument.  This  amplication  can 
partially  be  accounted  for  by  the  free  surface  amplication  factor  due  to 
reflection  (normally  taken  as  a  factor  of  2  for  pure  SH-waves),  and  partially  due 
to  the  effects  of  the  impedance  contrast  between  the  glacial  moraine  at  the  sur¬ 
face  and  the  underlying  bedrock  (approximately  a  factor  of  1.85).  Of  relavance 
to  the  present  study,  spectral  ratios  between  the  surface  and  bedrock  instru¬ 
ments  show  that  the  sediment  amplification  is  essentially  uniform  across  the  2.0 
to  50.0  Hz  band  and  does  not  shift  fc. 

The  effect  of  amplication  due  to  the  impedance  contrast  between  the  near 
surface  sediment  and  the  underlying  bedrock  at  the  various  sites  is  somewhat 
more  difficult  to  acess,  as  there  is  a  variety  of  sediment  types  and  thickness 
over  the  array  shown  in  Figure  1.  On  the  average,  the  effect  of  the  sediment- 
bedrock  impedance  contrast  should  be  approximately  a  factor  of  2,  irrespective 
of  the  variation  of  the  sediment  thickness  observed  in  the  Mammoth  Lakes 
region.  Thus  based  on  Archuleta’s  borehold  results,  the  long  period  spectral  lev¬ 
els  obtained  from  routine  analysis  of  the  seismic  spectra  of  locally  recorded  sur¬ 
face  data  in  the  Mammoth  Lakes  region  are  probably  overestimated  by  a  factor 
which  is  close  to  two  but  not  by  an  order  of  magnitude  or  more  (Barker,  1985; 
Peppin,  1985).  However,  there  appear  to  be  only  negligible  effects  of  near  sur¬ 
face  structure  on  the  estimated  corner  frequencies. 

The  spectra  have  been  normalized  to  a  hypocentral  distance  of  10  km.  Not 
all  events  were  located,  and  in  those  cases  the  S-P  times  were  used  to  determine 
distance.  The  hypocentral  accuracy  of  the  located  events  is  about  ±  3  km 
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(Cockerham,  personal  communication)  which  would  cause  as  much  as  30%  error 
in  the  normalization  of  any  single  spectra.  Systematic  errors  in  hypocentral 
location  could  affect  the  normalized  average  spectra,  however  there  is  no  reason 
to  suspect  systematic  location  errors. 

A  final  source  of  error  inherent  in  the  estimate  of  no  comes  from  the 
assumption  of  a  single,  average  radiation  pattern  correction.  If  sufficient  instru¬ 
ments  were  available,  or  if  the  source  parameters  (including  direction  and  velo¬ 
city  of  rupture)  were  known,  this  error  could  be  reduced.  A  value  of  0.6  was 
used  for  R ^  (Thatcher  and  Hanks,  1973)  in  equation  2.  At  some  azimuths,  the 
amplitude  could  be  1.67  times  larger  than  the  amplitude  determined  assuming 
this  correction,  while  near  a  node  in  the  radiation  pattern  the  amplitude  could 
be  quite  small.  Furthermore,  if  the  source  dimension  is  comparable  to  the 
wavelength  involved,  source  propagation  might  distort  the  radiation  pattern  by 
focusing  or  defocusing  energy,  thereby  changing  the  apparent  corner  frequency 
for  different  directions  and  possibly  reducing  the  calculated  moment. 

4.3  Surface  Wave  Analysis 

Determination  Surface  Wave  Moment:  Surface  waves  provide  estimates  of  fi0 
which  sample  a  lower  frequency  portion  of  the  seismic  spectrum,  and  which  are 
independent  of  those  obtained  from  S-wave  spectra.  We  have  determined  sur¬ 
face  wave  moments  for  the  earthquakes  by  fitting  the  general  features  of  the 
observed  Rayleigh  waves  recorded  at  regional  distances,  with  a  synthetic  surface 
wave  seismograms  (Figure  3).  Parameters  needed  to  calculate  the  synthetic 
seismograms  are  the  earth  velocity  and  attenuation  structures  and  the  fault 
parameters  of  the  individual  earthquakes.  Since  the  events  are  located  at  the 
western  edge  of  the  Great  Basin,  while  most  of  the  seismograph  stations  are 
located  in  western  North  America,  we  used  the  velocity  structure  determined  for 
the  Great  Basin  from  surface  wave  dispersion  measurements  (Priestley  and 
Brune,  1978).  The  Great  Basin  attenuation  model  by  Patton  and  Taylor  (1984) 
was  used,  with  a  slight  lowering  of  the  surficial  shear  wave  Q.  The  earth  model 
used  in  making  the  surface  wave  calculations  is  given  in  Table  3.  Since  for  the 
majority  of  events  we  have  studied,  we  do  not  have  fault-plane  solutions,  nor  do 
we  know  the  focal  depth  we  determined  an  average  relationship  for  the  Mam¬ 
moth  Lakes  area  in  the  following  manner.  Vetter  (1984)  has  determined  fault 
plane  sjlutions  for  a  large  number  of  events  in  the  Mammoth  Lakes  area,  and 
found  that  the  mechanisms  can  generally  be  subdivided  into  two  groups.  The 
first  group  consists  of  earthquakes  with  strike-slip  mechanisms  (average  plane 
of  strike  N14°E,  dip  80°ESE,  slip  N16°E  or  strike  N73°W,  dip  80°SSW,  slip  N76°E) 
which  are  most  commonly  at  depths  less  than  9  km.  The  second  group  consists 
of  earthquakes  with  oblique  or  normal  mechanisms  (average  planes  of  strike 
N03°E,  dip  64°E,  slip  N40°E  or  strike  N50°W,  dip  38°SW,  and  slip  N86°W)  usually  at 
depths  greater  than  9  km.  Synthetic  surface  wave  seismograms  were  calculated 
for  both  mechanisms,  for  2,  5,  8,  11  and  14  km  depth.  Surface-wave  amplitudes 
were  determined  from  the  10  synthetic  records  at  each  station  and  an  average 
conversion  factor  determined  for  both  focal  mechanisms  and  all  focal  depths. 
The  resulting  relationship  between  surface  wave  magnitude  Ms  and  seismic 
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moment  M0  for  the  Mammoth  Lakes  area  is  given  by  the  equation 

logM0  =  Ms  +  19.40(±0.06).  (3) 

Sources  of  Error:  Tucker  and  Brune  (1977)  used  surface  wave  spectra  to  make 
lower  frequency  estimates  of  seismic  moment  for  comparison  with  moments 
determined  from  body  waves.  They  made  a  careful  estimate  of  the  errors 
involved  in  estimating  the  seismic  moment  from  the  surface  wave  spectra. 
These  errors  included  procedural  errors,  path  structure  and  attenuation,  lateral 
refraction  (multipathing),  source  finiteness,  and  uncertainties  in  fault  parame¬ 
ters  (depth,  dip,  strike,  and  slip  vector)  in  correcting  for  the  radiation  pattern. 
They  concluded  that  these  effects  could  lead  to  a  maximum  error  in  the  final 
result  of  no  more  than  a  factor  of  4. 

Low-frequency  Spectral  Level  from  Surface  Waves:  Equation  (3)  allows  us  to 
rapidly  determine  the  seismic  moment.  M0,  from  the  surface  wave  amplitude 
measurements.  These  surface  wave  moments  were  converted  to  full-space  spec¬ 
tral  levels  (QJ  at  their  observed  frequencies  for  comparison  with  the  spectral 
levels  determined  from  the,  near  source  body-wave  recordings.  Using  a  shear 
wave  speed  (/?!  of  3.5  fcs-1,  a  density  (p)  of  2.9  kg  m-a,  0.6  for  the  radiation  fac¬ 
tor  (R$v),  and  the  seismic  moment  determined  by  fitting  the  surface  waves,  the 
spectral  level  from  the 

surface  wave  moments,  normalized  to  10  km  is  given  by 

n0 (full -space) ,ojan  =  =  3.84»  .  (4) 

This  gives  values  of  H0  at  surface  wave  frequencies  for  comparison  with  the  spec¬ 
tra  obtained  from  the  body-wave  accelerograms.  Values  of  fl0  derived  from  the 
surface  wave  moments  and  the  associated  error  in  fi0  are  given  in  column  14 
and  15  of  Table  1.  The  frequency  at  which  the  surface  waves  fi0  was  determined 
and  its  associated  error  are  given  in  column  16  and  17. 

5.0  Results 

5.1  Comparison  of  short-period  and  long-period  estimates  of  seismic  moment 

Figure  4  is  a  plot  of  no  determined  from  the  spectral  analysis  of  the  strong- 
motion  and  digital  data  recorded  at  near  distances  (A  <20  km),  and  H0  deter¬ 
mined  from  surface  waves  recorded  at  larger  distances  (A  >375  km).  Error  bars 
indicate  the  scatter  in  the  observations.  The  reference  line  through  the  data  is 
of  unit  slope  and  zero  intercept.  'ITiose  points  lying  above  the  line  have  larger 
local  body  wave  fl0  relative  to  the  observed  surface  wave  fl0,  while  points  below 
the  line  have  larger  surface  wave  Q0  compared  to  the  observed  local  body  wave 
n„.  As  shown  in  this  figure,  most  of  the  events  we  have  studied  show  no 
significant  differences  in  H0  measured  at  frequencies  near  1  Hz  compared  to  Q0 
measured  at  frequencies  near  0.1  Hz. 


Differences  in  Q„  over  this  band  can  arise  due  to  partial  stress  drop  (Brune, 
1970)  or  to  afterslip  on  the  fault  (Hartzell  and  Brune,  1979).  The  agreement 
between  the  low  and  high  frequency  determinations  of  00  indicates  that  the 
Mammoth  Lakes  events  we  have  studied  are  simple,  that  is  they  are  neither  par¬ 
tial  stress  drop  events  nor  events  with  a  large  long-period  heave  following  the 
breaking  of  the  initial  asperity.  Thus  our  data  suggest  that  seismic  radiation 
across  the  band  studied  is  produced  from  a  simple  mechanism,  possibly  the 
breaking  of  a  single  asperity.  This  may  arise  due  to  a  higher  overall  state  of 
stress  in  the  crust  in  the  vicinity  of  Mammoth  Lakes  as  suggested  by  Wyss  and 
Brune  (1971);  most  of  the  aftershocks  we  have  studied  occur  in  the  granitic 
block  where  strength  and  stress  may  be  higher  than  on  the  deeper,  through- 
going  fault  plane.  Alternatively,  the  granitic  block  may  be  now  so  brecciated 
that  stress  can  only  be  concentrated  over  planes  of  limited  extent  which  break 
as  a  single  unit. 

5.2  Relationship  between  m* ,  and  M, 

Figure  5  compares  the  mb  vs.  Ms  data  for  the  Mammoth  Lakes  earthquakes 
to  the  North  American  earthquake  and  explosion  data  of  Marshall  and  Basham 
(1972).  Magnitudes  in  that  study  were  determined  from  seismograms  recorded 
on  Canadian  Network  stations.  In  our  study,  Canadian  Network  data  were  used 
to  determine  mb,  while  data  from  long-period  stations  in  the  western  U.S.  were 
used  to  determine  Ms.  Though  it  is  possible  that  this  might  lead  to  a  systematic 
difference  between  Ms  determined  by  Marshall  and  Basham,  and  the  Ms  deter¬ 
mined  for  the  Mammoth  Lakes  earthquakes  by  us,  there  does  not  seem  to  be  a 
strong  effect.  This  is  evidenced  by  the  fact  that  our  paths  overlap  those  of 
Marshall  and  Basham  to  a  large  extent,  and  the  data  sets  in  Figure  5  overlap  as 
well.  Neither  the  Marshall-Basham  data  nor  the  Mammoth  Lanes  earthquakes 
data  have  been  corrected  for  depth  variations.  All  Mammoth  Lakes  events  fall 
within  the  earthquake  population  and  discriminate  from  explosions,  although 
some  events  are  offset  from  the  earthquake  mean  towards  the  explosion  popula¬ 
tion. 

5.3  Magnitude-Moment  relationships 

Log  M0  vs.  M^,  relationship'.  Figure  6  is  a  plot  of  log  seismic  moment,  M0,  derived 
from  the  spectra  of  the  near-source  digital  and  strong-motion  recordings,  vs.  Mi 
values  published  by  Chavez  and  Priestley  (1985).  They  found  that  the  data  are 
fit  by  the  line 

log-M0  =  (1.20  ±0.05)  Ml  +(17. 49  ±0.19)  ( 1  <  Aft  <  6)  .  (5) 

and  do  not  exhibit  curvature  of  the  type  noted  by  Bakun  (1984)  or  Hanks  and 
Boore  (1984)  for  other  California  earthquake  sequences. 

The  log  M0  vs.  Mi  relationship  given  in  equation  (5)  differs  from  that  given 
by  Archuleta  et  al.,  (1982)  which  used  data  primarily  in  the  range  3SA//,  §5  with 
only  one  event  above  Mi  5.  Our  data  span  the  range  1  £  Mi  §6  with  seven  events 
above  Mi  5.  Archuleta  et  al.,  (1982)  found  that  the  log  M0  vs  Mi  relationship  for 


the  40  events  for  which  they  had  University  of  California  -  Berkley  Mi  determina¬ 
tions  was 


iogM0  =  (0.96 ±0.06) ML  +  (18.14 ±  0.23)  (2.9 ZML  S8.2) 

but  for  the  larger  events  with  more  reliable  Mi  estimates  it  was 

log  M0  =(1.05  ±0.08)  Mi  +  (17.76*0.33)  (3.5  S  ML  S  6.2)  . 

The  greater  slope  in  equation  (5)  is  due  to  the  fact  that  we  used  larger  moments 
for  events  with  Mi  >5.0  than  did  Archuleta  et  al.,  (1982).  For  example,  for  the 
May  27,  1980  event  at  14:50,  Archuleta  et  al.,  (1982)  found  a  seismic  moment  of 
2.33xl024  dyne-cm  based  on  two  recordings  whereas  Priestley  et  al.,  (1985) 
obtained  a  value  of  7.24xl024  dyne-cm  from  nine  recordings.  The  slope  of  equa¬ 
tion  (5)  is  less  than  Archuleta  et  al.,  (1982)  found  for  the  data  of  Thatcher  and 
Hanks  (1973),  but  is  identical  to  the  relation  of  Bakun  and  Lindh  (1977)  for  the 
Oroville,  California  earthquake  sequence: 

logM0  =  (1.21  ±  0.03) if*  +  (17.02  ±  0.07)  . 

Oroville,  in  the  western  Sierra  Nevada,  is  locatec  ,n  a  structural  environment 
which  is  more  similar  to  that  of  Mammoth  Lakes  than  the  southern  California 
region  studied  by  Thatcher  and  Hanks  (1973). 

In  a  study  of  five  different  source  regions  (Parkfield,  San  Juan  Bautista,  the 
Sargent  fault.  Coyote  Lake,  and  the  Livermore  Valley),  Bakun  (1984)  detected  an 
inflection  in  the  \ogM0  vs.  Mi  plot  at  about  Mi  3.  For  small  events  he  found  the 
relationship  to  be 


\ogM0  =  l.Z5ML  +  17  (1.5  &ML  Z  3.5) 

which  is  identical  to  our  equation  (5)  which  is  based  on  data  in  the  range 
1§A//,  g  6.  However,  for  events  in  the  same  magnitude  range  as  our  data  Bakun 
(1984)  found  the  relationship  to  be 

logA/0  =  1.5ML  +  16  (3.5  £  ML  g  6.25)  . 

Bakun  (1984)  and  Hanks  and  Boore  (1984)  have  noted  that  log  M0  vs.  Mi 
plots  for  California  earthquakes  exhibit  positive  curvature  when  the  range  of 
magnitudes  considered  is  sufficently  large.  Hanks  and  Boore  (1984)  computed 
theoretical  log  M0  vs.  Mi  relationship  based  on  the  assumption  of  a  constant 
stress  drop  of  a  100  bars.  They  thus  restricted  the  spectrum  scaling  law  to  be  a 
function  of  seismic  moment  alone.  Their  calculated  values  are  shown  by  the 
stars  in  figure  G.  They  show  only  slight  curvature  for  Mi  <  5  but  seismic 
moment  increasing  ever  more  rapidly  for  larger  magnitudes.  This  phenomenon 
is  not  observed  in  our  Mammoth  Lakes  data,  suggesting  that  assumption  of  con¬ 
stant  stress  drop  independent  of  magnitude  is  not  valid  for  this  earthquakes 
sequence.  However,  given  the  small  curvature  in  the  theoretical  values  below 
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Mik5,  together  with  the  scatter  in  the  observations,  the  data  are  consistent 
with  constant  stress  drop  for  the  smaller  events,  but  suggest  that  the  larger 
events  are  greater  than  100  bars  stress  drop,  This  does  not  necessarily  imply 
that  the  overall  stress  drop  of  the  larger  events  is  greater  than  100  bars,  but 
that  the  high  frequency  radiation  giving  rise  to  the  observed  Mi  may  be  coming 
from  the  breaking  of  a  high  stress  asperity.  Alternatively,  the  moments  from 
the  near-in  body-wave  data  may  be  giving  an  underestimate  of  the  total  seismic 
moment.  Given  et  al.  (1982)  found  from  long  period  teleseismic  data  that  the 
seismic  moment  of  the  801481450  event  was  1.1  x  1025,  a  factor  of  5  to  10  greater 
than  the  estimated  moment  from  the  near-in  body-wave  data. 

Log  M0  vs.  Ms  relationship:  Equation  (3)  is  the  log  M0  vs.  Ms  relationship  derived 
from  fitting  the  observed  long-period  surface  wave  seismograms  recorded  at 
regional  distances,  with  synthetic  surface  wave  seismograms  computed  for  the 
Great  Basin  structure  (Priestley  and  Brune,  1978;  Taylor  and  Patton,  1984).  We 
can  compare  equation  (3)  with  similar  relations  established  for  the  western 
United  States.  Wyss  and  Brune  (1968)  originally  established  a  moment-local 
magnitude  relationship  for  events  less  than  magnitude  6  in  the  western  U.S. 
Seismic  moments  were  determined  for  272  events  which  averaged  over  a  wide 
range  of  tectonic  regions  (California,  Nevada,  Arizona,  Utah  and  Baja  California), 
and  were  found  to  be  best  fit  by  the  curve  given  by  the  equation 

log!l/0  =  1.71^  +  15.1  (3<Ml<6). 

Using  the  Ms  vs.  ML  relationship  given  by  Wyss  and  Brune  (1968),  this  becomes  in 
terms  of  MS{C)  (the  surface  wave  magnitude  defined  by  Gutenberg  (1945)) 

log  M0  =  M$(c)  +  19.2  (3<A/s-<8). 

This  definition  can  be  compared  with  the  relationship  derived  directly  from 
Gutenberg’s  (1945)  definition  of  surface  wave  magnitude.  According  to  that 
definition,  a  magnitude  6  earthquake  produces  a  far-field  displacement  of 
100pm  at  distance  22°  for  surface  waves  of  period  20  seconds.  In  terms  of 
moment  Gutenburg’s  relationship  becomes 

log M0  -  Ms{g)  +  19-3  . 

Thus  this  point  represents  an  average  of  the  numerous  observations  on  which 
the  surface  wave  magnitude  was  based.  As  pointed  out  by  Richter  (1958),  the 
scale  was  adjusted  to  agree  with  the  local  magnitude  values  of  6  to  7.  Tucker 
and  Brune  (1977)  determined  M0  from  spectral  analysis  of  regionally  recorded 
surface  wave  seismograms  of  the  San  Fernando  aftershock  sequence.  They  also 
used  the  Marshall  -  Basham  definition  (equation  1)  to  determine  Ms  for  these 
events.  The  values  they  obtained  were  related  by  the  equation 

logM0  -Ms(xB}  +  19.3  , 
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identical  to  the  relationship  postulated  using  the  Gutenburg  definition  of  surface 
wave  magnitude.  Comparing  these  various  relationships  between  Ms  and  log  M0 
we  conclude  that  the  use  of  the  Great  Basin  velocity  and  attenuation  models 
(Priestley  and  Brune,  1978;  Patton  and  Taylor,  1984)  has  increased  the 
estimated  moment  of  the  events  by  26%  over  that  determined  for  the  Gutenburg 
continental  model. 

6.0  Comparison  of  magnitude  estimates  with  source  spectra 

We  empirically  determined  formulae  which  relate  the  various  types  of  mag¬ 
nitude  estimates  to  the  source  spectrum  by  computing  constants,  C,  which 
satisfy 


log  0  =  M  -  C 

where  Q  is  the  spectral  amplitude  normalized  to  10  km  at  the  frequency  which 
the  particular  magnitude  (i.e.,  Mi ,  mb,  Ms)  was  determined  for  some  reference 
event.  We  chose  the  event  at  801490516  as  our  reference  since  (l)  both  the 
source  spectrum  (from  the  close-in  body-wave  seismograms)  and  the  magni¬ 
tudes  are  well  determined,  (2)  it  is  one  of  the  smaller  events  and  consequently 
the  effects  of  fault  finitness  are  minimized,  (3)  the  frequencies  at  which  mb  and 
Ms  were  determined  are  lower  than  the  corner  frequency,  and  (4)  the  event  is 
shallow.  The  relationships  which  we  obtained  are: 

\ogQ  =  Mi  -6.25  , 

logG  =  mb  -  5.95  ,  (5) 

log  Q  =  Ms  -  5.00  . 

These  values  compare  well  with  those  found  by  recasting  equations  (3)  and  (4)  in 
terms  of  the  normalized  amplitude  G: 


logG0  =MS  -  5.01  , 


logfl0  =  1.2 Mi  -  6.93  . 

Figure  7  shows  the  result  of  applying  equation  (5)  to  the  magnitude  esti¬ 
mates  of  four  representative  Mammoth  Lakes  events,  which  sample  the  range  of 
depth  and  focal  mechanism  observed.  In  each  plot  the  source  spectrum  deter¬ 
mined  from  analysis  of  the  close-in  body-wave  data  is  shown  by  the  intersecting 
heavy  lines.  Error  bars  at  the  intersection  indicate  the  scatter  observed  in  ll0 
and  fc .  The  other  points  with  error  bars  arc  the  spectral  amplitude  values  com¬ 
puted  from  the  magnitude  estimates.  In  each  case  the  points  are  plotted  at  the 
mean  frequency  at  which  the  individual  magnitude  estimates  were  made  ( Ms  to 
the  extreme  left.  Mi  to  the  extreme  right,  and  mb  between  them).  The  error 
bars  for  these  points  also  indicate  the  observed  scatter.  Figure  7a  is  for  the 
reference  event. 
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The  801481901  event  (Fig.  7b)  is  similar  in  magnitude  and  depth  to  the 
reference  event.  Both  events  have  a  similar  corner  frequency.  mb  fall  below  the 
corner  frequency  and  on  the  flat  portion  of  the  spectrum;  Mi  falls  above  fc  and 
near  the  high  frequency  asymptote.  The  Ms  value  falls  below  the  long-period 
level,  but  this  is  not  significant  considering  the  error  bars  on  Ms  and  no. 

The  event  801471857  (Fig.  7c)  is  almost  one  magnitude  unit  larger  than  the 
reference  event  801490516,  but  is  of  similar  depth.  The  Ms  value  for  this  event 
falls  close  to  the  long-period  level  while  both  mb  and  Mi  are  near  the  high  fre¬ 
quency  asympote.  These  latter  two  measurements  were  made  at  frequencies 
within  the  error  bars  fc. 

The  event  801471224  (Fig.  7d)  occurring  at  12.5  km,  is  the  deepest  event 
studied,  as  well  as  one  of  the  deepest  events  of  the  Mammoth  Lakes  sequence. 
For  this  event,  Mi  falls  near  the  high  frequency  asymptote  above  the  corner  fre¬ 
quency,  while  mb  is  at  a  frequency  below  the  corner  and  approximately  a  factor 
of  two  beneath  the  long  period  level.  The  Ms  error  bars  intersect  the  long  period 
level  although  the  average  value  is  slightly  below,  possibly  due  to  the  increased 
source  depth. 

The  corner  frequencies  of  the  801471224  and  801471857  earthquakes  are 
about  a  factor  of  2  lower  than  801481901  and  q01490516,  whereas  the 
corresponding  low  frequency  amplitudes  are  approximately  10  times  higher. 
Thus  the  corner  frequencies  and  moment  are  in  agreement  with  the  u~z  model 
for  events  of  approximately  the  same  stress  drop,  but  with  source  dimensions 
different  by  a  factor  of  2  (Aki,  1967;  Brune,  1970,  1971).  The  above  comparison 
confirms  that  the  spectral  shape  inferred  from  distant  measurements  of  magni¬ 
tude  (Mi,  mb,  and  Ms),  is  congruent  with  the  spectral  shape  observed  directly  in 
the  near  field. 

7.0  Comparison  of  the  Mammoth  Lakes  results  with  published  results  for  other 
regions  of  western  North  America 

Stress  drop  may  vary  greatly  according  to  the  local  tectonic  conditions, 
and  consequently  a  single  ampliute  parameter  such  as  magnitude  may  be  a  poor 
representation  of  source  conditions.  The  relationship  between  the  seismic  spec¬ 
tra  and  the  parameters  Mi.  mb  and  Ms  derived  for  the  Mammoth  Lakes  earth¬ 
quakes  does  not  necessarily  apply  to  other  regions. 

Figure  8  compares  the  Mammoth  Lakes  data  for  moment  estimates  from 
near  source  data  vs.  Ms,  with  similar  data  for  two  published  studies.  The  line 
log Ma  -  M3  +  19.4  determined  for  the  Mammoth  Lakes  data  form  the  upper 
bound  to  the  San  Fernando  aftershock  observations  of  Tucker  and  Brune  (1977). 
The  two  events  from  the  Victoria,  Baja  California,  Mexico  swarm  reported  by 
Munguia  and  Brune  (1984)  lie  above  the  line.  Only  events  for  which  Ms  was  avail¬ 
able  could  be  included.  Thus  of  the  43  events  published  by  Tucker  and  Brune 
(1977)  for  the  San  Fernando  aftershock  sequence,  Figure  8  only  includes  12.  The 
San  Fernando  events  plot  in  a  cloud  to  the  right  of  the  line  defined  by  the  Mam¬ 
moth  Lakes  events,  i.e.,  they  have  larger  surface  wave  magnitudes  for  a  given 
local  body-wave  moment  than  do  events  of  the  Mammoth  Lakes  aftershocks. 
Tucker  and  Brune  found  that  there  were  other  events  which  had  body-wave 
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moments  of  similar  magnitude  to  the  events  plotted  in  Figure  8  but  for  which 
surface  waves  at  regional  distances  were  too  small  to  determine  magnitude. 
These  would  presumably  be  more  consistent  with  the  trend  of  the  Mammoth 
Lakes  data.  In  contrast  to  the  San  Fernando  events,  the  Victoria  earthquakes 
fall  above  the  line  determined  for  the  Mammoth  Lakes  events,  i.e.,  they  have 
smaller  Ms  than  would  be  predicted  from  their  body-wave  spectra.  Munguia  and 
Brune  (1984)  have  only  determined  Ms  for  two  of  the  43  events  they  have  stu¬ 
died.  These  events  have  stress-drops  of  930  and  570  bars. 

The  flat  spectra  observed  for  the  Mammoth  Lakes  earthquakes  from  surface 
wave  frequencies  up  to  the  spectral  corner  determined  from  the  body  wave  data 
indicates  that  the  majority  of  the  Mammoth  Lakes  events  studied  are  simple 
events.  (1985).  Tucker  and  Brune  (1977)  have  suggested  that  the  larger  San 
Fernando  earthquakes  (those  plotted  in  Figure  8)  are  partial  stress-drop  events, 
indicated  by  the  fact  that  their  surface  wave  magnitudes  are  larger  than  would 
be  predicted  based  on  the  body-wave  determination  of  M0.  This  results  from  a 
rise  in  the  spectrum  at  frequencies  below  the  body-wave  corner  frequency.  The 
larger  events  of  the  Victoria  swarm  are  high  stress  drop  events  (Munguia  and 
Brune,  1984)  and  have  smaller  surface  wave  magnitudes  than  would  be  predicted 
from  their  body  wave  spectra. 

Hartzell  and  Brune  (1979)  have  proposed  an  alternative  model  for  explain¬ 
ing  the  difference  observed  in  low  and  high  frequency  moments.  They  suggest 
that  for  these  events,  faulting  initiates  with  a  rupture  over  a  relatively  small 
area  with  a  rapid  dislocation  rate.  This  first  phase  of  rupture  accounts  for  the 
smaller  body-wave  («  1  Hz)  estimate  of  moment.  If  the  rupture  then  continues 
to  expand,  but  at  a  much  slower  dislocation  rate,  most  of  the  energy  will  be  at 
surface  wave  frequencies  («  0.1  Hz)  and  thus  account  for  the  larger  moment 
based  on  low  frequency  waves.  In  this  case,  Figure  8  should  be  looked  at  in  a 
different  light.  Rather  than  Ms  being  over  estimated  for  a  given  body-wave 
moment,  the  body-wave  moment  is  underestimated  for  these  event.  Das  and 
Kostrov  (1985)  have  proposed  a  similar  model  to  that  of  Hartzell  and  Brune 
(1979),  in  which  such  earthquakes  consist  of  the  fracture  of  a  single  asperity  fol¬ 
lowed  by  slip  on  a  finite  fault.  They  find  that  for  such  a  model,  the  seismic 
moment  is  ( R/r )2  larger  and  the  corner  frequency  ( r/R )  lower  (where  r  is  the 
radius  of  the  asperity,  and  R  is  the  radius  of  the  circular  finite  fault)  than  for  a 
simple  dislocation  earthquake  of  the  same  magnitude.  Since  the  asperity  is  a 
one-sided  pulse,  the  envelope  of  its  spectrum  must  have  the  same  general  form 
as  the  spectrum  of  the  simple  dislocation  model;  that  is  its  corner  frequency  is 
determined  by  the  pulse  duration,  its  moment  is  determined  by  the  area  under 
the  pulse,  and  it  has  a  high  frequency  slope  of  -2.  The  major  difference  between 
the  pulse  due  to  an  asperity  fracture  on  a  finite  fault  and  a  simple  crack  are  the 
anomalously  long  pulse  duration  and  the  anomalously  large  seismic  moment  of 
the  former  compared  to  the  latter. 

We  might  now  speculate  as  to  why  differences  such  as  those  observed 
between  the  Mammoth  Lakes  and  San  Fernando  sequences  arise.  Lide  and  Ryall 
(1935)  have  shown  that  in  the  weeks  following  the  main  shocks  of  the  Mammoth 
Lakes  sequence  (the  time  period  for  the  events  of  this  study),  the  majority  of 
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the  aftershocks  were  located  in  the  granitic  Sierrian  block  and  were  distributed 
in  a  triangular  shaped  zone  above  a  steeply  dipping  plane  on  which  the  Mi>  6 
main  shocks  possibly  occurred.  In  comparison,  a  north-south  cross-section 
(perpendicular  to  the  surface  rupture)  of  the  San  Fernando  events  analysed  by 
Tucker  and  Brune  (1977)  lie  in  a  granitic  block  above  a  zone  dipping  at  33°  (Allen 
et  at.  (1973).  In  both  the  case  of  Mammoth  Lakes  and  San  Fernando,  the  events 
which  exhibit  anamolously  large  surface  wave  excitation  compared  to  body  wave 
moments  are  concentrated  near  the  bottom  along  the  dipping  zones  which  may 
be  considered  as  a  through-going  fault  plane,  while  events  which  would  be  con¬ 
sidered  as  full  stress  drop  or  simple  dislocation  events  are  concertrated  in  the 
shallow  crust.  These  data  suggest  that  partial  stress  drop  events  are  more  com¬ 
mon  on  developed  fault  systems,  whereas  full  stress  drop  events  are  more  com¬ 
mon  in  the  surrounding  country  as  it  readjust  to  the  stress  field  generated  by 
the  movement  on  the  developed  fault  plane.  Alternatively,  the  granitic  blocks 
above  both  the  Mammoth  Lakes  and  San  Fernando  main  fault  plane  may  be  so 
highly  shattered  that  stress  can  only  be  consentrated  on  planes  of  limited 
extent.  These  "intrablock"  earthqukes  occurring  in  the  granitic  block  above  the 
through-going  fault  at  Mammoth  Lakes  and  San  Fernando  may  be  higher  or  full 
stress  drop  events  whereas  the  "interblock”  events  occurring  on  the  through 
going  fault  are  partial  stress  drop  events  or  fracture  of  asperities  on  finite 
faults. 

8.0  Conclusions 

In  comparing  short  period  versus  long  period  excitation  of  events  in  the  Cal¬ 
ifornia  region,  Wyss  and  Brune  (1971)  found  that  events  occurring  to  the  north  of 
Bishop,  California  indicated  higher  apparent  stress  than  events  along  the  San 
Andreas  fault.  In  this  study  we  have  looked  in  detail  at  the  relative  excitation  of 
low  (0.05  -  0.125  Hz)  and  high  (0.5  -  10.0  Hz)  frequency  seismic  energy  in  this 
region  from  aftershocks  of  the  1980  Mammoth  Lakes  earthquake  sequence. 

We  have  compared  the  long  period  spectral  level  Q0  computed  from  body 
waves  recorded  on  strong-motion  and  digital  instruments  at  near  epicentral  dis¬ 
tance  (~  1  Hz)  with  Q0  computed  from  regionally  recorded  surface  waves  (~  0.1 
Hz)  and  find  these  two  values  to  agree  within  the  observational  error  for  almost 
all  events  of  the  Mammoth  Lake  sequence. 

We  have  compared  the  body  wave  magnitude,  mb,  and  surface  wave  magni¬ 
tude,  i\fSl  and  find  that  Mammoth  Lakes  events  discriminate  from  nuclear  explo¬ 
sions  occurring  at  the  Nevada  Test  Site  approximately  150  km  to  the  southeast. 

Body-wave  moment  (Af0)  has  been  related  to  local  magnitude  {\f ^ )  by  the 
relation 

log  Ma  =  (1.2  ±  0.05)  +  (17.49  -0.19)  (l  <  AfL  <6) 

and  surface  wave  moment  has  been  related  to  surface  wave  magnitude  (Afs )  by 
the  relation 


log  M0  =MS  +  (19.40  ±  0. 00)  . 


Mi,  mb,  and  Ms  for  four  representative  events  were  compared  with  the 
spectral  amplitude  at  their  appropriate  frequencies.  The  three  shallow  events 
(801471857,  801481901.  and  801490516)  and  the  one  deep  event  (801471224)  are 
consistent  with  the  cj-2  model.  The  stress  drops  of  these  events  are  nearly  simi¬ 
lar  but  the  source  dimension  of  the  801471857  and  801471224  events  are  nearly 
twice  that  of  the  other  two  events,  which  explains  their  larger  moment  and 
greater  surface  wave  excitation. 
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Figure  Captions 

Figure  1.  Location  maps  showing  the  events  studied  and  the  stations  from  which 
seismograms  were  analysed.  In  the  lower  right-hand  map,  the  solid  circles 
denote  locations  of  events,  triangles  denote  location  of  strong-motion 
accelerographs,  and  digital  seismographs.  The  upper  left-hand  map  shows 
the  locations  of  the  long-period  seismographs  used  for  the  surface  wave 
analysisd  (denoted  by  squares),  and  Canadian  network  short-period  sta¬ 
tions  used  to  determine  body-wave  magnitude  (denoted  by  stars). 


Figure  2.  Representative  analysis  of  a  short-period  digital  seismogram.  The 
bracket  over  the  seismogram  indicates  the  S-wave  window;  the  heavy  solid 
lines  indicate  the  long-period  and  high-frequency  asymptotes  fit  to  the 
spectrum. 


Figure  3.  Represe'  ative  fit  of  computed  surface  wave  seismogram  to  the 
observed  seismogram  for  the  WWSSN  station  at  Dugway.  Utah  (DUG). 


Figure  4.  Comparison  of  the  body-wave  and  surface-wave  long  period  spectral 
levels  for  the  events  studied.  The  line  is  for  reference  and  has  unit  slope 
and  zero  intercept. 


Figure  5.  Comparison  of  the  body-wave  and  surface-wave  magnitude.  The  solid 
dots  (earthquakes)  and  pluses  (explosions)  are  taken  from  Marshall  and 
Basham  (1972);  the  solid  squares  are  the  Mammoth  Lakes  events  from  this 
study.  All  body-wave  magnitudes  are  from  measurements  from  Canadian 
Network  seismographs.  Surface  wave  magnitudes  taken  from  Marshall  and 
Basham  are  from  Canadian  Network  long-period  seismographs;  surface 
wave  magnitudes  for  the  Mammoth  Lakes  earthquakes  are  from  long- 
period  WWSSN  stations  in  the  western  United  States. 


Figure  6.  Log  M0  vs.  A fi  for  earthquakes  in  the  1930  Mammoth  Lakes  earthquake 
|  sequence.  One  standard  deviation  of  the  mean  error  bars  are  given  for 

events  with  multiple  recordings.  The  line  is  from  equation  (4),  and  the 
stars  are  the  theoretical  values  from  Hanks  and  Boore  (1934).  (After 
Chavez  and  Priestley,  1985.) 

I 

Figure  7.  Comparison  of  the  source  spectra  and  the  amplitudes  determined  from 
Mi ,  m6 .  Ms  at  their  appropriate  frequencies.  The  error  bars  assigned  to 
each  spectral  value  represent  the  estimated  standard  error  of  the  meas¬ 
urement:  (a)  801490516  (event  #40),  (b)  801481901  (event  #34),  (c) 

I 
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801471857  (event  #22).  (d)  801471224  (event  #18). 


Figure  8.  Comparison  of  logAf0  vs.  Ms  for  Mammoth  Lakes  aftershock  data,  and 
log M0  vs.  Ms  for  several  other  western  North  American  earthquake 
sequences.  The  line  plotted  is  from  equation  (3). 
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ABSTRACT 


We  have  determined  an  average,  frequency  dependent,  apparent  Lg 
Q  (f )  function  for  Great  Basin  paths,  by  measuring  the  spatial 
decay  of  spectral  amplitudes.  We  find  that 

{?(/)  =  212/ 0  83  (0.2  Hz  ^  /  g  4.0  Hz). 

This  result  is  in  excellent  agreement  with  Singh  and  Herrmann’s 
(1983)  Q  model  for  this  same  area,  based  on  observations  of  coda 
decay. 


Introduction 

Lg  waves  are  short-period  surface  waves  which  propagate  in  the  continental 
crust  with  typical  group  speeds  of  3.0  to  3.5  km/sec  (Press  and  Ewing,  1952; 
Bath,  1954).  Knopoff  ef  al.  (1973)  have  identified  the  transverse  component  of  Lg 
with  higher  mode  Love  waves  while  Panza  and  Calcagnille  (1975)  have  shown  that 
the  vertical-radial  component  can  be  identified  with  higher  mode  Rayleigh 
waves. 

Currently,  most  estimates  of  seismic  moment  for  Great  Basin  earthquakes 
have  been  made  using  body  waves  recorded  by  local  networks  or  fundamental 
mode  surface  waves  recorded  at  regional  distances.  Numerous  moderate  size 
earthquakes  (3<M £,<4)  occur  in  the  Great  Basin  outside  of  any  local  network  and 
are  too  small  to  be  seen  at  teleseismic  distances.  They  are,  however,  recorded 
by  stations  at  regional  distances  where  Lg  is  almost  always  a  clearly  recognized 
arrival.  Street  et  al.  (1975)  and  Herrmann  and  Kijko  (1983)  determined  empiri¬ 
cal  formulas  which  relate  the  long  period  asymptote  of  the  Lg  spectrum  to  the 
seismic  moment.  They  used  data  from  the  central  and  eastern  United  States 
where  the  low  attenuation  allowed  them  to  neglect  any  Q  effects.  There  is  sub¬ 
stantially  greater  attenuation  in  the  Great  Basin  however,  so  a  frequency  depen¬ 
dent  Q  model  must  be  obtained  in  order  to  correct  regional  Lg  recordings  prior 
to  determining  the  seismic  moment. 

Nuttli  (1973,  1980,  1981,  1984)  has  used  Lg  to  estimate  mb  for  earthquakes 
and  explosions,  estimate  apparent  crustal  attenuation,  discriminate  between 
small  earthquakes  and  explosions,  and  estimate  the  yield  of  nuclear  explosions. 


In  the  latter  case,  he  found  that  if  Lg  amplitudes  are  to  be  used  to  estimate 
explosion  yields  to  an  accuracy  of  30%,  then  average  values  of  spatial  attenua¬ 
tion  over  large  areas  (such  as  the  western  U.S.)  are  inadequate.  Instead, 
smaller  scale  regional  estimates  must  be  used  in  order  to  account  for  lateral 
variations  in  the  attenuation  structure. 

While  Lg  has  been  shown  to  be  a  very  useful  phase  for  determining  a  variety 
of  seismological  quantities,  its  utility  is  diminished  in  the  presence  of  appreci¬ 
able  attenuation.  In  an  effort  to  quantify  the  attenuation  affecting  Lg  propaga¬ 
tion  for  Great  Basin  paths,  we  examine  recordings  of  Lg  along  several  paths  and 
measure  the  spatial  decay  of  spectral  amplitudes.  The  data  provide  an  estimate 
of  the  average  Lg  Q  in  the  Great  Basin. 

Data  and  Analysis 

The  data  set  for  this  study  consists  of  regional,  digital  seismograms  for  12 
earthquakes  in  and  around  the  Great  Basin  recorded  at  four  broad-band  seismo¬ 
graph  stations  operated  by  the  Lawrence  Livermore  National  Laboratory.  Figure 
1  shows  the  locations  for  the  earthquakes  and  seismic  stations  and  the  travel 
paths  between  them.  We  selected  earthquakes  which  were  recorded  by  at  least 
three  of  the  four  stations  and  which  covered  a  large  range  in  epicentral  distance 
(le.  no  two  stations  were  at  a  similar  distance).  As  indicated  by  Figure  1,  the 
travel  paths  used  provide  a  fair  sampling  of  the  Great  Basin  crust. 

The  digital  seismographs  are  velocity  flat  from  0.7  to  10  Hz  and  operate  at  a 
sample  rate  of  40  samples/second.  However,  due  to  memory  limitations  in  the 
computer  used  to  process  the  data,  we  were  forced  to  decimate  the  data  down 
to  a  sample  rate  of  10  samples/second;  consequently  we  have  no  resolution 
above  5  Hz.  Figure  2  is  the  magnification  curve  for  the  seismograph  system 
used. 

A  typical  set  of  recordings  is  plotted  in  Figure  3,  which  shows  the  original 
(decimated)  data  and  the  same  data  after  high-pass  filtering  by  convolving  the 
mstrument  corrected  signal  with  the  response  of  a  WWSSN  short-period 
seismometer.  The  Lg  phase  is  indicated  on  the  seismograms,  and  can  be  seen  to 
be  a  prominant  feature  on  the  short-period  record. 

Wc  have  measured  Lg  group  speeds  for  250  regional  seismograms  with  Great 
Basin  paths,  including  those  used  in  this  study.  Figure  4  is  a  histogram  of  the 
group  speeds  obtained,  and  shows  that  Lg  in  the  Great  Basin  typically  propagates 
with  a  speeds  between  3.3  and  3.6  km/scc.  The  most  common  speed  measured 
is  3.5  km/scc. 

To  determine  the  Lg  attenuation,  we  computed  instrument  corrected,  dis¬ 
placement  amplitude  spectra  of  the  vertical  component  Lg  phase.  We  inverted 
these  spectra  for  the  Q(f)  function  in  the  following  manner.  The  observed  Lg 
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spectrum  is  modeled  as 

a(j.r)  =  e  °(J/T  (0 

where  A(f  ,R)  is  the  spectral  amplitude  observed  at  a  distance  R ,  /  is  fre¬ 
quency,  t  is  the  travel  time,  S(f )  is  the  source  term,  and  Q  is  the  quality  factor. 
This  simple  model  neglects  radiation  pattern  and  scattering  effects,  therefore  it 
provides  a  measure  of  the  apparent  rather  than  the  intrinsic  Q .  Since  Lg  has 
been  succesfully  modeled  as  surface  wave  (Knopoff  et  al.,  1973;  Panza  and  Cal- 
cagnille,  1975)  we  have  assumed  that  the  frequency  domain  geometrical  spread¬ 
ing  scales  with  the  square  root  of  distance.  Taking  the  logarithim  (base  10)  of 
both  sides  of  (1)  yields 

log  i0A  (J.R)  +  log  10  R  05  =  log  10  S  (/  )  -  j  1 1  (2) 

which  is  the  equation  for  a  straight  line,  where  the  source  term  is  the  intercept 
and  the  Q  term  controls  the  slope.  Fixing  f,  we  know  A(f  ,R),  R.  and  t  for  each 
of  the  stations  and  we  solve  for  5  and  Q  using  least-squares,  dy  looping  over  all 
the  frequencies  we  obtain  the  source  and  Q  spectra. 

We  initially  determined  the  travel  time  for  use  in  equation  2  by  taking  the 
travel  time  to  the  center  of  the  window  which  was  used  to  compute  the  spec¬ 
trum.  However,  we  found  that  we  obtained  larger  correlation  coefficients  overall 
if  we  computed  the  travel  time  using  a  speed  of  3.5  km/sec.  The  results 
presented  here  were  obtained  using  the  latter  method,  but  there  is  no 
significant  change  in  the  results  if  the  former  method  is  used.  This  suggests 
that  errors  in  travel  time  due  to  errors  in  the  earthquake  locations  or  origin 
times  do  not  strongly  influence  the  conclusions  presented  here. 

Since  for  any  given  earthquake  we  are  fitting  only  three  or  four  points  at 
each  frequency,  the  results  obtained  can  be  expected  to  have  significant  random 
error.  This  will  be  particularly  true  if  the  data  contain  any  variations  in  ampli¬ 
tude  due  to  radiation  pattern.  We  have  attempted  to  compensate  for  these  pos¬ 
sible  errors  by  rejecting  any  results  for  which  the  linear  correlation  coefficient 
was  below  0.8.  For  several  events,  this  eliminated  the  data  at  lower  frequencies 
which  were  affected  by  the  6-second  microsetsms.  In  general,  the  data  at  higher 
frequencies  gave  more  stable  results. 

Figure  5a  shows  all  the  Q  (/ )  functions  obtained  with  correlation 
coeffiecients  of  0.3  or  greater.  The  average  Q  (J  )  function  is  shown  in  Figure  5b. 
A  linear  least-squares  fit  to  the  data  iri  Figure  5b  yields 

Q(f  )  -  212/°  S3  (0.2  Hz  §  /  &  4.0  Hz).  (3) 
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Discussion  and  Conclusions 

It  is  interesting  to  compare  the  result  obtained  above  with  other  published 
Q  functions  for  the  Great  Basin.  Cheng  and  Mitchell  (1991)  used  long  period  Lg 
waves  to  estimate  a  constant  Q  value  of  85.  Their  data  indicate  a  predominate 
period  of  around  5  seconds  and  so  we  can  consider  their  Q  to  be  the  average 
crugtal  value  at  D.2  Hz.  This  is  in  agrees  well  with  the  value  of  90,  predicted  by 
equation  (3). 

Our  Q  {f )  function  is  also  in  excellent  agreement  with  that  found  by  Singh 
and  Herrmann  (1983),  who  determined  regional  frequency  dependent  Q  func¬ 
tions  for  the  United  States  using  measurments  of  coda  decay.  Their  result  for 
the  Great  Basin  in  the  vicinity  of  our  profile  is  f?  (/ )  as  200/  (°-45±0-05).  Singh  and 
Herrmann  were  able  to  relate  their  observations  of  coda  decay  to  the  Lg  spatial 
attenuation  coefficient  and  concluded  that  the  same  attenuation  mechanism  was 
responsible  for  both  Lg  and  coda  decay.  Our  results  lend  support  to  this  conclu¬ 
sion. 
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Figure  captions 

Figure  1.  Map  showing  locations  of  the  earthquakes  (stars)  and  recording  sites 
(squares)  used  in  this  study  and  the  travel  paths  between  them. 

Figure  2.  Magnification  curve  for  the  seismograph  system  used. 

Figure  3.  Examples  of  typical  seismograms  recorded  at  each  of  the  stations. 
Each  record  has  been  individually  normalized  to  the  same  peak  plot  ampli¬ 
tude. 

Figure  4.  Histogram  of  group  speed  measurments  for  250  recordings  of  Lg  with 
Great  Basin  paths. 

Figure  5.  a)  Summary  of  Q  (/ )  functions  obtained  for  each  earthquake.  Values 
whose  associated  linear  correlation  coefficient  was  below  0.8  have  not  been 
plotted,  b)  Average  Q  (/ )  function  for  the  data  in  fa).  One  standard  devia¬ 
tion  error  bars  are  included,  as  is  the  least-squares  power  curve  fit  to  the 
data  (e  juation  3). 
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Introduction 

On  25  and  27  May  1980  three  ML  6+  earthquakes  occurred  in  the  western  Great 
Basin  near  Mammoth  Lakes,  California,  about  200  km  NW  of  the  Nevada  Test  Site.  As 
part  of  our  investigation  of  these  earthquakes  we  quested  seismograms  from  several 
seismic  stations  in  the  Soviet  Union,  including  a  station  at  Semipalatinsk,  approximately 
100  km  NE  of  the  eastern  Kazakh  test  site.  This  note  summarizes  mk  determinations 
for  five  events  on  the  Soviet  records  —  the  three  large  Mammoth  Lakes  events,  a  smaller 
(mk  5.3)  Mammoth  Lakes  shock,  plus  an  earthquake  at  Tonga  (mt  6.0)  that  occurred 
during  the  time  frame  of  the  recordings.  The  records  were  searched  for  eight  more 
events  on  the  ISC  list,  but  none  of  these  were  recorded  at  the  Semipalatinsk  station  and 
most  were  not  recorded  by  the  other  stations.  Another  strong  (m*  5.4)  Mammoth  Lakes 
shock  at  16:49  GMT  on  25  May  and  an  earthquake  in  the  Kurile  Islands  (25  May,  23:22 
GMT,  mt  4.8)  were  recorded  but  not  analyzed.  Table  1  summarizes  information  on  the 
events  used  in  this  study. 


Table  1.  Llat  of  Event* 

Date 

Time,  GMT 

Latitude 

Longitude 

Depth 

m4 

M. 

800525 

16:33:44.7 

37.596 '  N 

118.830  *  W 

6.5 

6.1 

6.1 

800525 

19:44:51.1 

37.547  *N 

118.826*  W 

5.0 

5.6 

6.0 

800525 

20:35:48.5 

37.616  ’  N 

118.847*  W 

6.1 

5.3 

5.7 

800527 

14:50:57.1 

37.472  •  N 

118.807*  W 

10.8 

5.7 

6.0 

800527 

13:01:37.9 

18.610  *S 

174.700  *  E 

55 

6.0 

4.8 

Data 

The  data  included  recordings  from  eight  seismic  stations  on  a  northeast-trending 
profile  from  Kzyl-Agach  in  the  southern  part  of  the  Kazakh  Fold  System  to  Seymchan  in 
the  Northeast  Siberian  Fold  System.  Table  2  lists  the  station  coordinates  and  Figure  1 
is  a  polar  projection  of  Asia  showing  the  station  locations.  Table  3  gives  epicentral  dis¬ 
tances  and  azimuths  for  the  five  events. 


Table  3.  Station  Location* 

Code 

Name 

Coordinates 

Region 

BOD 

Bodaibo 

57  *  51’N,114  *  ll’E 

Baikal 

ELT 

El’tsorka 

53  *  15’N,86  *  I6’E 

Altai-Sayan 

IRK 

Irkutsk 

52  *  16’N,104 ' 19’E 

Baikal 

KZL 

Ksyl-Agach 

45  *  25’N,78  *  45’E 

Kasakhstan 

SEI 

Seymchan 

62  ’  53’N,152  *  26’E 

Northeast 

SEM 

Semipalatinsk 

50  *  24’N,80  ’  15’E 

Kazakhstan 

UST 

Ust’-Kan 

50  *  57’N,84 ' 45’E 

Altai-Sayan 

YAK 

Yakutsk 

62  *  01’N,129  *  43’N 

Yakutiya 

*  —  Work  done  at  Center  for  Seismic  Studies  and  partly  supported  by  Science  Applica¬ 
tions  International  Corporation. 


Figure  2  is  a  polar  plot  of  the  world  centered  on  a  point  (81.5  *  N,  162.0*  W)  about 
halfway  between  the  Mammoth  Lakes  earthquakes  and  station  SEM,  showing  the  epi¬ 
center  and  recording  stations.  Figure  3  shows  the  position  of  the  Soviet  stations  on  a 
lower-hemisphere,  equal-area  projection  of  the  focal  sphere  for  the  three  largest  Mam¬ 
moth  Lakes  earthquakes,  together  with  projections  of  the  fault  and  auxiliary  planes  for 
these  events  from  a  study  of  long-period  P-  and  surface-waves  (Given  et  al .,  1982). 
Based  on  point-source  theory  (Keilis-Borok,  1950)  and  the  position  of  the  stations  on  the 
focal  sphere,  P-waves  to  the  Soviet  stations  should  have  about  80%  of  the  maximum 
radiated  amplitude. 

Figure  4  is  a  polar  projection  centered  on  a  point  (23.0  *  N,  147.0  *  E)  about  halfway 
between  station  SEM  and  the  Tonga  earthquake.  The  Kurile  event  was  in  a  poor  dis¬ 
tance  range  (16-47*)  for  this  study,  and  and  had  large  scatter  in  mk  values;  it  will  not 
be  considered  in  the  detailed  discussions  that  follow.  For  the  Tonga  earthquake  a  plot 
similar  to  Figure  4  indicates  that  the  Soviet  stations  are  close  to  the  null  axis  on  the 
fault-plane  solution,  although  the  latter  is  not  well-constrained. 

Table  4  lists  instrument  parameters  for  the  86  recordings  that  were  supplied. 
Values  of  the  maximum  magnification  Vm  and  the  period  range  Tm  (corresponding  to 
V  =  0.9  Vm )  were  written  on  the  records.  For  station  KZL  the  period  range  Tm  was 
not  supplied  and  the  range  given  in  the  table  was  taken  from  the  Soviet  publication 
Earthquakes  in  the  USSR  in  1980  (Akademiya  Nauk  SSSR.  983).  For  station  BOD  the 
values  of  Vm  marked  on  the  recordings  —  "2900"  for  the  ho  component  and  "2700"  for 
the  Z  and  EW  —  appeared  to  be  an  order  of  magnitude  too  low,  and  the  magnification 
(52,000)  for  that  station  was  also  taken  from  Earthquakes  in  the  USSR  in  1980.  Use  of 
the  larger  magnification  for  BOD  is  supported  by  Shishkevish  (1974)  who  lists  Vm  as 
45,000-49,000  for  this  station  in  1970. 

For  this  study  only  vertical-component  recordings  from  SKM-3  seismographs  were 
used.  The  SKM-3  system  consists  of  a  seismometer  and  galvanometer  with  instrument 
constants  designed  to  produce  magnification  of  ground  motion  in  the  range  30,000- 
100,000  over  the  period  range  0.3-1. 5  seconds.  According  to  Aranovich  et  al.  (1974)  a 
complete  system  consists  of  two  SGKM-3  (horizontal)  and  one  SVKM-3  (vertical) 
seismometers,  three  GK-VIIM  galvanometers  and  a  PS-3M  drum  recorder.  The  latter 
registers  the  light  beams  from  the  three  galvanometers  on  a  29-cm  wide  by  90-cm  long 
photographic  recording.  Time  marks  from  a  chronometer  are  printed  once  per  minute. 
Drum  speed  was  60  mm/minute  for  all  of  the  SKM-3  records  used  in  this  study  except 
KZL,  which  recorded  at  120  mm/minute.  Figure  5  shows  relative  amplitude  response  for 
the  eight  stations,  determined  from  equations  given  by  Aranovich  et  al.  Note  that  the 
Soviet  short-period  seismographs  peak  at  somewhat  lower  frequency  (ca.  1  Hz)  than  sys¬ 
tems  commonly  used  in  the  US,  and  that  for  some  stations  (t.e.,  SEM  and  IRK  in  this 
study)  the  system  gain  is  reduced  by  about  a  factor  of  ten  at  10  Hz. 
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Seismometer  and  galvanometer  constants  (period,  damping)  for  stations  BOD,  IRK 
and  SEM  were  taken  from  Shishkevish  (1974),  and  the  resulting  response  was  checked  to 
insure  that  the  range  of  Tm  was  the  same  as  that  indicated  on  the  recordings.  For  ELT 
the  constants  given  by  Shishkevish  were  modified  to  give  the  appropriate  upper  limit  of 
Tm ,  and  the  same  constants  were  used  for  UST.  For  YAK  Shishkevish  lists  instrument 
constants  only  for  the  horizontal-component  seismographs,  and  these  give  a  different 
range  of  Tm  than  that  specified  on  the  records.  For  KZL  and  SEI  Shishkevish  does  not 
list  instrument  constants.  As  a  result,  for  KZL,  SEI  and  YAK  various  combinations  of 
instrument  constants  were  tried,  based  on  tables  of  standard  setups  given  by  Aranovich 
et  al. ,  until  response  curves  matched  the  range  of  Tm  marked  on  the  records. 

Except  for  the  SKM  recordings  at  station  KZL  and  the  long-period  recordings  at 
UST  all  of  the  records  were  12  hours  long.  The  KZL  record  has  twice  the  time  resolu¬ 
tion  of  the  other  short-period  stations  and  appears  to  be  changed  three  times  a  day;  the 
SKD  record  for  UST  for  25  May  was  24  hours  long.  Unfortunately  the  records  of  pri¬ 
mary  interest  to  this  analysis  —  from  station  SEM  —  were  changed  at  different  times 
than  those  at  other  stations  and  overlap  with  the  latter  for  only  two  six-hour  periods. 


Table  4.  Data  Received  for  25  and  27  May  1980 


BOD 

SKM-3 

ELT 

SKM-3 

SKM-3 

IRK 

SKM-3 

SKM-3 

SKM-3 

SKD 

KZL 

SKM-3 

SKM-3 

SKM-3 

SKM-3 

SEI 

SKM-3 

SKD 

SEM 

SKM-3 

SKM-3 

SK 

SK 

UST 

SKM-3 

SKM-3 

SKD 

YAK 

SKM-3 

SKM-3 

SKM-3 

SK-KPCh 

SK 

SK 

SK 

Comp 


NS.Z.EW  52,000 


NS.Z.EW 


0.2  -  1.2 


NS.Z.EW 

NS.Z.EW 


NS,EW 

Z 

NS.EW 

Z 


NS.Z.EW 

E-W 

NS.Z.EW 


Z 

NS 

EW 

NS 

NS 

Z 

EW 


1.1  - 

1.6 

1.1  - 

1.6 

1.1  - 

1.6 

0.2 

-  20 

0.08 

-  1.6 

0.08 

-  1.6 

0.08 

-  1.6 

0.08 

-  1.2 

0.2  - 

1.2 

0.2  ■ 

■  20 

0.84 

-  1.5 

0.8  - 

1.5 

0.41  - 

10.9 

0.5  ■ 

■  11 

0.2  - 

1.4 

0.2  - 

1.4 

0.2  ■ 

•  22 

Record  Quality.  Copies  of  the  records  were  on  35  mm  film  of  poor-to-good  qual¬ 
ity.  For  many  of  the  short-period  records  the  contrast  between  the  trace  and  the  back¬ 
ground  was  poor,  and  attempts  to  make  enlargements  of  the  waveforms  using  a  reader- 
printer  were  not  successful.  As  a  result  the  records  had  to  be  projected  on  a  viewing 
screen  and  traced  by  hand;  the  tracings  were  then  digitized  for  computer  analysis. 


'•  .'V  J 
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Quality  of  the  KZL  records  was  especially  poor,  and  the  P-waveform  for  only  one  event 
was  traced  for  that  station. 

Timing.  Two  or  three  time  corrections  were  marked  on  each  of  the  records,  usu¬ 
ally  corresponding  to  the  beginning  and  end  of  the  recording  period.  Clock  corrections 
changed  at  most  stations  by  less  than  a  second  per  day.  A  few  records  were  mislabeled 
as  to  sense  of  the  time  correction,  and  in  one  case  the  time  corrections  at  the  two  ends  of 
a  recording  period  were  interchanged  on  different  records  of  the  same  station.  The  larg¬ 
est  clock  drift  was  0.5  second/hour  for  the  SEM  station,  but  this  rate  was  constant. 
Tracings  of  the  P-waves  were  enlarged  such  that  the  average  time  scale  was  about  270 
mm/minute.  Station  KZL  had  a  drum  speed  twice  that  of  the  other  stations,  so  the 
enlarged  traces  were  viewed  at  a  scale  of  540  mm/min.  For  most  of  the  events  the 
beginning  of  the  P-wave  (PKP  for  stations  at  A  >  110  * )  was  identified  on  records  of  one 
of  the  better  stations  (SEI,  YAK,  BOD)  and  the  same  point  was  marked  and  timed  on 
traces  for  the  other  stations  by  overlaying  the  recordings. 

Traveltime  residuals  were  calculated  for  the  Mammoth  Lakes  earthquakes  using 
hypocentral  coordinates  and  origin  times  determined  by  the  University  of  Nevada  from 
local  network  recordings  and  traveltimes  calculated  from  the  Herrin  et  al.  (1968)  tables. 
The  Soviet  stations  had  average  delays  of  0.5-2. 2  seconds  for  the  Mammoth  Lakes 
events,  relative  to  the  Herrin  tables.  These  delays  are  not  considered  to  be  meaningful 
for  the  present  study,  since  the  larger  earthquakes  of  this  sequence  appear  to  have  been 
multiple  events,  initiated  by  small  ruptures  that  recorded  at  regional  stations  but  not  at 
teleseismic  distance  ranges  ( e.g .,  Given  et  al .,  1982).  On  the  other  hand,  traveltime  resi¬ 
duals  for  individual  stations  relative  to  the  average  delays  for  all  the  stations  (Table  5) 
are  fairly  consistent  for  the  Mammoth  Lakes  and  Tonga  earthquakes.  Note  that  the  two 
stations  north  of  Lake  Baikal  (BOD,  IRK)  have  early  arrivals  for  waves  travelling 
southwest  across  the  Siberian  platform  from  the  Mammoth  Lakes  events,  but  are  late  for 
paths  from  Tonga  that  cross  the  Baikal  graben.  Station  SEM  has  small  traveltime  resi¬ 
duals  that  average  to  almost  zero. 


Table  5.  Relative  P-wave  Residuals,  Seconds 

Sta 

251633 

251944 

252035 

271450 

271301 

Average 

BOD 

-1.0 

-0.3 

-1.2 

-0.8 

+0.5 

-0.6+0.7 

EXT 

sa 

-1.5 

-i.3+0.4 

IRK 

-0.9 

-0.3 

+0.2 

-0.5+0.5 

KZL 

0.0 

^1 

^  ^ 

SEI 

+0.8 

+0.5 

+0.6 

+  1.4 

+0.8 

+0.8+0.3 

SEM 

+0.5 

-0.5 

-0.2 

+0.5 

+0.5 

+0.2+0.5 

UST 

+0.5 

-0.5 

-0.2 

+0.5 

-0.7 

-0.1+0.6 

YAK 

mm 

+0.6 

-0.1 

+0.6 

-0.0 

+0.2+0.4 

Magnitude  Determination 

Figure  6  shows  an  example  of  the  digitized  P-waves  for  the  first  large  Mammoth 
Lakes  event,  at  16h  33m  on  25  May  1980.  All  traces  tire  from  SVK-M3  instruments,  and 
are  normalized;  up  on  the  traces  corresponds  to  ground  motion  up.  With  the  exception 
of  KZL  most  of  the  recordings  had  good  signal-to-noise  ratios.  Stations  YAK  and  BOD 
appear  to  have  higher  frequency  content  than  some  of  the  other  stations,  including  SEM, 
but  this  could  be  a  result  of  differences  in  instrument  response.  P-wave  spectra  for  this 
event,  corrected  for  instrument  response,  have  peaks  at  about  0.5  and  0.7  Hz;  for  station 
BOD  the  peaks  are  about  equal  in  amplitude,  but  for  the  other  stations  the  amplitude  at 
0.7  Hz  is  less  than  half  that  at  0.5  Hz.  In  general  the  character  of  the  P-wave  at  the 


various  stations  is  similar,  consisting  of  a  small  first  arrival  followed  6  seconds  later  by  a 
larger  phase.  The  time  interval  between  these  phases  is  too  consistent  for  the  second 
arrival  to  be  PcP,  since  for  the  distance  range  of  the  Soviet  stations  t(PcP-P)  should 
vary  from  60  seconds  to  zero.  The  second  phase  could  be  another  earthquake,  but  this  is 
unlikely  because  all  of  the  Mammoth  Lakes  events  have  a  large  second  arrival  at  the 
Soviet  stations.  A  more  likely  explanation  is  that  the  phase  includes  the  reflected  waves 
pP  and  and  sP.  Based  on  the  known  crustal  structure  in  the  Mammoth  Lakes  area,  a 
pP-P  time  of  6  seconds  would  be  consistent  with  focal  depth  of  about  16  km,  and  for 
that  depth  the  sP-P  time  would  be  about  8  seconds.  For  the  second  and  third  Mam¬ 
moth  Lakes  earthquakes  in  Table  1,  measured  X-P  times  .  4.9  and  2.5  seconds  would  be 
consistent  with  focal  depths  of  about  13  and  6.5  km,  respectively,  if  these  arrivals  were 
pP.  The  fourth  Mammoth  Lakes  event  has  a  complex  P-signature  and  the  timing  of  the 
second  phase  is  not  consistent  from  one  station  to  another. 

Magnitude  was  determined  from  computer  plots  of  the  digitized  P-waves  using  the 
standard  formula 

m4  =  logl0  +  B  (A). 

where  A  is  the  amplitude  of  vertical  ground  motion  in  nanometers  and  T  is  the  period 
corresponding  to  amplitude  A.  From  the  range  in  mk  values  given  in  the  Bulletin  of  the 
International  Seismological  Centre ,  it  appears  that  some  observatories  measured  the 
amplitude  of  the  initial  P-wave  for  the  Mammoth  Lakes  events,  while  others  measured 
the  amplitude  of  the  second  phase.  The  IAS  PEI  Commission  on  Practice  Concerning 
Amplitude  and  Period  Measurement  recommended  in  1979  that  "the  P  wave  amplitude 
measured  should  be  that  of  the  maximum  trace  deflection,  usually  within  the  first  25 
seconds  of  the  first  onset  or  before  the  arrival  of  the  next  clear  phase."  Other  workers 
have  recommended  measuring  A  and  T  within  the  first  few  cycles  (Zavadil,  1980)  or 
even  within  the  first  3/4  cycle  (the  "b"  amplitude)  of  the  P-wave  (Eisenhauer,  1980). 
Because  of  the  discrepency  between  maximum  amplitude  in  the  P-coda  and  amplitude  of 
the  P  onset,  we  determined  mt  for  the  Mammoth  Lakes  earthquakes  from  the  max¬ 
imum  amplitude  in  the  first  two  cycles  and  from  the  maximum  amplitude  in  the  first  12 
seconds  --  in  both  cases  taken  as  one-half  the  peak-to-peak  amplitude.  The  resulting 
individual  and  average  values  are  listed  in  Tables  6  and  7,  together  with  magnitudes 
determined  by  the  International  Seismological  Centre  (ISC)  and  Moscow  (MOS),  and 
deviations  6  m4  from  the  ISC  values. 

A  third  set  of  magnitudes  was  determined  from  measurements  of  A  and  T  on  traces 
produced  by  deconvolving  the  instrument  response  from  the  digitized  recordings.  The 
deconvolution  was  accomplished  using  a  program  written  by  W.  Peppin  to  calculate  the 
complex  transfer  function  of  a  seismograph,  given  the  period  and  damping  of  the 
seismometer  and  galvanometer,  together  with  the  maximum  magnification  of  the  system 
and  the  period  at  which  the  maximum  magnification  occurs.  The  records  were  tapered 
and  bandpass-filtered  (0.25-8.0  Hz)  before  deconvolution,  and  filtered  again  (0. 1-5.0  Hz 
bandpass)  after  deconvolution.  Figure  7  gives  examples  of  the  P-waves  for  the  first  large 
Mammoth  Lakes  earthquake  after  correction  for  instrument  response  and  filtering.  For 
these  records  m4  was  determined  only  using  the  maximum  amplitude  in  the  large  phase 
following  the  initial  P-wave.  Magnitudes  and  residuals  are  listed  in  Table  8. 


Table  6.  Magnitudes  and  Residuals,  Uncorrrected  Data,  First  Two  Cycles 


Average 
6  mh 


-.63  4.96 


IRK 

5.53 

-.57 

KZL* 

5.58 

-.52 

SEI 

5.67 

-.43 

SEM 

5.72 

-.38 

UST 

5.42 

-.68 

YAK 

6.16 

+.06 

-.63  I  4.97 


5.17 


.62  4.84 


5.34 

-.36 

5.78 

+  .08 

5.76 

+  .06 

-.49±.16 


+  .08±.0  I 


-.40±.31 


-.52±.0 


-.31±.14 


-.39±.40 


-.50±.39 


+  .1S±.44 


Table  7.  Magnitudes  and  Residuals,  Uncorrrected  Data,  First  12  Seconds 


05251633  05251944  05252035  05271450  Average 

Sta  m4  |  6  m*  m4  |  6  mi  m4  I  6  m4  m4  I  6  mb  6  mb 


+.16  I  5.46  +.16  5.66 


BOD 


ELT* 


IRK 


KZL* 


-.06  5.76 


SEI 

6.48 

+.38 

5.81 

+  21 

5.62 

+.32 

6.03 

+  .33 

+.31±.07 

SEM 

6.34 

+.24 

5.83 

+.03 

5.48 

+.18 

5.76 

+  .06 

+.13±.10 

lire™ 

6.18 

+.08 

5.52 

-.08 

5.22 

-.08 

5.71 

+.01 

-.02±.08 

YAK 

6.60 

+.50 

8.10 

+.50 

5.90 

+.60 

6.04 

+.34 

+  .49±.ll 

+.06±.12 


-,06±.0 


Table  8.  Magnitudes  and  Residuals,  Corrrected  Data,  First  12  Seconds 


Sta 


BOD 


ELT* 


IRK 


KZL* 


05252035 
ti4  I  6  mk 


05271450 
I  6  mh 


Average 
6  mb 


5.61 

+.01 

1  5.38 

i  +.08 

+  .46 

6.00 

+  .23 

5.80 

+  11 

5.64 

+  .80 

6  35 

5.86 

+.14±.0 
+  .38±.07 


+  .17±.07 
+  .02±.07 
+  .68±.08 


5.80 

5.6 

5.7 
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For  the  Tonga  earthquake,  values  of  mb  were  determined  from  the  maximum 
amplitude  in  the  P-wave,  from  traces  uncorrected  and  corrected  for  instrument  response. 
Four  of  the  stations  for  this  event  were  beyond  100*  and  distance  corrections  were 
taken  from  a  curve  developed  by  Ringdal  (1985). 


Table  9.  Magnitudes  and  Residuals  for  Tonga  Event 

Uncorrected  Data 

Corrected  Data 

Sta 

mb 

6  mb 

S  mb 

BOD 

5.85 

-.14 

6.12 

+  .12 

6.09 

+  .09 

6.22 

+.22 

IRK 

5.S8 

-.62 

5.42 

-.58 

6.19 

+.19 

6.17 

+.17 

SEM 

6.07 

+  .07 

6.06 

+.06 

UST 

5.79 

-.21 

5.78 

-.22 

YAK 

6.16 

+  .16 

6.15 

+.15 

Avg 

5.93 

5.99 

ISC 

6.0 

6.0 

MOS 

5.9 

5.9 

Discussion 

In  this  study  we  have  determined  magnitude  mb  for  five  earthquakes  on  25  and  27 
May  1980,  from  recordings  at  eight  Soviet  seismic  stations  on  a  4,300  km-long  profile 
from  eastern  Kazakh  to  eastern  Siberia.  In  general,  our  average  magnitudes  for  the 
Mammoth  Lakes  earthquakes  (Tables  7  and  8)  are  about  0.2  mb  unit  higher  than  those 
determined  by  the  ISC,  but  for  the  Tonga  event  our  m4  is  close  to  the  ISC  value.  As 
might  be  expected  our  values  are  also  about  the  same  as  magnitudes  attributed  to  Mos¬ 
cow  in  the  ISC  Bulletin.  For  the  data  uncorrected  for  instrument  response  our  average 
values  differ  from  the  Moscow  mb ’s  by  0.00±0.06  unit;  for  the  corrected  data  our  values 
are  higher  than  the  Moscow  mb ’s  by  0.06+0.06.  Magnitudes  determined  from  A  and  T 
measured  in  the  first  two  cycles  of  the  P-wave  for  the  Mammoth  Lakes  earthquakes 
(Table  6)  averaged  0.33  and  0.50  mk  unit  smaller  than  the  ISC  and  Moscow  values, 
respectively. 

The  largest  positive  mb  residuals  were  for  station  YAK  (average  6  mb  =  0.33  for 
the  uncorrected  data  in  Tables  7  and  9;  0.42  for  the  corrected  data  in  Tables  8  and  9). 
Yakutsk  is  located  on  the  central  Siberian  platform  in  an  area  where,  according  to 
Potap’ev  et  aL  (1974),  the  depth  to  crystalline  basement  is  shallow,  0  5-10  km.  About 
1C  km  west  of  Yakutsk  a  major  north-south  fault  offsets  basement  rocks,  with  the 
western  side  downdropped  by  about  4  km.  P-wave  velocities  in  the  crust  are  high,  6.2- 
6.4  km/sec  (Vol’vovsky,  1973),  the  crust  is  about  40  km  thick,  and  the  Pn  velocity  is 
about  8.0  km/sec  (Vol’vovsky  and  Vol’vovsky,  1975;  Ryall  et  al.,  1980). 

The  largest  negative  mk  residuals  were  observed  for  station  IRK  for  the  Tonga 
event  (-0.62  uncorrected,  -0.58  corrected),  for  which  the  raypath  crosses  the  southern 
part  of  the  Baikal  rift  zone  in  a  WNW  direction.  The  rift  zone  is  characterized  by  com¬ 
plex  geologic  structure,  including  deep  crustal  inhomogeneities  and  velocity  anomalies  in 
areas  of  recent  rifting.  Within  the  rift  zone  P-wave  velocity  in  the  crust  increases  from 
5.8  km  / sec  near  the  surface  to  6.4  km/sec  at  the  c.rust/mantle  interfare,  crustal  thick¬ 
ness  is  34-36  km,  and  the  Pn  velocity  is  7. 7-7. 8  km/sec.  (Puzyrev  et  al.,  1975;  Ryall  et  al., 
1980).  Puzyrev  et  al.  (1977)  attribute  the  low  Pn  ''elocity  to  partial  melting  of  upper 
mantle  material,  and  compare  the  Baikal  region  with  the  East  African  rift  system,  the 
Rhine  graben  and  the  Basin  and  Range  province.  It  is  interesting  that  IRK,  which  is 
located  on  the  southern  edge  of  the  central  Siberian  platform  just  north  of  the  Baikal 


rift,  has  smaller  values  of  6m *  (+0.03±0.10  uncorrected,  +0.03+0.07  corrected)  for  the 
Mammoth  Lakes  earthquakes,  for  which  the  P-waves  approach  the  station  in  a  SSW 
direction  across  the  platform  and  do  not  appear  to  be  affected  by  the  structure  of  the 
rift  zone.  According  to  Puzyrev  et  al.  (1977)  the  platform  is  characterized  by  uniform 
layering,  with  distinct  reflecting  horizons  at  depths  of  18-21  and  23-27  km,  crustal  thick¬ 
ness  of  about  40  km,  and  "normal"  Pn  velocity  of  about  8.1  km/sec. 

The  m*  residuals  for  station  SEM  are  intermediate  between  IRK  and  YAK.  SEM 
is  near  the  Irtysh  River,  on  the  boundary  between  the  west  Siberian  platform  to  the 
north  and  the  Kazakh  fold  system  to  the  south.  In  this  area  the  crustal  thickness  is 
about  45  km,  Pn  velocity  is  8.2-8.4  km/sec,  and  P-velocity  in  the  crust  is  high 
(Vol’vovsky  and  Vol’vovsky,  1975).  The  station  is  in  the  Zaysan  fold  belt,  which  is  rela¬ 
tively  simple  in  terms  of  stratigraphy  and  structure.  Sediments,  principally  of  Carboni¬ 
ferous  age,  contain  thick  limestone  deposits  and  lesser  amounts  of  interbedded  volcanics. 
Folding  and  faulting  are  relatively  minor  in  this  area,  compared  with  more  extensive 
folding,  faulting  and  intrusion  in  the  Chingiz-Tarbagatai  geanticlinal  zone  to  the 
southwest,  in  which  the  eastern  Kazakh  test  site  is  located  (Peyve  and  Mossakovsky, 
1982).  For  the  Mammoth  Lakes  earthquakes  £m*  for  this  station  is  +0.13+0.10  for  the 
uncorrected  data  (Table  7)  and  +0.17+0.07  for  the  corrected  data  (Table  8).  For  the 
Tonga  earthquake  it  is  -(-0.07  for  the  uncorrected,  +0.06  for  the  corrected  data.  As 
noted  above,  SEM  is  near  a  maximum  on  the  radiation  pattern  for  the  Mammoth  Lakes 
earthquakes,  and  near  the  null  axis  for  the  Tonga  event. 

Several  published  works  treat  magnitude  residuals  for  stations  of  the  Soviet  net¬ 
work,  and  provide  for  a  comparison  with  our  results  (Table  10).  First,  in  a  study  of 
magnitude  and  global  network  detection  capability  Ringdal  (1985)  recomputed  m*  for 
about  70,000  earthquakes,  using  A  and  T  values  given  in  the  ISC  Bulletin  and  a 
maximum-likelihood  estimation  technique  (Ringdal,  1976).  In  another  study,  North 
(1976)  calculated  mean  station  magnitude  bias  from  m*  values  given  in  the  ISC  Bulletin 
for  nearly  40,000  events  from  1964  to  1973.  The  biases  were  computed  for  the  "best"  (in 
terms  of  events  reported)  72  stations  with  respect  to  the  mean  magnitude  of  observa¬ 
tions  reported  by  this  set  of  stations,  with  the  requirement  that  an  event  be  reported  by 
more  than  15  of  these  stations  before  a  bias  was  calculated. 

In  work  by  Soviet  authors  Vanek  et  al.  (1978,  1980)  determined  magnitude  correc¬ 
tions  (Am* )  for  32  reference  stations  of  the  Unified  System  of  Seismic  Observations 
(ESSN)  of  the  Soviet  Union,  following  a  recommendation  of  the  KAPG  Conference  at 
Prague  in  1972  to  create  a  uniform  system  for  determining  magnitude  for  the  Eurasian 
continent.  The  station  corrections  were  determined  separately  for  phases  PV,  PVS  and 
PH  —  respectively  the  P-wave  recorded  on  vertical  mid-band,  vertical  short-period  and 
horizontal  mid-band  seismographs.  Calculations  were  also  made  separately  for  five  dif¬ 
ferent  source  regions  around  the  USSR  —  Alaska,  Japan,  thu  Phillipines,  Asia  and  the 
Mediterranean.  The  number  of  earthquakes  in  each  source  region  was  not  given  in  tiie 
1980  paper,  but  in  the  earlier  work  it  ranged  from  35  events  for  Asia  and  the  Mediter¬ 
ranean  to  93  events  for  Japan.  The  corrections  were  with  respect  to  one  of  the  reference 
stations,  OBN,  selected  at  least  in  part  for  its  bias  towa  d  large  mj  values.  Based  on  a 
comparison  of  Am*  values  for  the  various  source  regions  the  reference  stations  are 
grouped  according  to  the  number  of  corrections  needed  for  the  different  source  regions. 
Thus,  a  station  for  which  Am*  was  the  same  for  all  five  source  regions  was  classed  as  a 
reference  station  of  the  I  Kind,  a  station  with  the  same  correction  for  four  of  the  five 
source  regions  would  be  a  station  of  the  II  Kind,  etc.  Of  interest  to  our  study,  station 
SEM  has  the  same  PVS  correction  for  all  of  the  source  regions,  Am*  =  +0.32,  making  it 
a  station  of  the  /  Kind  and  indicating  that  it  has  a  6  m*  bias  of  -0.32  relative  to  station 
OBN.  Standard  deviations  are  not  given  by  Vanek  et  al.  (1980)  but  in  t,h<-  earlier  work 
they  average  ±0.05  m*  unit.  For  comparison  with  magnitude  bias  given  by  other 
authors,  we  reversed  the  sign  of  Am*  and  increased  all  of  the  resuiting  values  by  0.39  — 
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RingdaJ’s  (1985)  bias  for  station  OBN,  which  Vanek  et  aL  use  as  a  base  station. 


Table  10.  Comparison  of  Station  Residuals  and  6  mi  (OB2-NV) 


Ringdal  North  Vanek* 


*  -  Given  by  author*  as  ZWlj  corrections  relative  to  base  station  OBN.  Sign  reversed  and  all  values 
increased  by  0.39  for  comparison  with  Ringdal  (1985). 

1  ~  Mean  of  S  for  Tonga  earthquake  plus  average  S  for  four  Mammoth  Lakes  events  from  Table  7. 
Unconnected  data. 

2  ••  Mean  of  S  m*  for  Tonga  earthquake  plus  average  6  m^  for  four  Mammoth  Lakes  events  from  Table  8. 
Corrected  data. 

3  ••  Average  residual  (-0.10±0.35)  for  OB2-NV  with  respect  to  ISC  Bulletin  minus  average  station  residual. 


With  a  couple  of  exceptions  the  station  residuals  listed  in  Table  10  from  the  work 
of  Ringdal  (1985),  North  (1976)  and  Vanek  et  oL  (1978,  1980)  agree  to  within  a  few  hun¬ 
dredths  of  a  unit  of  mk .  Our  station  residuals  are  based  on  a  very  limited  data  set  and 
show  more  scatter  when  compared  to  those  of  the  other  authors.  However,  they  are  in 
general  agreement  with  the  published  values,  and  two  of  the  stations  —  SEM  and  YAK 
have  values  of  6  mk  that  are  in  excellent  agreement  with  those  of  Vanek  et  aL  and  Ring¬ 
dal,  respectively. 

As  a  final  step  we  computed  the  mt  bias  between  a  granite  site  at  the  Nevada  Test 
Site  and  the  Soviet  stations  listed  in  Table  10.  In  this  calculation  we  used  A  and  T 
measurements  listed  by  Der  et  aL  (1978)  for  83  seismic  events  for  the  period  September 
1976  to  March  1977,  recorded  by  a  digital  seismic  system  (SDCS)  at  station  OB2-NV 
(Climax  stock  at  north  end  of  Yucca  Valley  on  NTS).  In  the  Der  et  al.  report,  the  A 


values  are  peak-to-peak  amplitudes  in  nanometers,  and  magnitudes  are  computed  using 
distance  corrections  of  Veith  and  Clawson  (1972).  For  consistency  with  the  work  by 
Ringdal  (1985)  and  North  (1976),  magnitudes  were  recomputed  using  the  Gutenberg  and 
Richter  (1956)  corrections  and  dividing  the  peak-to-peak  amplitudes  by  two  to  obtain 
zero-to-peak  amplitudes.  Residuals  for  the  83  events  were  calculated  as 
6 mi  -  mi(0B2-NV)—mt(ISC),  and  these  were  averaged  to  obtain  an  average  Smt 
of  -0.09±0.39.  Five  values  that  fell  outside  the  2<r  limits  (-0.87  ^  ^  0.69)  were 

dropped  and  the  average  6mt  recalculated  to  obtain  -0.10±0.35.  This  number 
represents  the  average  bias  of  the  OB2-NV  site  with  respect  to  network-averaged  mi 
values  listed  in  the  ISC  Bulletin.  To  determine  the  bias  of  OB2-NV  with  respect  to  the 
Soviet  stations,  the  average  residuals  in  Table  10  for  those  stations  were  subtracted  from 
-0.10.  The  resulting  bias  values  are  listed  in  the  right-hand  column  of  Table  10. 

Of  particular  interest  to  questions  of  yield  verification,  the  mt  bias  of  the  NTS 
granite  site  with  respect  to  station  SEM  at  Semipalatinsk  is  -0.20,  with  a  range  of  -0.17 
to  -0.22.  The  smaller  of  these  figures  (-0.17)  is  based  on  the  study  by  Vanek  et  aL 
(1980),  and  the  larger  (-0.20,  -0.22)  are  from  our  measurements  of  Soviet  records  of  the 
Mammoth  Lakes  and  Tonga  earthquakes.  It  should  be  noted  that  the  bias  values  in 
Table  10  represent  only  the  bias  due  to  attenuation  in  the  upper  mantle  and  crust  under 
the  respective  seismic  stations;  they  do  not  include  other  effects  such  as  differences  in 
coupling  for  explosions  at  the  two  test  sites,  effects  due  to  tectonic  release,  or  those 
related  to  focusing  and  defocusing  of  seismic  waves  in  the  vicinity  of  a  given  explosion. 
The  reader  should  also  be  reminded  that  the  Semipalatinsk  station  is  about  100  km  from 
the  East  Kazakh  test  site,  and  that,  according  to  Peyve  and  Mossakovsky  (1982),  crustal 
rocks  under  the  test  site  have  been  subjected  to  greater  folding,  faulting  and  intrusion 
than  those  under  the  seismic  station. 
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Figure  Captions 


Figure  1.  Polar  projection  of  Eurasia  showing  stations  used  in  this  study.  Radius  of 
map  is  40  * . 

Figure  2.  Polar  projection  showing  location  of  Mammoth  Lakes  earthquakes  (solid  cir¬ 
cle)  and  stations  used  in  study  (triangles).  Radius  of  map  is  52  *. 

Figure  3.  Fault-plane  solutions  for  three  of  the  largest  Mammoth  Lakes  earthquakes 
(from  Given  ct  at,  1982).  Lower-hemisphere,  equal-area  projection.  Solid  circles 
—  location  of  raypath  to  Soviet  stations  used  in  study. 

Figure  4.  Polar  projection  showing  location  of  Tonga  earthquake  (solid  circle)  and  star 
tions  used  in  study  (triangles).  Radius  of  map  is  60  * . 

Figure  5.  Response  curves  for  SVK-M3  short-period  vertical  instruments  for  stations 
used  in  this  study. 

Figure  6.  P-waves  digitized  from  Soviet  recordings  for  the  Mammoth  Lakes  earthquake 
on  25  May  1980  at  16:33  GMT. 

Figure  7.  P-waves  corrected  for  instrument  response  for  the  25  May  1980,  16:33  GMT, 
earthquake. 
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CRUSTAL  STRUCTURE  IN  THE  VICINITY  OF  YUCCA  MOUNTAIN,  NEVADA 
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ABSTRACT 


We  have  interpreted  data  from  two  seismic  refraction  profiles 
recorded  by  the  (J.S.  Geological  Survey  in  the  vicinity  of  the  pro¬ 
posed  Nuclear  Waste  Repository  at  Yucca  Mountain  in  southern 
Nevada.  The  first  shot  consists  of  two  profiles,  the  longer  of  which 
has  33  recordings  of  a  nuclear  explosion  at  the  Nevada  Test  Site 
and  extends  from  45  to  110  km  and  the  other  has  26  recordings 
extending  50  to  90  km  from  the  shot  point.  The  second  consists  of 
71  recordings  of  a  chemical  explosion  near  Beatty,  Nevada,  and 
extends  «70  km.  east-west  across  Yucca  Mountain.  Gravity  and 
well  log  data  are  used  to  correct  the  refraction  travel-time  data  for 
delays  due  to  lateral  variations  in  near-surface  structure.  The 
travel-times  for  the  two  profiles  were  combined,  reparameterized 
in  terms  of  delay  time  and  ray  parameter,  and  inverted  for  aver¬ 
age,  extermal  velocity-depth  bounds  for  the  region.  These  show  a 
steep  gradient  in  the  shallow  crust  (depth  less  than  5  km)  with 
speeds  increasing  from  «4.5ks-1  to  6.2  ks_1  and  a  shallow  gradient 
below  5  km  where  speeds  increase  from  «6.2ks-1  to  6.5  ks-1  at  25 
km.  A  specific  velocity-depth  model  consistent  with  the  extremal 
bounds  has  been  found  by  two-dimensional  ray  tracing.  The  struc¬ 
ture  of  the  upper  crust  is  highly  variable  with  the  Paleozoic  base¬ 
ment  varying  between  the  surface  and  3.5  km  depth.  A  thick,  low- 
speed  section  is  indicated  by  large  delays  across  the  Crater  Flat- 
Prospector  Pass  caldera  complex.  A  prominent  arrival  1.5  to  2.0 
seconds  after  the  initial  P-wave  time,  observed  in  the  distance 
range  48  to  65  km,  is  interpreted  as  a  mid-crustal  boundary  at  13 
to  17  km  depth.  The  profiles  are  of  insufficient  length  to  see  Pn, 
however  a  prominent  arrival  3  to  4  seconds  after  Pg,  observed  at 
distances  greater  than  80  km,  is  interpreted  as  PmP  from  the  Moho 
at  approximately  35  km  depth. 


Introduction 

Since  1978  the  U.S.  Geological  Survey  has  conducted  a  series  of  seismic 
refraction  experiments  to  elucidate  the  shallow  crustal  structure  in  the  vicinity 
of  Yucca  Mountain,  on  the  western  boundary  of  the  Nevada  Test  Site  (NTS) 
(Hoffman  and  Mooney,  1984).  The  purpose  of  these  experiments  is  to  aid  in  site 
evaluation  for  the  proposed  nuclear  waste  repository  at  Yucca  Mountain.  How¬ 
ever,  much  of  the  data  collected  is  suitable  for  delineating  deeper  crustal  struc¬ 
ture  in  the  vicinity  of  NTS.  We  have  obtained  the  digital  data  from  the  U.S. 
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Geological  Survey  and  have  embarked  on  a  program  of  interpretation  of  the  data 
in  terms  of  the  deeper  crustal  structure.  Figure  1  shows  the  location  of  the 
refraction  profiles.  To  date,  our  analysis  have  been  confined  to  two  profiles, 
those  from  shot  point  3  ,  a  nuclear  explosion  on  Pahute  Mesa,  and  shot  point  4,  a 
chemical  explosion  near  Beatty,  Nevada.  Data  from  shot  point  3  was  collected 
along  two  profiles.  The  first  profile  (3A)  is  along  the  axis  of  Yucca  Mountain,  in 
the  distance  range  of  40  to  110  km,  and  the  second  profile  (3B)  is  across  the 
Crater  Flats  Caldera  Complex  with  a  distance  range  of  50  to  90  km.  Data  from 
shot  point  4  was  collected  along  a  70  km,  east  trending  profile  from  near  Beatty, 
across  Bare  Mountain.  Yucca  Mountain,  to  Skull  Mountain  on  the  NTS.  All  the 
seismograms  are  vertical  component  velocity  records,  digitized  at  200  samples 
per  second. 

Data  Analysis 

Figure  2  is  record  sections  of  the  data  recorded  along  the  profiles  3A  and 
3B.  These  record  sections  have  been  normalized.  Delays  of  the  order  0.3  to  0.5 
seconds  are  observed  in  the  profile  3A  at  distances  of  approximately  45  to  60 
km.  A  reflective  phase  is  identified  at  distances  between  80  to  110  km  and  is 
interpretated  as  reflections  from  the  Moho  (PmP).  The  record  section  for  profile 
3B  shows  delays  of  0.4  to  0.5  seconds  in  the  distance  range  of  52  to  60  km. 
These  can  be  attributed  to  the  large  amount  of  volcanic  tuffs  and  sediments  in 
the  caldera  complex,  which  erupted  approximately  10  m.y.  ago  (Carr.  1984). 
Figure  3  is  a  record  section  of  the  data  collected  along  the  profile  from  shot 
point  4.  At  distances  between  48  to  65  km,  the  first  arrivals  are  attenuated  and 
emergent,  however,  there  is  a  strong  reflected  phase  1.5  to  2.0  seconds  after  the 
first  arrival  coming  in  at  these  distances.  This  phase  is  identified  as  reflections 
from  an  intra-crustal  boundary.  The  length  of  this  profile  is  insufficiently  long  to 
observe  any  P„  arrivals.  The  average  error  in  picking  arrival  times  for  profiles 
3A  and  3B  are,  0.01  seconds,  and  0.08  seconds,  with  maximum  errors  of  0.02  and 
0.18  seconds  at  distances  of  70  km  and  109  km,  respectively.  Profile  4  has  an 
average  timing  error  of  0.03  seconds  with  maximum  of  0.11  seconds  at  47  km. 

The  travel-time  data  for  these  profiles  have  been  corrected  for  the  near 
surface  delays  using  the  delay  times  from  the  shallow  structure  interpretation 
of  Hoffman  and  Mooney  (1984),  and  inferences  of  the  near  surface  structure 
from  gravity  data  (Healey  et  al.,  1978).  The  two  data  sets  from  shot  points  3  and 
4  were  combined,  and  reparameterized  in  terms  of  the  delay  time,  t  and  the  ray 
parameter,  p.  There  are  several  advantages  of  reparameterizing  the  travel-time 
data  in  this  manner.  One  is  that  the  r(p )  function  and  the  earth  structure  arc 
linearly  related,  allowing  the  application  of  linear-inverse  theory  to  invert  the 


-3- 


r(p)  data  in  order  to  construct  and  evaluate  consistent  earth  models.  Another 
advantage  is  that  the  r(p)  function,  which  is  a  single  value,  monotonically 
decreasing  function,  will  show  a  discontinuous  in  the  presence  of  a  low-velocity 
zone.  The  r(p)  data,  along  with  time-distance,  cross-over  and  triplication  con¬ 
straints,  were  inverted  for  bounds  on  the  velocity  structure  using  the  pro¬ 
cedures  outlined  by  Garmany  et  al. ,  1979,  Orcutt  et  al. ,  1980,  and  Orcutt,  1980. 
The  latter  constraints  tend  to  smooth  the  inversion  of  velocity-depth  function. 
The  resulting  bounds  on  the  velocity  depth  function  are  shown  in  Figure  4. 
These  bounds  indicate  a  steep  gradient  in  the  shallow  crust  for  depths  less  than 
5  km.  The  speed  corresponding  to  this  depth  is  approximately  4.5  ks-1  to  6.2 
ks-1.  At  deeper  depths,  the  velocity  gradient  shallows  to  speeds  which  increase 
from  6.2  ks-1  to  6.5  ks-1  at  25  km  deep.  The  velocity-depth  bounds  of  Figure  4 
contain  all  the  possible  models  compatible  with  the  travel-time  data,  but  do  not 
in  themselves  represent  the  actual  velocity-depth  function.  To  develope  a 
specific  velocity-depth  model,  we  used  an  average  velocity-depth  function  as  a 
starting  model,  and  a  two-dimensional  ray  tracing  routine.  (McMechan  and  Moo¬ 
ney,  1980).  By  fitting  the  obse-'red  travel-time  data  to  the  calculated  travel- 
times  based  on  the  specific  velocity-depth  model,  we  have  constructed  a  prelim¬ 
inary  model  which  is  consistent  with  the  velocity-depth  bounds  (Figure  5). 

We  have  derived  a  preliminary  earth  model  using  the  delay  time  method 
which  is  consistent  with  the  Paleozoic  basement  that  varies  from  surface 
outcrops  to  depths  of  3.5  km.  Large  delays  observed  over  the  Crater  Flats  Cal¬ 
dera  Complex  is  explained  by  a  3  to  4  km  thick,  low-speed  section,  underlying 
the  complex.  We  have  interpreted  the  reflected  phase  in  seen  in  profile  4  as  an 
intra-crustal  boundary  whose  depth  ranges  from  13  to  17  km.  From  the  profiles 
3A  and  3B,  we  estimate  the  depth  to  the  Moho  to  be  approximately  35  km.  Our 
model  is  further  constrained  using  gravity  data.  Figure  6  shows  the  results  of 
gravity  modeling  with  our  model.  Our  analysis  suggests  the  following  density 
versus  speed  relationship:  a  density  of  2.7  g/cc  corresponds  to  speeds  of  6.1 
ks-1  and  higher,  a  density  of  2.5  g/cc  corresponds  to  speeds  of  4.0  ks-t  to  5.7 
ks-1,  a  density  of  2.45  g/cc  corresponds  to  speeds  of  3  ks-1  to  4.5  ks-1,  and  a 
density  of  2.1  g/cc  corresponds  to  speeds  of  2.0  ks-1  to  3.1  k.,  We  have  also 
incorporated  well  log  data  in  our  gravity  modeling  which  constrains  the  surface 
layer  at  Yucca  Mountain  (Spengler  and  Rosenbaum,  1980). 

Discussion  and  Conclusion 

Our  preliminary  model  correlates  well  with  that  obtained  by  Hoffman  and 
Mooney  (1984),  derived  by  trial  and  error  fitting  of  the  observed  travel  times. 
Our  data  indicate  that  the  deeper  layers  are  not  very  well  constrained,  however. 
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We  combined  data  from  profiles  3At  3B,  and  4  for  our  inversion.  Combination  of 
these  data  results  in  relatively  high  scattering  which  affects  the  r  analysis. 
Better  resolution  of  the  near  surface  structure  would  allow  us  to  reduce  this 
scatter.  We  have  relatively  good  control  for  the  lateral  structure  in  the  east- 
west  line,  however  most  of  the  scattering  effects  can  be  attributed  to  the  near¬ 
surface  lateral  variations  between  shot  point  3  and  its  closest  station,  45  km 
away.  We  have  no  control  on  the  shallow  structures  at  this  distance.  The  shal¬ 
lowest  structure  we  can  define  is  approximately  5  to  6  km  deep,  based  on  the 
velocity  structure  of  profile  4.  Since  the  line  passes  through  Timber  Mountain 
caldera,  we  estimate  that  there  exist  a  thick,  low-speed  section  in  Timber  Moun¬ 
tain. 

The  high  speeds  seen  in  the  Crater  Flats  area  suggest  that  the  caldera  con¬ 
tains  a  high  density  fill.  This  is  not  implausable,  considering  that  highly  welded 
tuffs  are  known  to  exist  there. 
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Figure  Captions 

Figure  1.  Map  showing  the  locations  of  Shot  points  3  and  4  (closed  circle),  and 
station  locations  (open  triangles). 

Figure  2.  Record  section  for  profiles  3A  and  3B. 

Figure  3.  Record  section  of  profile  4 

Figure  4.  Extremal  velocity-depth  bounds  for  profiles  3A,  3B,  and  4. 

Figure  5.  Two-dimensional  ray  tracing  of  the  specific  velocity-depth  model  that 
is  consistent  with  the  velocity  depth  bounds  in  Figure  4. 

Figure  6.  Comparison  of  observed  and  calculated  gravity  data  for  the  model  in 
Figure  5. 
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